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B vuroHe 1959 roma Hay4Ho-KOHCYIbTaTMBHbIM KOMU TeT Mex yyHapOyAHOrO areHTCTBa MO aTOMHOM SHEPrMu pe- 
KOMeHOBan ATeHTCTBY MPMHATb Mepbl, cnocoOcTByW ume OOmMeHy MHdopMauMen no uccreqOBaHMAM B OONacTU 
cbu3uKM M1a3MbI “ KOHTpONMpyemoro cuHTe3a. B KavecTBe KOHKpeTHOrO mpezyomReHNA KOMMTeT PCKOMCHOBalI 
BbINycK MexXyHapoyzHOrTO Hay4Horo wKypHata B STON OOnacTu. B cooTBeTCTBMM c 3TMM NpezNomeHNem AreHT- 
CTBO B HacTOALIee BPCMA IIPMcTynNaeT K M3TZaHMIO MyPHaNa «APepHbIM CUHTE3>». 

MbI cuutTaeM, 4TO STOT HOBbIM KypHaN SGyzeT cOMeUCTBOBATL YIY4UWIEHMIO COTPYAHMYNeCTBA Mey cTpaHaMu 
B PeWUICHUM TPYAHOM 3aqa4uu HAaxORKAeHMA TyTeu NONyueHuA 9HEPrMM M3 KOHTPONMpyemoro AepHOTO CUHTe3A. 
Jlaxe ecim cepbe3HbIe MU COriacOBaHHble yCUJIMA, HAMpaBNeHHbIe Ha pellieHve STON 3azauM, cpa3y uM He NOZ- 
BewyT K MpaKTM4ecKMM MCTOYUHMKAaM SHEPIMM, 9TA MOMbITKa, MOMMMO yruyONeHMA OCHOBHBIX HayUHbIxX 3HaHMit, 
3Ha4MTeIbHO PacCWIMPMT NO3HAHMA B OONacTM GbM3MKM MNa3MbIl. DTO yceumMe HeECOMHEHHO MpuBezeT Take K 
HOBbIM M BaxKHbIM TEXHNUYECKMM U3Z00PeTCHUAM MMPHOTO xapaktTepa. 

HeoOsiKHOBeHHbIM CbakT 3aKTIOUAeTCA B TOM, 4TO, XOTAH KOCMOC — KaK OH MOHMMAaeTCA B HACTOALIee BPeMA 
— MOUTM NOJHOCTbIO COCTOMT U3 MWia3Mbl, YYCHbIeC JIMUWIb CpaBHUTeNbHO HeEAaBHO AOcTursM mporpecca B MOHM- 
MaHMM CBOMCTB 9TOTO «4YeTBEPTOTO COCTOAHMA MaTepuM». Kak moKa3bIBaloT ONyOJMKOBaHHbIe B MKypHase «fApep- 
HbIM CMUHTe3» CTAaTbM, Mia3Ma OONaqaeT OTHIOQb HE MpOCTbIMM CBOMCTBAMM, UM HeNerKO MOAMAWTCA YTM CBOIi- 
cTBa UCCIe_OBaAHUW. 

Bce uutTaTeiu *ypHasa «AepHbIM CMHTe3» MpursiawialwTcCA yuacTBOBaTb B ero u3qaHuM. Ux 3ame4aHna, 
IpeqJIO#KeHMA, a TAK*KEe cTaTbu OyzyT MpMBeTcTBOBaTbcA. OT uMeHwu MexmyyHapogzHoro areHcTBa MO aTOMHOM 
9Heprmuu A xOTe ObI BbIPASMTb MOKO MUCKPeHHOW GJlaroyqapHOCcTb BCeM TeM, KTO MOMOr HaM B CO3ZaHMM 9TOTO 
*KypHadga. 


En junio de 1959, el Comité Consultivo Cientifico del Organismo Internacional de Energia Atomica recomend6 
que el OIEA fomentase la difusi6n y el intercambio de informacién en el terreno de las investigaciones sobre la 
fusi6n nuclear controlada y la fisica del plasma. Ademas, el Comité recomend6 concretamente que se iniciase la 
publicacion de una revista cientifica internacional que tratase de estas cuestiones. Como consecuencia de esta 
propuesta, el Organismo lanza hoy el primer numero de “FUSION NUCLEAR”. 

Creemos que esta nueva revista contribuira a acentuar la cooperacion entre las diversas naciones en la dificil labor 
de descubrir procedimientos para la obtencién de energia partiendo de la fusi6n nuclear controlada. Aun cuando los 
mas serios y armonizados esfuerzos que se dediquen a esta labor no conduzcan inmediatamente a disponer de fuentes 
energéticas, la empresa extendera muy considerablemente el conocimiento de la fisica del plasma ademas de ampliar 
el caudal de conocimientos cientificos basicos. Estos esfuerzos se traduciran sin duda alguna en la invencién de 
nuevos e importantes dispositivos utilizables en tiempo de paz. 

Es realmente sorprendente el hecho de que, aunque el universo — segun el concepto que hoy se tiene del mismo — 
esta formado casi enteramente por plasma, hasta hace relativamente poco no hayan logrado los hombres de ciencia 
progresar en el conocimiento de las propiedades de este “‘cuarto estado de la materia”. Como se refleja en las paginas 
de “FUSION NUCLEAR”, la investigacion de las propiedades del plasma no es labor facil ni sencilla. 

Quienquiera que lea y haga uso de “FUSION NUCLEAR” queda invitado a colaborar en la revista. Sus observa- 
ciones y sugerencias, asi como los originales que envie, seran bien acogidos. En cuanto a quienes nos ayudaron a 
lanzar esta revista, deseo expresarles aqui, en nombre del Organismo Internacional de Energia Atomica, mi sincero 
agradecimiento. 
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PLASMA OSCILLATIONS (I)* 
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This article is the first part of a review of the theory of wave phenomena in plasmas. The basic kinetic theory is 
developed taking into account coulomb collisions. From this magnetohydrodynamic equations are derived. It is 
demonstrated that it is legitimate to deal with a closed subset of these equations either in the limit of collision- 
dominated phenomena, or in the limit where the effective phase velocities of the phenomena of interest are much 
greater than thermal speeds. In case of collision-dominated plasmas the theory of transport coefficients is discussed. 

These equations are then applied to an extensive treatment of small amplitude wave phenomena in plasmas. 
A discussion of the dissipative effects on hydromagnetic waves is given. Hydromagnetic waves are also considered 
from the Chew, Goldberger and Low theory. Longitudinal and transverse oscillations in current-carrying plasmas 
are also discussed. 

Oscillations of a cylindrical plasma are considered and the phenomenon of ion cyclotron resonance is discussed. 
The possibility of radiation by plasma oscillations by a uniform sphere is exhibited. Some general results on the 
stability of longitudinal electron oscillations in non-uniform plasmas are given. 

A brief treatment of large amplitude electron oscillations is given and the breaking of these oscillations as a 
dissipative mechanism for the organized plasma motion is discussed. 

Part II of this paper, to appear in a future issue of this journal, will be devoted to the discussion of piasma 
oscillations directly from the kinetic theory. 
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I. Introduction 


In this article we shall attempt to present a systematic 
survey of the theoretical knowledge of wave phenomena 
in plasmas. We shall use the term “plasma” to describe 
a fully ionized gas composed of positive and negative 
ions of effectively zero net charge density. For the sake 
of convenience we limit our attention to a plasma com- 
posed of only two species (positive ions and electrons); 
the extension of the results to gases consisting of more 
than two species is straightforward and will not be con- 
sidered here. Some of the results are also applicable to 
conducting liquids and gases not fully ionized. 

Astrophysical examples of such systems are the ion- 
osphere, the solar corona, stellar interiors and the at- 
mosphere of hot stars, the H,,; regions of interstellar space, 
and cosmic radio sources. Plasmas are also of fundamental 
interest in the proposed schemes for controlled thermo- 
nuclear fusion, plasma propulsion of space ships, micro- 
wave generators, and gas discharge devices. 

One distinctive feature of plasmas is the special nature 
of collisions, in contrast to the two-body collisions 
familiar in the conventional theory of neutral gases. 
Because of the long range nature of the Coulomb forces 
which govern the encounters of charged particles, co- 
operative phenomena are of major importance in describ- 
ing the dynamics of a plasma. Also these long range col- 
lisions lead to a diffusion in velocity space and are most 
appropriately described by the Fokker-Planck equation. 
The elastic or inelastic two-body encounters of ions or 
electrons with neutrals or perhaps with the core of an 
incompletely ionized atom are usually negligible in 
plasmas and will be considered, effectively, only when 
dealing with dissipative effects in collision-dominated 
plasmas. 

There are two levels of description of plasmas which 
have received major consideration during the past dec- 
ade. The first of these is based on the direct solution of 
the Boltzmann transport equation —usually referred to 
as the kinetic theory description. This method has the 
intrinsic merit of providing a more complete description 
of the dynamics of a plasma. However, this procedure is 
mathematically more difficult to handle. The second 
method is based on using a closed set of moment equa- 
tions to characterize the behaviour of a plasma. This de- 
scription, usually referred to as the hydrodynamic or 
hydromagnetic description, though of limited applica- 
bility (in a sense to be discussed in detail later) has the 
merit of relative mathematical simplicity and is quite 
helpful in the physical understanding of the various 
phenomena involved. This method is applicable either 
when (1) collisions dominate so that one has local thermo- 
dynamic equilibrium and can apply the Chapman-Enskog 
expansion technique, or (2) the phase velocity of the 
wave is much larger than the characteristic thermal 
speed in a collision-free plasma. The latter will be re- 
ferred to as the “Low Temperature Approximation.” 

For plasmas of interest under astrophysical and lab- 
oratory conditions quantum mechanical effects are un- 
important in considering most of the dynamics of plasmas. 

This work, therefore, divides itself naturally into two 
parts. In part I the kinetic theory description is formulated 
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and its relation to the hydrodynamic description est- 
ablished. The resulting hydromagnetic equations are then 
applied to a variety of wave phenomena. Part II of this 
article, to appear in a later issue of this journal, will be 
devoted to the discussion of wave phenomena based 
directly on the kinetic theory. 

The theory described here is concerned mainly with 
the small motions about a steady state of the plasma. The 
general case of large motions is considered only briefly 
in Section V. 

Section III deals with oscillations in unbounded 
homogeneous plasmas, admittedly a severe idealization. 
However, these treatments do serve to give some in- 
sight into the basic physical processes and the mathe- 
matical methods. The neglect of boundaries would imply 
that the description is valid only for wavelengths much 
smaller than the characteristic length of the system and 
for time intervals short enough that the disturbance has 
not reached the boundaries and been reflected and trans- 
mitted. 

Unfortunately, because of the inherent difficulty of 
establishing in the laboratory plasmas uncomplicated by 
boundary effects, impurities, etc., there has not been, as 
yet, much experimental verification of the theoretical 
investigations. 


II. The basic equations 


By and large one is interested only in the knowledge 
of some macroscopic properties of the plasma like the 
mean particle density, mean velocity, etc., rather than 
in the detail afforded by the description of the individual 
motion of the particles. It is, therefore, appropriate to 
characterize a given class of particles by a distribution 
function /(x,v,/) such that fd*x d*v represents at time / 
the probable number of particles with velocities between 
vand v + dv, and with positions between x and x +-dx. 
It is possible to relate f to the more recondite density in 
the phase space of all the particles, but a discussion of this 
matter is beyond the scope of this article; suffice it to say 
that the introduction of / is compatible with the follow- 
ing notions about particle interactions in plasmas. 

The particles, being charged and possessing no in- 
ternal degrees of freedom at the range of energies of 
interest here, interact primarily through the electro- 
magnetic fields. In particular, when two particles are 
close together it is their Coulomb interaction which 
predominates. Consider then such two-body Coulomb 
encounters. They can be characterized, in their center 
of mass system, by the impact parameter and the relative 
velocity. The collisions can then be conveniently divided 
into three classes. 

The first of these corresponds to impact parameters 
less than the mean impact parameter for 90° deflection, 
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where 7, and 7, are the charge numbers of particles in- 
volved, m, and m, are their masses, <v*> represents their 
mean square relative velocity, and e¢ is the electronic 
charge. It is usually legitimate to describe such en- 
counters by the familiar two-body collision integrals of 
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the kinetic theory of gases when writing the kinetic 
equation for the distribution function. In most cases of 
interest these collisions are negligible. 

The second class comprises the many-body encounters 
which can be regarded as composed of a succession of 
uncorrelated small angle Coulomb deflections. These 
correspond to elementary Coulomb scattering events 
characterized by impact parameters greater than the dis- 
tance of closest approach b, but less than the electron 
Debye length 


. (1/3) w <v*>7] 3 . 
4a Ne | (2) 
where m is the electron mass, <»”) is the electron mean 
square velocity, and N is the electron number density. 
These encounters are most appropriately described by 
the Fokker-Planck form for the collision integral in the 
Boltzmann equation. 
When there are many particles in a sphere of radius 
4, (i.€., 4,.V > 1), density fluctuations over distances 
4,, ate suppressed by macroscopic electric for- 
ces. Thus the third class of uncorrelated two-body en- 
counters, those with impact parameters greater than /,,, 
do not occur per se.These are replaced in the kinetic 
eqation by the correlated effect of many particles through 
the macroscopic electromagnetic fields described by the 
Maxwell equations; this leads to co-operative phenomena. 


1. BoLTZMANN EQuaATION 
The distribution function for each species of particles 
is governed by the Boltzmann equation 
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where V and V,, represent the gradient operators in co- 
ordinate (x) and velocity (v) space respectively, F de- 
notes the force per unit mass, and (0//0/)con. represents 
the time rate of change of / resulting from collisions. It 
may be remarked here that there is a separate Boltzmann 
equation for each species of particle. We now discuss 
the Fokker-Planck form for the collision term. 


1.1 Fokker-Planck Form of the Boltzmann Equation 


It has been demonstrated by Jeans that when pazsticles 
interact through long range forces, the cumulative effect 
of weak deflections resulting from the relatively distant 
encounters is more important than the effect of occasional 
large deflections. The cumulative effect of weak de- 
flections is most appropriately taken into account by 
writing the collision term in the Fokker-Planck form. 
This reduction is possible only if there exist time inter- 
vals A¢ long enough for the particles to suffer a large 
number of collisions but small enough for the net mean 
square change in velocity, <|Av|*>, to be small com- 
pared with the mean square velocity, and to be of 
order A/. 

Let P(v; Av) denote the probability that in the time 
interval A+ a particle with velocity co-ordinate v under- 
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goes a displacement Av. 
then given by 


The distribution function is 


f(x, v, 7) 


\f (x, Vv Iv, 7 I*)P(v Iv; Av) d?Av. (4) 


We may remark here that in expecting this integral 
equation to be true we are actually supposing that the 
course which a particle will take depends only on the 
instantaneous value of its physical parameters and is 
entirely independent of its past history. In general prob 
ability theory a stochastic process which has this char- 
acteristic, namely that what happens at a given instant 
of time ¢ depends only on the state of the system at 
time /, is said to be a Markoff process. 

Since A¢ and 
Taylor’s expansion we can write (4) in the form 


lv are both assumed small, using 
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Eq. (5) gives for the time rate of change of / resulting 
from the cumulative effect of small deflections: 
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2 dv 
where we have retained terms only to lowest significant 
order and 


Iv \ f(x, v, 4 P (wv; Av) Av Av, (8) 


Iv lv \f(x,v,A P(v; Av) AvAvdAv. (9) 
Eg. (7) represents the collision term in the Fokker- 
Planck form. 

Consider the two-body encounters of particles of type / 
with those of type / and let o,,; (9,v) denote the differen- 
tial cross section for scattering of their relative velocity 
vector v = v, —v, through an angle §. Then the as- 
sociated momentum transfer cross section is given by 

OQ; (v) = 22 \6;;,(0,7) (0 cos 9) sinOd@. (10) 
In terms of O, retaining only the dominant terms in the 
mean square scattering angle, one can write for the con- 
tribution of such encounters to stochastic 
averages, when polarization effects are negligible, 
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1.2 Fokker-Planck, Equation for Coulomb Force 


For the case of Coulomb collisions, the differential 
scattering cross section is given by 


Z, 2,0? (m,-+- m;))? 1 
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(13) 
The integral occurring in O will diverge at the lower 
limit of integration (0 = 0). However, the small angle 
deflections correspond to scatterings with very large 
impact parameters, and the divergence arises from the 
long range nature of the Coulomb force. This divergence 
is eliminated when we take into account the shielding 
that arises from the polarization of charge surrounding 
the scatterer. The polarization screens the scattering 
particle and provides a natural cutoff on the maximum 
impact parameter of the order of the Debye length 4,. 
We thus take for the minimum scattering angle 


(m,+- m;) e 
4m; m; <v *» An (14) 


Using Eqs. (13) and (14) into Eq. (10) we obtain 
2,2, 0 (m+ “I 
n 


Mm; M2; v" 


2m, m; v* 


Mm; m;) e* 
OC = 7 ( a J . —, 15 
Qi=4 2m, m;<v*> hp (15) 
Using the foregoing results in the Fokker-Planck equa- 
tion we obtain after an integration by parts what is called 


the Landau form for the collision term: 


QE) =, (REA 


m;* m* 


(m,; + mj) e 
‘ail ; 2m, m; <v®> Ap 
’ 9 
vI—vv sal 
Vv 3 yy ™ 
x | da? 9; 3 (M4500 Si — SiN o Sif: (16) 

e a } 

It may be noted that, apart from the slowly varying 
logarithmic term, in order of magnitude the collision 
terms are inversely proportional to 9) 55 the mean 
square relative kinetic energy. Thus in the limit where 


this factor is appropriately large, collisions are negligible. 


1.3 The Moment Equations 


The particle density and the local mean velocity are 

defined as 

N = \f (x, v, 4 dv (17) 
and 

if. ' 

u yy \ £66 v4) vay. (18) 
In terms of these moments we can define some of the 
other relevant macroscopic quantities *: 
o=—>mN 


— 


mass density, 


(19) 
1 ’ : ; : 
Uy 2 2 m Nu = velocity of thecenter of mass, (20) 


* We shall adopt the convention that whenever m appears after a 
summation sign, it is the mass of any of the species of the plasma; 
otherwise m# denotes the electron mass and M, the proton mass. 
The electron charge is taken to be —e while that of the proton, + e. 
Similarly, w,. will denote the electron or ion cyclotron frequency 
whereas Ww, and @, will denote the electron and ion cyclotron fre- 
quency respectively. 


6 


€= > e N=charge density, 


(21) 


J p e Nu = current density, (22 
where the summation is to be carried over all the species 
of the system. The random velocity with respect to the 
center of mass is defined to be 

w =v— U,. 


(23) 
It is convenient to define the other relevant quantities 
with respect to the center of mass velocity. Thus 


U 7 | wf (x, v,/) d8v mean velocity of (24) 

ied each of the species 
with respect to the 
center of mass, 


- pressure tensor, 


(25) 


P—w» \ wwf (x, v,/) dv 
and 


Q=m | w w w/ (x, v, 4) dv = heat flow tensor, (26) 


where Q is a tensor of rank 3 completely symmetric in 
all its indices. We further define 


P—»P and Q=)DQ. (27) 
It is clear that 
> mNU=0. (28) 
Note that 
| “1 
> Trace P \ a” nw? f(x, v,/) dv (29) 


represents the mean kinetic energy relative to the center 
of mass and is referred to as the internal energy, or the 
thermal energy of the gas. This suggests that we can 
define the generalized temperature by the relation 


3 
: Trace P= 5 NO, 


(30) 


where 90 — AT, & is the Boltzmann constant, and TJ is 
the temperature. We shall always use @ to denote the 
temperature in energy units. 

This definition of temperature is natural when colli- 
sions dominate, in which case the various particle dis- 
tribution functions are very closely locally Maxwellian, 
relative to the center of mass velocity, corresponding to 
the same temperature. Occasionally one encounters an 
alternative definition of temperature, where it is defined 
separately for each species of the plasma, as the mean 
kinetic energy of that species relative to its own mean 
mass velocity. The other moments can be similarly de- 
fined. This latter definition is particularly convenient 
when collisions are negligible and will be used in part II 
of this review. The two notions, of course, coincide if 
we consider the motion of a single species of the plasma, 
a picture which is reasonable in investigating high fre- 
quency waves where only theelectron motionis significant. 
In general, however, the two definitions are quite dif- 
ferent and care must be exercised to make sure which is 
being employed in any given context. 
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The vector 
“I 
\ 5m ww f (x, v,¢) dv (31) 


] 
q 71:0 


represents the flow of internal energy relative to the 
center of mass and will be called the heat flow vector. 

Since x and v are independent vectors, it is con- 
venient to write the Boltzmann equation (3) as 


‘a of : of 
Mf=a,tV M+ FN=(3,) 62) 
Quite generally, we can write 
e ] 
F = . (x, 7) + =o? B (x, d| —Vq, (33) 


where E and B denote the electric and magnetic fields 
respectively, and @ is the gravitational potential. The 
electromagnetic fields, of course, satisfy Maxwell’s 
equations: 


42 10E i 10B 
Vv xB : J+ > VxE > a (34) 
V-E=42¢e, and V-B=0; (35) 


Gaussian units will be used throughout in this work. 
The gravitational potential satisfies the Poisson equation 


V2¢ —42Go, (36) 


where G denotes the gravitational constant. 

In order to derive the equations governing the’ macro- 
scopic quantities, it is convenient to write the Boltzmann 
equation in terms of the random velocity w. This can 
be readily effected by the transformations 


f(x, v, 4) > f (x, w, 4), (37) 

:.s. 

we ee (38) 

vV-V— (Vv Uy): Foie (39) 
V.>Vo- (40) 


The Boltzmann equation now takes the form 


of , , 
2/=F+V-[wt wos |+%-b=(F) an 


where 


e 3 
b = [E (x, 7) + - U, » B(x.) 


U a - 
—Vgo— _ -(w-V)U,+ awe B (x, /), (42) 
and 
d rs) 
ae t eV (43) 


denotes the time derivative as we follow the motion of 
an element of the fluid. 

Multiplying Eq. (41) by 1, ww, and mww, and as- 
suming that f vanishes sufficiently strongly in velocity 
space so that all the surface integrals vanish, we obtain 
for the moment equations 


oN 


N+. [vU+NU9]| =0, (44) 


PLASMA OSCILLATIONS (1) 


dU, 


miN ry 


-V: \P } mNUV| 


MNU+NelE +. U, x BI 
mNVp—mN(U-V)Uy+UxB, (45) 

and 

°.P+V-(Q+U,P)+P-VU,+(P-VU,) 

+ —" (Bx P—P xB) 


| l 
~m Nin (E+ : U, x B) —Vq | U 


mNU|5 E+ | Uy xB) V9) 
[dU, dU, 0 
-mN | dt U+ U dt | lay 4 (40) 


where (A)” denotes the transpose of the matrix A. We 
now sum Eqs. (44)—(46) over all the species of the 
system and, making use of Eqs. (19)—(22), (27), and 
(28), we obtain 


00 


52 +V-(@U)=0, (47) 
1 
9 =—V-P+eE+ J xB—oVg, (48) 
and 
oP +V-(Q+UP)+P-VU + (P- VU)! 
1 


+ (Bx P—P x B)—(J—e U) (E+ -U xB) 


1 y) 
(E+ UxB)J—eU)=(.P) 49) 
é coll 
where in writing Eq. (49) we have made use of the fact 
that 
>deNU=)>eN(U—U,)=J—eU,, (50) 
and we have dropped the subscript on U, in Eqs. (47)—(49) 
as it is no longer necessary. If we take the trace of Eq.(49) 
and make use of the relation (31), we obtain 


(570) +-V-(q- >,0U) +P:VU 
J eU).(E+ 1 UxB)—0, (1) 


where the right hand side vanishes since collisions con- 
serve the internal energy of the system. 

We recognize Eqs. (47) and (48) as the well-known 
hydrodynamical equations of continuity and motion. 
Eqs. (47)—(49) are rather useless in their present form; 
they do not constitute a closed set of equations as they 
contain more unknowns than the number of equations. 
One could obtain an equation governing the heat flow 
tensor Q by taking the fourth moment of the Boltzmann 
equation; however, this equation would involve the next 
higher order moment of the distribution function and 
the system of equations is not closed again. 
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THE CASE OF STRONG COLLISIONS 


We now consider the case where collisions dominate. 
Then, to a high degree of approximation, the particle 
distribution functions tend to be locally Maxwellian and 
we may set 


f(x, v, 4) = fo (x, v, 4) [1 -+ ® (x, v, 4], (52) 
where 
m j3 2 
fo (Xs Vy 4) N (x, /) 220 (x /)| 


feo fy WV DEL 53 

exp ys) (x) ff’ (53) 
and @ is considered small compared to unity. The density 
VV, temperature 0, and mass velocity U, are assumed 
to be slowly-varying functions of space and time such 
that their spatial variations over one mean free path, 
and temporal variations over one mean free time, are 
very small. This limiting situation forms the basis of the 
Chapman-Enskog theory of gases and leads to a trun- 
cation of the hierarchy of moment equations. 

The rigorous Chapman-Enskog theory for a plasma 
composed of two or more species, interacting with 
electromagnetic and gravitational fields, is quite com- 
plex. However, the smallness of the electron to ion mass 
ratio and the inability of a plasma to support large 
electric tields—referred to as the condition of quasi- 
neutrality —results in a substantial simplification of the 
Chapman-Enskog expansion technique. 

The condition of quasi-neutrality implies that the 
ratio of the net charge density to the charge density of 
each of the species is very small compared to unity. We 
shall systematically neglect both these ratios compared 
with unity in our discussion of the transport co-effi- 
cients. 

‘Consider for simplicity a plasma composed of .N elec- 
trons and an equal number of protons of masses 7 and M 
and charges —e and +e respectively. Following the 
Chapman-Enskog procedure we can show that, within 
the approximation outlined earlier, the small quantities ® 
are determined by the integro-differential equations (cf. 
Eqs. (3), (16), (52), (53)): 


Mn? 5 M 1. 
| 20 5) w-Vin 04 | @VY—3 w?1):VU, 
5 \E : U, ; B| + Vin (ve) Ww) ft 





(54) 


= (JXB)-wfj + 47 Bewx (VP) f,—VoP 
and 
(5 3) vm +g ww —twer)-00, 
' l6 \E + : U, x B) + Vin (N 0)|-wh 7, 

~ Bewx (V.P)f5-Vy-T>. (55) 


The currents in velocity space I’, obtained by an appro- 
priate expansion of the Fokker-Planck Eq. (16), are 
given by 


8 


2zeélnA | 


_— gi — 
Me \\42wS5 (w)Ss (wr EB. 


V.®* (w) —V.@ | 


r 
‘ 


A 
: \ a3 wif. (w/o (w’ 


u 


w7l Ww w’ . r 
w’ VD (w )} (56) 
and 
. 2ePian iit... . gl 
| OM paws, (wf, (w) oS Be 
V,® (w)—V,,® (w) 
; 2 I , , 
fa (w) Sa Ve (Ww) (57) 
where 
: e (4aNe\! 2 pom 
g-w—w,A=36(-6 | (98) 


Equations (54)—(57) yield, correct to first order in 
@ , the following closed system of moment equations: 


0 
= V-(0 U,) = 9, (59) 
ra) os 
e\,, U,-V)U, r ves oVgq 
as < B—eE V-(P— pb, (60) 
and 
ra) 5 2 
(5+ Uy: Vip- 3PV-U, -3(P pl):VU, 
2 z 1 
3¥-4+30—eU,): (E+ ~ U, x B), (61) 


where p = (NN N_)O = 2 N@ is the isotropic part 
of the stress tensor. In the foregoing equations the 
quantities on the right hand side are determined in terms 
of ® by the equations: 


J > \ Bwwf,D-, (62) 


P—pIl=— Ym J Bwwwl, PD (63) 


and 


b i @:. (64) 


» m\ EA. 
q - \ d?w wn/ 
— . 

The most comprehensive calculations of transport co- 
efficients for a fully ionized gas are due to Marshall using 
effectively the model outlined here. We shall simply 
state his results. Introducing the auxiliary electric field 

1 1 


D=E+—U,xB+5y7,VP, (65) 


and n, the unit vector along the magnetic field, we obtain 


jJ—eU, 
+z nn-VO+x,n 


o,nn-D+o,n x (Dx n)+o,n xD 


(VO) xn] +%,nx VO, (66) 


where the various co-efficients are given by 


( ’ 
O71 OB 305 | 
0, = a0, (w,” tT? + 1.80) . (67) 
0; - 40), T (w, tT? + 4.38) | 
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%, = 1.55% | 
Hs 1.55 ax, (w,? tT? — 0.97) } , (68) 
Xs —4.30 ax), T 
and 
Ne T Net 
% m ’ “0 MW 
( 

; (6 )) 

a w,' t4 4+- 6.28 w,? t? + 0.93 


Here m, denotes the electron cyclotron frequency and 7, 
the effective collision time, is given by 
3 1 ) F (2 (70) 
Fr 10 
(22): \o 2Né@InA 
where, it may be recalled, that 2.N is the total number 
density. 
The heat flow vector is given by 


q 4nn-VO—A,n x [(VO) x n]—A,nx VO 


y,nn-D—yr,nx (D x n)—»,n x D, (71) 
where 
A, = 3.59 dy | 
Ay = dg la (.458 w,? tT? + 3.01) + 4.46 5} 


he 1.25 A, |[bw, t aM, T (w,? tT? + 6.2)] 


¥, = 3.19 m% | 

Yo = avy (0.5 w,? T? + 2.98) (73) 
Ve a Vy (1.25 w,” t? + 7.63) | 

and 

, 2NQt 2NeOr 
4. ’ v 

e nl ° nl | (74 

; 7 
;, 1 | ) 
b —— = 
@,? tT? +- 12.72 


It is clear from Eqs. (66) and (71) that, in the presence 
of a temperature gradient, not only will there appear a 
flow of heat q, but also, a flow of electric current. Simi- 
larly, an electric field produces a flow of heat. According 
to the thermodynamics of irreversible processes, these 
four co-efficients are not independent as three of them 
are linked by the Onsager relations 


0 
*), (75) 


é 


vy, = O\x, 


hm UI 


' 
Ww 


where / takes the value 1, 2, or 


In a representation in which the magnetic field is 
along the z-axis, Eqs. (66) and (71) take the form 


J—«U,=6«-D+x-VO (76) 
and 
q a-VO— vy-D, (77) 
where 
O» o, 0 
o | Os O» o f; (78) 
0 0 oO 


and the dyadics x, A, and y also have similar structure. 
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\ particularly interesting situation results when, in 
the equilibrium state, ¢ 0, and 


j,=0, Vp=0, VO=0, and Ds =0. (79) 


Eq. (78) then yields 


dD, D, (80) 
ee  * 
and 
/; o, Dg, (81) 
where 
o*, a, 
. (82) 
O> 
~. By when @, T So (83 
~6,/1.93 when w,T # (84) 


In presenting the expressions for the components of 
the stress tensor, it is convenient to introduce the co- 
efficient of viscosity, in the absence of a n.agnetic field, 
5 M \"): 


| it, @& 


where T. - 
ere 2 3 (2)' i WI 


1 
Mu 3 NOT, 


and the symmetric traceless tensor 


l 


A : |V Up L(V U,)"| ,V-U,I. (86) 


Then, in a co-ordinate system in which B is along the 
z-axis, we have 
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l zy ij I yr 
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1 oO? T,* / i 2 ( x vy) 
where w, denotes the ion-cyclotron frequency. 
2.1 Weak Magnetic Fields 
In the limiting case of a weak magnetic field wm, ; 7<<1, 


and it is to be observed from the foregoing relations 
that in this limit, to the lowest significant order, 


0, ~ 92, 0,~ 0 % ~ Hs 


A, ~he; 


Expressions (76), (77), and (87) now take the simple 
form 


J—<e«U,=0,D 
q 4,VO—y»,D, (90) 


x,V0, (89) 
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and 


? 
P pl+ zuV-U I -")V Up + (V Up)”| - (91) 


On substituting Eqs. (89) 
state (61), we obtain 


(91) into the equation of 


d 2 . 
(Pe) =3V-4,VO+D) 
s 2 2 a 
Fog AAT 30,78 (92) 
- 
(51),¥9 ' 3V0)-(E re U, > B), 
where y = 5/3. From the Maxwell equationcV = E 


—dB/d¢ we obtain 


oB 


37 = VX (Uy x B)— 5 VxJ—eU,). (93) 


Thus our complete set of hydromagnetic equations 
are: 


e 7+ V. (0 U,) QO, (94) 
| . | | I VY 
o(5,+ U,: ae ‘Vp +eE+—JxB- -0Vq 
-}- uv? Uy {- - bt VV ° Up, (95) 
oB e 
57 ~ ¥ X (Uy x B)—, Vx J—eUy), (96) 


1 


0 2 
@” (5, + WV} (pe ”) 3 VA VO + % D) 


2 


4 2 
FRA 36 J—e«U, (97) 


1 
: Vp+2 x ve). E + —-U, x B), 
K Wi 3 c 

where the charge and current densities are determined 
by the Maxwell equations. For an ideal, perfectly con- 
ducting fluid, 0, — 00, 

1 

E+ U, x B=0,u=0,q = 0, 

and if the characteristic frequencies of interest are much 
less than the ion cyclotron frequency and the foregoing 
relations reduce to the well-known hydromagnetic 
equations: 


se +V-(@U)=0, (98) 


/0 
e(3,+U-V)U 


1 
{ - (J—«U) x B—o Vo, (99) 
V x (U x B), (100) 
and 
d | 
“(Pe | 0, (101) 


where in writing these equations we have dropped the 
subscript on U,, as it is no longer necessary. 
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2.2 The A. C. Conductivity 


There is one further important transport co-efficient 
which has not yet been worked out for a plasma. Since 
one is often interested in frequencies of the electro- 
magnetic fields comparable to the effective collision 
frequency, the usual Chapman-Enskog theory needs 
some modification. This can be effected readily by intro- 
ducing terms in 0@ * /0¢ on the left hand side of Eqs. (54) 
and (55). Now it is the electrons which primarily de- 
termine the conductivity. For the purpose of an order- 
of-magnitude estimate, we can neglect electron-electron 
collisions. Then, in the absence of temperature and 
density gradients, magnetic fields, or mass flows, Eq. (55) 
reduces to 


oD é 
fo) (97 + @E-¥] 
O\':1 3 [#I—ww = 
‘ 7 
(=)"= TOW | ny Jo 9 (wv) ’ (102) 
where 
mW 2 
"a 6 2aNeélnA (103) 


If E varies in time like exp (— / w/), then so must D, 
and the solution of Eq. (102) is 


~ [(m w)?/20]"/ 
q , 1 104 
‘2074 + Z 1 MT, (mn?/O)': (20%) 
The associated current is then given by 
J=—c\@wwfo(w)D, (105) 


which, on using (O/m)": as the unit of velocity, can be 
written in the form 


J=cE, (106) 
where 
+9 eo dwwiewh 
e : N é Ty 1 1 1 ; (107) 
m (22)'l23 1+ 5 OTUW 
If wt, <1, Eq. (107) yields to the lowest significant 
order 
Net, 1 1 945 . 
_ yt): 0 
o m (nyh3 (32 4 fO%); (108) 
whereas if wt, > 1, we obtain 
iNet, 1 3a 1 
_— | . Q 
ss mo wT, 3 a es (105) 


3. THE CAsE oF WEAK COLLISIONS 


If the plasma is so dilute, or the temperature so high, 
that the direct interparticle collisions of the first and 
second class described earlier are negligible, then the 
particles interact mainly via the macroscopic electro- 
magnetic fields. We may thus drop altogether the col- 
lision term in the Boltzmann equation (3) in treating the 
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dynamics of the plasma. Firstly, we wish to emphasize 
that this description of a plasma via the collisionless 
Boltzmann equation is entirely equivalent to the solution 
of the equations of motion for the particle orbits. This 
can be shown as follows. 


3.1 Equivalence of Collisionless Bolzmann 
Equation and Particle Orbit Theory 
Consider the collisionless Boltzmann equation or 


what is referred to as the kinetic equation in the Russian 
literature: 


os t(v-V)f/+F-V,f=0. (110) 
The Lagrangian subsidiary equations are: 
dv 
F qt’ (111) 
and 
dx 
? 
dt- (112) 
Let the solutions to these equations be 
V = V (41, %q°-- XG, ¢) (113) 
and 
X = XK (0 4 gy +> Xe, f), (114) 


where the «’s are the integration constants. We further 
suppose that Eqs. (113) and (114) can be solved for the 
x’s. Now the general solution of Eq. (110) is an arbitrary 
function of the constants of motion, i.e., 


F(& v4) = £ (0, 4° +g); (115) 
for substituting Eq. (115) into (110) we obtain 
6 
of da; , 
a," a7 be v-V x; + F-V,0;| 
i=1 
6 
Ot dx, i 
da, qt 0 (116) 
i=1 


since the a’s are constants along a particle trajectory. 
Thus we have shown that the most general solution of 
Eg. (110) is an arbitrary function of the integrals of the 
Lagrangian subsidiary equations (111) and (112). But 
these are just the equations of motion for the particle 
orbits. The equivalence of the particle orbit theory and 
the collisionless Boltzmann equation was first demon- 
strated by Jeans, and in the astronomical literature is 
referred to as Jeans’ theorem. When the «’s are not true 
constants but adiabatic invariants, the resulting / satis- 
fies the Boltzmann equation to within the accuracy of 
these invariants. The considerations of the orbit theory 
provide a better insight into the physical factors involved 
in a particular problem. However, the Boltzmann 
equation provides a rigorous treatment of complex pro- 
blems which may not be easily accessible to the orbit 
theory. 
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3.2 The Low Temperature Approximation 


Neglecting, then, the right hand side of Eq. (49), the 
equation for the stress tensor now becomes 


o.P4V-(Q+UP)+P-VU+(P-V Uy 


+ “(Bx P P x B)—JJ—eU) (E+ Ux B| 


MC 


] 
-|E+—U x B|(J—eU)=0. (117) 
In order of magnitude |0P/d¢| ~ w»,,3 and |\V-Q 
~ kv,,* where v,, = (O/m)!'2, and w and & denote the 
characteristic frequency and reciprocal length respec- 
tively. Therefore, it is clear that if 


@ 


hm "tn? (118) 


the heat flow tensor can be neglected in Eq. (117). For 
the linearized equations of motion w/& is the phase 
velocity of the wave, while for a non-linear system it 
might be termed the effective phase velocity. Thus if the 
phase velocity of the wave is much larger than the 
thermal speed of the ions or electrons, Eq. (117) re- 
duces to 


oP +V-(UP)+P-VU+(P-VU)' 


é 


4. (B 


MC 


P P x B)—(J—eU) (E+. U xB 


—|(E+ "Ux B) (J éU)=0. (119) 

Thus, by virtue of the approximation (118) we are 
able to obtain a closed set of moment equations. This 
scheme will be referred to as the Low Temperature Ap- 
proximation and will yield results correct to first signi- 
ficant order in the parameter of smallness 4v,,/@. 

Clearly, one can decouple the sth moment equation 
from its dependence on the (” + 1)th moment in an 
analogous fashion, and obtain the corresponding higher 
order thermal corrections. In all such calculations care 
must be taken to discard any spurious roots which may 
occur, that is, roots which are incompatible with the 
basic assumption (118). 


3.3 Chew, Goldberger and Low Approximation 


If we now assume that the magnetic field is so strong 
that the electron or ion cyclotron frequency is very large 
compared to the characteristic frequencies of interest, 
we can make a systematic expansion of the Boltzmann 
equation (110) in inverse powers of B. This has been 
carried out by Chew, Goldberger, and Low. An equi- 
valent procedure is to carry out this expansion in the 
equation of motion (119) for the stress tensor. If we 
assume that the Joule dissipation is negligible, i.e., 


(E+ =U xB)J—eU) | - BP), (120) 


then the leading term in Eq. (119) gives 
BxP=PxB. (121) 
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It immediatel» follows from this equation that 


P—P.1+(P,—P)nn, 


(122) 
where n denotes a unit vector along B. Eg. (122) im- 
plies that the pressure tensor is diagonal in a local rec- 
tangular co-ordinate system, one of whose axes points 
along n. In the plane perpendicular to n the pressure 
is a scalar of magnitude P,, and the pressure along B is 
P. which, in general, is not equal to P,. Equation (119) 
now reduces to 


d 


77 (PLL+ (P,P) nn] + [PI 


(P,—P,) nn] V-U + [P,1+ (2, —P,) anl- VU 


+ {P,1-+ (P,—P) an]-VU}"—0, (123) 
where we have now set 
E + lu X B—0, (124) 


Eq. (124) implies that the medium has very high elec- 
trical conductivity so that it is unable to support, in its 
rest frame, any significant electric fields. Taking the 
trace of Eq. (123), we obtain 


d 1 \ : 1 . ; 
a7 Pi tz Puls \?, +5 P,|V-U LP,V-U 
| (P,—P,)nn:VU=0. (125) 
The nn component of Eq. (123) is 
/ 
7 P+P,V-U+2P,nn:VU=0. (126) 


From Eqs. (125) and (126) we readily obtain 


4p 42P,V-U 


A -P,nn:VU=0, 
dt , 


(127) 


Using the relation (124) in the Maxwell equation 


cV xE 0B/07, we obtain 
oB 
\ ° 28 
ry, V x (U x B). (128) 
It is convenient to write Eq. (128) as 
dB re) 
a= last U-V)B=(B-V)U—B(V-U). (129) 


From Eq. (129) and the equation of continuity 


do 
-~+toV-U=0, 


Q 
dt (130) 


it can be shown readily that 


od 4 


nn: VU Bdt\o 


(131) 
where B denotes the magnitude of B. On using this 
result in Eqs. (126) and (127), we obtain the two 
adiabatic relations 


d |P,, B 

ai 08 0 (132) 
and 

ad; P,\ 

Zi\oB) ™ (133) 
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Starting from the Boltzmann equation (110), Chew, 
Goldberger, and Low have derived the same adiabatic 
relations on the assumption that the heat transport along 
the lines of force is negligible. 

Equation (133) implies that, if one moves along with 
an element of the fluid, the average transverse kinetic 
energy per particle is proportional to B. This is com- 
patible with the well known constancy of the magnetic 
moment of a charged particle in the adiabatic approxima- 
tion. It can be shown that Eq. (132) implies the con- 
stancy of the longitudinal invariant of Chew, Gold- 
berger, and Low. Thus for a low density plasma in the 
presence of a strong magnetic field, in the framework 
of the low temperature approximation or on the as- 
sumption of negligible heat transport along the lines 
of force and when the conductivity is very high, we 
obtain the following complete system of equations: 


00 eS 
5, + V-@U)=9, (135) 
dU 0 
Oar 2 \ar +U-V)U 


V-P+ : (J—«U) x B—oV9q, (136) 


oB ~ 
rv, V x (U xB), (137) 
P—P,14+(P,—P)nna, (138) 
d (P,, Bt 
dt 0° 0, (139) 
and . 
d{P,\ 
qi\oB) 0, (140) 


where the charge and current densities are determined 
by the Maxwell equations. 


III. Wave Phenomena in Plasmas 


1. GENERAL REMARKS 


We shall now discuss the behaviour of plasmas for 
states initially close to a steady state. This implies the 
existence of a time-independent, self-consistent solution 
of the appropriate equations of motion, jointly with 
the Maxwell equations, subject to proper boundary con- 
ditions. The notion of closeness implies the existence of 
some dimensionless parameter of smallness, in terms of 
which one can expand the time dependent quantities of 
the equations governing the system. In most situations 
the appropriate parameter is the ratio of the macro- 
scopic velocity to the phase velocity. The lowest order 
system of the expanded equations is just that governing 
the steady state. To next order one obtains a set of linear 
equations, whose coefficients do not in general depend 
on time, which describe the small motions about the 
steady state. 

If there exists no small initial perturbation whose 
amplitude grows monotonically or sinusoidally in time, 
then the system is said to be stable; otherwise, unstable. 
In the absence of any dissipative mechanisms the small 
motions of a stable system may tend to persist in- 


36) 
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38) 


39) 


40) 
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-Oon- 
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definitely; whereas in the presence of a dissipative 
mechanism they may tend to die out, the system re- 
turning to the steady state. 


If the system is unstable, after a certain interval of 


time the amplitude of the motion will become so large 
that the linearized equations will cease to be valid. The 
subsequent behaviour of the system will be governed by 
non-linear equations, whose analysis, in general, is dif- 
ficult to carry out. 

There are two plausible situations into which these 
instabilities might possibly carry the system. The first 
of these might be termed a turbulent state, characterized 
by randomly fluctuating electromagnetic fields and 
motions, which could manifest itself by an enhancement 
of diffusion and resistivity over and above that resulting 
from collisions. The second situation is associated with 
instabilities whose exponentiation time is short, can 
affect significantly the heavy particles and may lead to 
coherent macroscopic motion of the plasma, resulting 
in a radically different final state. The former of these 
has been conjectured as a cause of the loss of contain- 
ment in the stellarator type of controlled fusion devices. 
\n example of the latter has been observed in pinch 
effect experiments where the discharge has destroyed 
itself upon achievement of a critical current. 

Since the coefficients of the linearized equations do 
not depend explicitly on time, one can seek a complete 
set of characteristic solutions whose time dependence is 
of the form 


exp (—/w/). (1) 


This leads to a characteristic value problem for @ which, 
in general, may be complex. If there exists an w for 
which 


Ino>0, (2) 


then the system is unstable. If the coefficients also 
happen to be independent of one or more of the space 
variables, one can seek solutions of the form exp (ééx). 
The solution of the characteristic value problem, sub- 
ject to the appropriate boundary conditions, then leads 
to a relation between m and & which may be written 
implicitly in the form 


iin, =o. (3) 


A relation of the form (3) is usually referred to as the 
dispersion relation. The solution of Eq. (3) for @ in terms 
of & would immediately disclose whether the system is 
stable or unstable. 

Sometimes it is convenient to characterize the family 
of steady states by a certain parameter, say #. As one 
changes f continuously, the system may pass from 
stability to instability for some determinate value of /, 
say Pp = f,; 1.e., w» passes from the real axis into the 
complex @ plane (with Imw > 0) at § = p,. The state, 
characterized by 8, which marks the boundary between 
the stable and the unstable regions is referred to as the 
marginal state. 

In the case of static equilibria for non-dissipative 
systems governed by the hydromagnetic equations, it 
has been shown by Bernstein, Frieman, Kruskal, and 
Kulsrud that the linearized equations of motion are self- 
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adjoint; consequently the characteristic values m? are 
real. Therefore, 
o- > 0 (4a) 


corresponds to stability and 
o- <0, (4b) 


to instability. However, for more general types of 
equilibria, the characteristic values «m may be complex. 
In case the imaginary part of @ is positive, the system 
will be unstable; following Eddington, this type of in- 
stability is referred to as overstability. 

In addition to questions of stability which are usually) 
posed for isolated systems, there is the important class 
of driven problems. In these situations a system in a 
steady state is excited by an external agency, perhaps at 
a fixed frequency. When the amplitude of excitation is 
small, the problem can again be treated by the same set 
of linear equations as those employed for the analysis of 
stability, but with appropriate modification of the 
boundary conditions. In the cases where a Fourier 
analysis in space is in order, one is again led to a dis- 
persion relation, but now it is to be viewed as determin- 
ing the wave number & (in general complex) as a func- 
tion of the prescribed frequency wm, assumed real. The 
fact that & can be complex permits waves to grow or die 
out in space, giving rise to the phenomena of amplifi- 
cation and evanescence. Problems of reflection, trans- 
mission, scattering, and amplification of waves by plas- 
mas fall into this category. It must be emphasized that 
when any particular problem is properly posed with due 
regard to boundary conditions, there will be no con- 
fusion about the interpretation of the resultant dispersion 
relation. 


2. UNBOUNDED PLASMAS IN THE ABSENCI 


or EXTERNAL FIELDS 


2.1 Loneitudinal F:lectron Plasma Oscillation 


We first consider the simplest case of electron oscil- 
lations in a uniform plasma, neglecting the thermal 
motion of the electrons. The ions will be assumed to 
be at rest and to form the uniform neutralizing back- 
ground for the electron charge in equilibrium. Let N 
denote the electron density (which is equal to the ion 
density) in equilibrium. In the perturbed state let the 
density be denoted by N + 2, where n is considered 
small so that all terms quadratic or higher in # can be 
neglected. The equations of motion now are (cf. Eqs. 


(II-44) and (II-45)): 


On y . 
NV-v=0, (5) 
o7 
and = 
ay 
nl cE, (6) 
0s 
where —e and m denote the charge and mass of the 


electron respectively. For longitudinal oscillations the 
electric field is determined simply by the Poisson equation 


V-E 47é 


-4 en. (7) 
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Eliminating v and E from Eqs. (5)—(7) we obtain for 
the equation of motion 


O72 4aNen P 3 
= = = D,2n, 
or ml sa liad (8) 
where 
42 Ne*\'/: - 
o 22 X 8980 N ’*sec-!. (9) 
? m 


It follows from Eq. (8) that the density fluctuations 
oscillate sinusoidally with the characteristic frequency 
w,. The frequency mw, is referred to as the plasma 
frequency. Formula (9) was first obtained by Tonks and 
Langmuir in 1929. 

The absence of space co-ordinates in Eq. (8) shows 
that these waves are non-dispersive —their group veloc- 
ity vanishes — and there is no tendency for a wave 
packet of this type to propagate through the plasma. 
As a result we can specify the phases of the electron 
displacement in such a way as to obtain a traveling wave 
—a wave, however, which moves continuously through 
a fixed region without ever progressing beyond. 


2.2 Transverse Electron Plasma Oscillations 


We shall now consider the transverse (or electro- 
magnetic) plasma oscillations, neglecting the thermal 
motions of the electrons. The equations governing the 
problem are the linearized equations of motion (cf. 
Eqs. (II-44) and (II-45)) and the Maxwell equations 


On : 
ry, +NV-v=0, (10) 
ov 
m >, eE, (11) 
4a Ne 10E 
4 ‘ wits ? 
VxB Ca ak TZ (12) 
and 
10B 
VxE wt Th (13) 


From Eggs. (10)—(13) we readily obtain 
er 


dz —w,*E —AVxVxE. 


(14) 


Taking E to be of the form exp/(—m/+k-r), 


Eq. (14) leads to 
(w? — w,? — c? k*) kxt&E 0, (15) 


(16) 


For transverse oscillations k « E does not vanish and 
Eq. (15) gives the dispersion relation 


and 


(w? — o,2) (k-E) = 0. 


w? = w,? + ck; (17) 
while for longitudinal oscillations k - E does not vanish 
and Eq. (16) gives the dispersion relation already dis- 
cussed. In contrast to the longitudinal oscillations, the 
transverse oscillations propagate in the plasma with the 
group velocity. 


Ow ck 


“9 dk &[1+(,2/ 2)" 
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2.3 Effect of Thermal Motions on Electron 
Plasma Oscillations 


We shall now treat electron plasma oscillations, taking 
into account the thermal motions of the electrons using 
the low temperature approximation developed in 
Section II. In the equilibrium state, we shall assume the 
plasma pressure to be a scalar; however, we do allow 
for anisotropy of the material stress tensor during the 
perturbations. The perturbed state is characterized by 


pPI+ Pr (19) 
where » = N@ denotes the electron pressure in equili- 
brium, and p,, the perturbation in the stress tensor. The 
linearized equations of motion now are (cf. Eqs. (I-44) 


to (II-)46): 


v, N+4a4, 


Ow. ., 
+ NV-¢ 0, (20) 
OV ; 

m IN =< -V-p,—NeE, (21) 


and 
0 
oP =—P@-vI—plVv+ (Vv). (22) 


Equations (20)—(22) must be solved in conjunction 
with the linearized Maxwell equations 


42aNe 10E 
VxB=— - v+ . e (23) 
and 
10B 
VxB=——7Z- (24) 


Introducing &, the Lagrangian displacement of an ele- 
ment of the electron fluid from its equilibrium position, 
to a first approximation, we can write 


(0 
“Ot 


Vv 


(25) 
Eqs. (20) and (22) can then be readily integrated to give 


n=—NV-E (26) 
Pi = —p(V-&)I—p[VE+(V)’], (27) 


where the constants of integration have, by definition, 
been set equal to zero. Assuming the space time behaviour 
of all quantities to be of the form exp / (—w/ + k-r), 
the foregoing equations yield 


0 6. e 
2,,kk-6+7@5—_ VE, (28) 


and 


we 


and 


@ekxkxE= oe 42a New? E —w'E > (29) 


From Eqs. (28) and (29) we obtain for longitudinal 
oscillations the dispersion relation 


(30) 


while for transverse oscillations we obtain 


. & 
(wo? — c2 &?) |? = is) + ww,? = 0. 


Mm 


(31) 
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Eq. (31) leads to a quadratic in mw? whose roots are 
given by 


2__! ar rn. : we. ore, RO 
oe 2 \o,? | ¢2 £2 4 - + , lo, c2 R2 — 
42° RE O|m \ 
fot + 2+ POlmpPs (32) 


In our approximation m? >> £2 O/m. Retaining terms to 
the lowest order in 4? O/m, we obtain for the two trans- 
verse modes of oscillation (distinguished by the sub- 
scripts 1 and 2): 


P . on , oe o,” 
C4? = Wy! + CRE + m o,2+ ¢c2k (33) 
and 
P RO 1 34 
2 m 1+ (w,?/c? k*) (34) 


The root w,” corresponds to the transverse oscillation; 
however, the root @,” is incompatible with our basic 
approximation that w?> &? O/m and therefore must be 
discarded. 

It may be remarked here that if, instead of using the 
equation for the stress tensor, we had used the classical 
adiabatic relation 


d ° 
dt (po-*) =0, (35) 


we would have obtained for longitudinal oscillations the 
dispersion relation 
w=o,?+ y — ke (36) 

instead of the relation (30). However, in the absence of 
collisions, it is necessary to allow for anisotropy in the 
pressure during the oscillations. 

The effect of thermal motions on electron plasma 
oscillations was first considered by Thomson and later, 
by Bailey. They obtained the dispersion formula 


wo? = w,* + — & (37) 


for longitudinal oscillations. It can be seen readily that 
the appearance of the factor O4?/m instead of 3 O4?/m 
or y OR? /m results from neglecting the change in tem- 
perature of the gas resulting from the density fluctua- 
tions. This amounts to considering the oscillations to 
be isothermal and not adiabatic. This is incorrect for the 
high frequency electron plasma oscillations; for a con- 
sistent approximation it is necessary to use the complete 
set of equations (20)—(22) to treat the problem. 

In part II of this survey this problem will be discussed 
from the standpoint of kinetic theory, and it will be 
seen that expressions (30) and (33) determine correctly 
to first order in O the real part of w. The higher order 
thermal corrections can be recovered by going to higher 
order moment equations. However, an essential con- 
sequence of the treatment from the kinetic equation in 
the absence of collisions, which is not recovered on 
using the truncated set of moment equations devoid of 
collision terms, is the phenomena of Landau damping. 
That is, the kinetic equation yields for isotropic equilibria 
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a negative imaginary part effectively proportional to the 
number of particles moving with the phase velocity of 
the wave. This number, of course, tends to be ex- 
ponentially small since, by assumption, the thermal 
spread is much less than the phase velocity in question. 
It must be emphasized here that the longitudinal and 
transverse oscillations are strictly uncoupled only in the 
case of a non-relativistic plasma and in the absence of 
any external magnetic fields, temperature or density 
gradients. The presence of an external magnetic field or 
inhomogeneities in plasma density and/or temperature 
result in the coupling of longitudinal and transverse 
modes, and the behaviour of the plasma, in general, is 
quite complex. This will be discussed in detail later. 


2.4 Ion Oscillations 


So far we have considered only electron plasma oscil- 
lations on the assumption that these oscillations are too 
rapid for the heavy ions to follow, which implies that the 
ions, therefore, may be considered at rest. Another class 
of oscillations, which is possible in a collision-free 
plasma, is the so-called ion oscillations. These are so 
slow that electrons see them as quasi static and con- 
sequently are distributed according to the Boltzmann 
distribution. During the oscillations it is appropriate to 
treat the ions and electrons as having different tempera- 
tures. However, we shall restrict ourselves to the case 
of longitudinal oscillations only. 

The linearized equations of motion now are (cf. 


Eqs. (I1-44)—(II-46)): 
On; 
ap th V-v=0, (38) 
OV — 
MNy,, V-p,+NeE, (39) 
") 
P= —AV-vi—pVv+(Vv)], 40) 
V-E=42e(n,;—n,), (41) 
and 
n, = N\exp A -1) (42) 


where #,; and ”, denote the perturbation in the ion and 
electron equilibrium densities (assumed equal) respec- 
tively, O, is the electron temperature; p, = NO, where 
Q, is the ion temperature, p, is the perturbation in the 
ion material stress tensor, and M is the ion mass. g 
denotes the electrostatic potential defined as E V¢. 
Following the normal mode procedure outlined in § 2.3 
of this section, we obtain for the dispersion relation 
— 


= TH (kAy) 2008 +3 oF 


where w,,,, the ion plasma frequency, is given by 
4a Ne*\": 
Oy; = |p 


\ 


(44) 


and we recall that 4,, the Debye shielding distance, is 
given by 

—0, "s 45 
Ao (3 > nal (45) 
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Physically the reciprocal of /, acts like a screening radius 
for electrical forces in ionized media. Thus the potential 
in the neighborhood of a sphere charged with ¢ units of 
electricity is effectively given by 
6 A 
p= Lexp(—Ayr) (46) 
in an ionized medium near thermal equilibrium, com- 
pared to q/r in a dielectric. 


If kA, 1, ie., the wavelength of the disturbance 

is large compared to the Debye length, formula (43) 
reduces to 

0 


wo —w,2k4,2 +384 


I . 
2 = O-+ 30) &. 


(47) 
It is clear that if the ion temperature is very small 
compared to the electron temperature, the effect of the 
ion thermal motions can be neglected. We shall show 
in part II, using the kinetic description, that if O; is not 
small compared with 9,, the ion oscillations are strongly 
damped. It may also be observed from Eq. (44) that the 
ion plasma frequency is much smaller than the electron 
plasma frequency. 
The low frequency ion oscillations differ from the sound 
waves in the mechanism responsible for the organized 
motion. lon oscillations are produced by long range 
Coulomb forces whereas sound waves are produced by 
short range collisions. 


3. UNBOUNDED CoLpD PLASMAS IN THE PRESENCE OF AN 


EXTERNAL MAGNETIC FIELD 


We shall now consider oscillations in a uniform 
plasma in the presence of a constant external magnetic 
field B,. We shall assume the plasma to be at zero tem- 
perature; this amounts to neglecting the thermal motions 
of the ions and electrons. The effect of thermal motions 
will be taken into account quite rigorously via the 
Boltzmann equation in part II]. We shall further assume 
that each constituent of the plasma has in the equili- 
brium state a constant velocity u. The motion of the 
constituents of the plasma will, in general, give rise to 
a current distribution J. We shall, however, assume that 
the magnetic field which results from J is negligibly 
small compared with the external magnetic field. 

The linearized equations of motion for each of the 
species of the plasma are*: 
oY 4 (u-V)v ‘(e+tuxB+ iv 


Wi c c 


B, |, (48) 


On 


Th V-(Nv- 


eu) = 0, (49) 

where the perturbations in the electric (E) and magnetic 

(B) fields are given by the Maxwell equations 
4% 1O0E 


VxB=—J+ 73>) (50) 


* The notation used here (u for the velocity in the steady state 


and v for the perturbation in the velocity) is different from that 
used in section II. This should, however, cause no confusion. 
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10B 


VxE > >?’ 


(51) 
and 
J=D>e(Nvinn). (52) 
The summation Eq. (52) is to be carried over both the 
species of the plasma. The cyclotron frequency is defined 
to be 
eB, 53 
Oe . J: 
’ MC (93) 
We shall now assume that all quantities have space 
and time dependence of the form exp/ (—w/ + k-r). 
Eqs. (48) and (49) then yield 
— ) 
‘Vv 
(o k- u) 
and 
] 


. ( A i c 


k-u)v u > B) w,Xv- (55) 
The appearance of the factor m —k- u instead of w is 
typical of a current-carrying plasma; the term k - u ac- 
counts just for the Doppler shift in the frequency m due 
to the motion of the electrons or the ions. From Eqs. (51) 
and (55) we obtain 


€ 


i (w—k-u) (w,- v) = (we: E) + 
é ~ 
5, we [a x (k x E)] (56) 
and 
-i (ow —k-u) (w, X Vv) 5) (We < E) 4 : 
Wyte X [w X (k x E)] —w, (w.-¥) +2. (57) 


Using the results (56) and (57) into Eq. (55), we obtain 


e ] 
3 MW,” —(w—k -u)* 
ln! . RE) 
i(o—k-u) e rou x (k x E); 
1 
w, X E- “we * {e x (k > E)| 


af = -u) Kw, -E)+ * we [u x (k x E)|}| (58) 


From Eqs. (52) and (54) we obtain for the current density 
due to oscillations 
k-v 


wo—k-u 


J=SeN(v+u (59) 


Equations (50), (51), and (59) lead to 


(w* — c*? k*) E + c?k (k- E) 
. k-v ( 
w—k-u . 


L tzio DeN(v+u (60) 


where v is given by Eq. (58). On substituting for v 
accordance with Eq. (58) into Eq. (60) we obtain, 
general, an equation of the form 


R-E-—0, 


n 


_- 


n 


(61) 


is 


is 


:) 


where R is a 3 x 3 Hermitian matrix. In order that 
Eq. (61) have a non-trivial solution, we must demand 
that the secular equation 


R 0, (62) 


be satisfied. This results in the required dispersion re- 
lation. We shall now discuss the implications of Eq. (62) 
in several special cases. 


3.1 Oscillations in a Static Plasma 


We shall first discuss the case of oscillations in a static 
plasma, i.e., where the mean fluid velocity vanishes 
everywhere. Setting u = 0 in Eq. (58) and substituting 
the resulting expression into Eq. (60), we obtain an 
equation of the form (61). Taking k to be along the 
z-axis of a Cartesian system of co-ordinates and assuming 
that the magnetic field lies in the xz plane making an 
angle @ with the z-axis, the elements of the matrix R are 
given by 





. et os Bn Conall 
R 2___ 242 2 S Mp: to¢" COS" O 
ua o ~s T / w,2 2 
Laut wo 
- > @o,cos 6 
Ry y= Yias—e 
= - , oO," — w* 
— 9 
« oO” 
Ris + Rs) Bs, 0," —z 5 sinO cos 
‘ , hand oO." ow ( 3) 
6 
ie — = 
| ow -¢2 & 3 >” @ s 7 S @ = 9 9 
= ° e fe ah ow 
XX , wo,sin$ 
Rog R39 i> Oy" 2. 2 
is , wo 
- “ 2 cin? 
R a a a, es... ] 
33 @ - Jf DB T / Dy 2 2 
os kommt , 7) 
where w,,” = 4 a Ne?/m 


3.11 No Magnetic Field Case 

As the simplest illustrative example of Eq. (62), con- 
sider the case of electron oscillations. Setting w, = 0 in 

Eq. (63), we obtain Eq. (64), see below on this page, 
where w,, now denotes the electron plasma frequency. 
We thus obtain for longitudinal oscillations (k along E) 
r=, (65) 
and for transverse oscillations (k perpendicular to E) 


wo o 3 |. ¢2 £2. (66) 
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3.12 Propagation along the Magnetic Field 


For propagation along the magnetic field, 6 — 0, 
and Eq. (61) takes the particularly simple form (cf. Eqs. 
(61) and (63)) of Eq. (67), see bottom of this page. 


For longitudinal oscillations we thus obtain the relation 


" 1 
wo oo oo,8 == 9,511 +4 7 |= Oh." s (68) 
y 


where w,, denotes the electron plasma frequency. The 
transverse oscillations are circularly polarized and are 
described by the equation 


(69) 


M 


where 


ee oe (70) 


the plus or minus sign corresponding to the right or left 
handed circularly polarized waves respectively. The 
oscillations are, therefore, governed by the dispersion 
relation 


uo c2 k2@ — ow, 2 7) 
Oni” (0). (71) 


for the right handed circularly polarized waves. Here 
the electron and the ion plasma frequencies are denoted 
by w,, and w,,, respectively, while the electron and ion 
cyclotron frequencies are denoted by w, and w, respec- 
tively. A relation similar to (71) applies to the left 
handed circularly polarized waves. 


a. High Frequency Oscillations 
If w > w,, the dispersion relation (71) reduces to 


a? — c2 k2 —w,? o,°=0. (72) 


In the limit when wm w,, Eq. (72) yields, correct to 


the lowest significant order in w,/a, 


ao ce T Moe I (wo, + c® &)! 2 : (73) 


It is clear that a necessary condition for the validity of 


which are in accordance with Eqs. (8) and (17). this relation is that w,? < w,,? + ch?. 
w? — w,? — ck 0 0 E, 
2 2__ 248 
0 w*? — ow,” k 0 5. 0 (64) 
0 0 w*—wo,* || I 
-_ 9 9 —_ 
oO," wr . , Ow, 
w? — ¢2 £2 + S = - i> w,” ——~» 0 E 
hag,” —~- OF doug @,* (* : 
he ’ OO, a er % : wo ~ 
7 o 2 ~ c " ot— ck ™~ o,2 - = 0 E. 0. (67) 
>? o,2 —a?* La *® ww? —o? ‘ 
0 0 ow > o,2 E. 
id 7 
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b. Oscillations near the Electron Cyclotron Frequency 


Another case of interest is the one when the frequency 
of oscillations is close to the electron cyclotron frequency. 


In this case we find that Eq. (72) leads to 


2 
1 ; D pe 


w,? — (c2 k8 + w,,?) 


7) -W, 


(74) 


It is clear that the necessary condition for this approxi- 
mation to be valid is that 
Wy << | w,? 


—c2 k2 |, 


We can write Eq. (72) in the form (for w > w,) 
1 + (@,/@) 


V\2 ; 
| c | 1+ (w,/@) pene (@y¢2/@?)’ (76) 


«w/k& is the phase velocity of the wave. The 
behaviour of the phase velocity as a function of the fre- 
quency of oscillation is illustrated in Figs. 1—3 for 
some values of the field strength. 


(75) 


where |” 





45r- 


-_—--- 








e7 
| 








eooe” 
Lee 








Fig.1 ‘The behavior of Z = (V/c)*, where V is the phase velocity 
of the wave, as a function of x = w/w, for electron oscillations for 
propagation along the magnetic ficld for the case when w, = 0.15. 
The negative values of Z correspond to imaginary values of &, 


which correspond to absorption of waves in a bounded plasma. 
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Fig. 2 The behavior of Z = (V/c)? as a function of x = w/w, 
for electron oscillations for propagation along the magnetic field 
for the case when w, = w,. 
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Fig. 3 The behavior of Z = (V/c)? as a function of x = w/a, for 


electron oscillations for propagation along the magnetic field for 
the case when w = 10 @,. 


c. Oscillations much below the Ion 


Cyclotron Frequency 


Now consider the low frequency oscillations such that 
the condition m <w, is satisfied. Then to the lowest 
significant order we can write Eq. (71) as 


@ 


“ 2 £2 > O 7) . 7 
w* — c* R* + Wo; | eo Woe %, 0. (77) 


Using the fact that w,,?/a, 
c2 4a N M/B,?, Eq. (77) leads to 


2) ; 2 
Wve /@, and Opi WO; 


P c2 B2 " 
w i + (@/A%’ (78) 
where 
A By 


(42 M)'h (79) 
denotes the Alfvén speed. Eq. (78) is the dispersion 
relation of Astrém for the extra-ordinary hydromagnetic 
wave. It must be emphasized here that the hydro- 
magnetic wave, in principle, is just a special case, in the 
appropriate frequency region, of the well-known trans- 
verse electromagnetic waves. 

In the limit when ¢/A > 1, Eq. (78) reduces to 


“i 

la} ° 
that is, the phase velocity of the wave is just the Alfvén 
speed. In the other limit when ¢/A < 1, we get 

= ¢2 £2 , (80a) 


which corresponds to the usual electromagnetic modes. 


fB; 


(80) 


wo 


d. Oscillations near the Ion Cyclotron Frequency 


Let us now consider the case when wm ~@,. This 
implies that w <w,, and the dispersion relation (71) 
reduces to 


s OO; 


—_— 
pe Wi 


w? — c2k2 —@w 0. 


(81) 
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This leads to the frequency of oscillation The behaviour of the phase velocity as a function of the 
" frequency corresponding to the roots of this equation 
Ws” . . . si 2 = 
wo~o,|)1+-s (22 | (81a) are illustrated in Figs. 4—6 for some values of the 
w2 —(c2 #2 +- aw 
i pi 


field strength. 


The condition for the validity of this approximation is 








+40 
that 
w2>w,2 + 2k. 30 b a 
‘sem phe — +20 | 7 
3.13 Propagation Transverse to the Magnetic Field 
For propagation perpendicular to the magnetic field, *° [ 7 
7] 2/2, and we obtain the secular Eq. (82) see below. 
Z-0 2 a - -- 





One of the modes of oscillation corresponds to the usual 











transverse oscillation with electric field along B, and -~ | 4 
satisfies the dispersion relation 
F -20 “ 
aw? — c2k2 + 3 w,? ~ 7k + w,,?. (83) 
r . : 5. 7 -30 F - 
The other two modes of oscillation are determined as 
roots of the equation (84), see bottom of this page. -40 





Consider the case when w < chk. Then to the lowest Fig.4 The behavior of Z (V/c)*?, where V’ denotes the phase 




















significant order, the 22 element of this determinant velocity of the wave, as a function of x = w/w, for electron 
must vanish. We thus obtain oscillations for propagation transverse to the magnetic field for the 
case when o, 0.1 w,. 
wo, ,," . 
" pi 1 pe rc 
. se eg 5=0. 85 
 o2— oo * w2— oe (85) 
12 
In the limit when w,; < w < w, < w,, Eq. (85) reduces to 
9 
9 9 
WO, 7 | Wu? 
pi pe 
me 3 0. 86 
w ' wo? (86) . 
The frequency of oscillation is then determined by 
. ” 3 
o,” Wy? . \1 ~— 
a) - 2 (m, @,;) !2 (87) 2-0 xa ae 
pe >xX 
This frequency of oscillation is referred to as the /ower __ 
hybrid frequency. 
Now consider the limiting case where the conditions a 
w? > w? and w? > @,/? are satisfied. Eq. (85) then 
a. 
leads to 
@ (ow a Wy") te (88) -l2 
This frequency of oscillation is referred to as the upper = eon ; 
2 f. Bi Fig.5 The behavior of Z (V/c)?, where V’ denotes the phase 
hybrid frequency. . s , ; 
, : e , cag velocity of the wave as a function of x = w/w, for electron oscil- 
For the high frequency oscillations, wm > w,, 


pi 


we may neglect the motion of the ions, and Eq. (84) 
leads to 


lations for propagation transverse to the magnetic field for the case 
when @, ®,, 


9 9 
@- i vy ° 
o+o,2- - | (co? + Wy —s — ke 3.14 Weak Magnetic Field 
: o,2 — or o,2 — a? 
wow, 0 89 We now consider the effect of a weak magnetic field 
-W ; ; i 3 m : , : 
Be w 2 — w? (89) on the high frequency electron oscillations described 
@ ~ S w,2 — 2k Q 0 
a? - 2 - 
es oi \ oO,” : >» WO, 
0 ow c2 ke + 3 _~—~e On 2 2 
L, w2 — oa? ? w,° — w* 0. (82) 
“ « 2 
; S >» WO, ‘ S ee 
mad o,~ F ‘ oe + @,~ ‘ 4 
0 ZZ”? w2—o Foe? wy 2 — ot 
9 = 
2 ca, oe a 
w—cPk+ Yaw? —s » w,* 
Dp 9 9 F p 9 9 
hag @,° —w* yp ,* wr 
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Fig.6 The behavior of Z (V/c)?, where V denotes the phase 


velocity of the wave as a function of x = w/w, for electron oscil- 
lations for propagation transverse to the magnetic field for the 


case when w, = 10, 


in § 2 for arbitrary direction of the magnetic field. The 
weak field approximation, specifically, implies that 
m. Retaining terms up to second order in @,/@, 
the dispersion equation now becomes (cf. Eqs. (61) 
and (63)) Eq. (90) see bottom of this page, where we 
have put 


@, 
é 


Wy? = w,? + c7k; (91) 
@, is the frequency of oscillation of the transverse wave 
in the absence of a magnetic field. Since we have assumed 
,/@ to be small, to a first approximation, for the 
transverse oscillations, @ can be replaced by wm, in all 
elements of the matrix R except R,, and R,,. Then it is 
clear that R,, is the leading term in R, being of order 
c?k* and Eq. (90) gives 


Pee, =9. (92) 


Since & +0, we must require that /, vanish. Thus to 
first approximation in w,/w, the longitudinal and 
transverse modes of oscillation remain uncoupled. 
Retaining terms up to first order in w,/w, Eq. (90) 
reduces to 


o 


«w? — w,2 -i 2 — cos 0 ee 
o 
0 , 
0, (93) 
° 9 We 9 2 
im,2—~cos? ow —wa,2 E 
pe 0 y 
Mp 
w,2 7) 
ww? —w—w,,2 —5 cos 0 iw,,2 — cos0 
we w 
o w,2 
° 9 e 6 
im, — cos 8 w?—a,2—w,,?—G 
70) pe a 
@ 2 @ 
w,,2 — sin® cos0 1 jm, 2—* sin® 
pe a —~ @ 


— 2k? —w,?2 cos?6 
i@,, @, cos 6 


c2 2 


1 Wy, W, Sin 7) 


2 
o, 


w,? sin0 cos@ 


In order that this equation have a non-trivial solution, 
we must demand that the determinant of the co-efficients 
vanish; this leads to the dispersion relation 

@ 


“cos 0. 
0 


2 24 2 
0) Wo" = ye 


= (94) 
These modes of oscillation are circularly polarized with 
Ei, = + 7E,,. It can be verified readily from Eqs. (90) 
and (94) that to first order in w,/m, E, vanishes. 

To carry through the calculations consistent to the 
second order in w,/m, we now substitute @ in accord- 
ance with Eq. (94) in all elements of the matrix R except 
R,, and R,, in Eq. (90). Retaining terms up to second 
order in w,/@, we obtain 
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w,” 
2 2) 2 Ve 2 f)) _ 
o @,* +- ~w,, 5 (1 + cos*@ 
0 2 pe We? ( | ) ame 
2 2 4q@2 
On” MO, P on ge Oh , Be &,;* : 
_— 4 cos? @) I F— — cos8 + —= — cos® g' 
2m \ (Wg? Wo Wo! We f 
° 2 
ee f . @),.. 
+ | 4 (1 —cos 6)? + sin? 6", |2 (1 + cos? 0) 
Vay | 1 Cc k2 I 
0 
2 1 
| | 2 
| pe (Oz 
+ sin? @ —™. . 95 
c? ke | 9) 





In particular, for 0 — 0, we obtain 


2 

9 , ai M, 
pe ” 
My 


9 
' 1 Dye 


>|» (96) 
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ow | + 

pe 2 +7 2 
Mo < Wy 


On" +w 


and for = 2/2, Eq, (95) gives 


1 ow,” w, 2 
WO" + o, 2 “x 1+ {1 + 2/2! . (97) 


ow ) 
? pe 2 
a Mo 





We now find that the longitudinal and transverse 
waves are coupled; their polarizations are given by 


ae 10,20 
E Fill+>5—"5 — cos0 (98) 
“y <= Mo" Mo 
and 
1 @ 
2 . 9 @: 

i= @,* sin ™ 99 

5. c2h2 ?* wy (99) 


Let us now consider longitudinal oscillations which 
are characterized by the plasma frequency w = w,, in the 
absence of a magnetic field. To a first approximation we 
can put wm = @,, in all elements of the matrix R except 
R,, and obtain the secular equation (100) below. 


wo,” . 
w,,2 —S sin@ cos@ E 
pe or 


iw,” We 
-_ t - 
pe w 


sin§ I 


9 
d, 


oe.” 
9 9 9 afi 
OP — Wy_” —Wyq" — g Sin 0 


-iW,,@, CoS) w,? sin# cos 


i@,,@, Cosf 0. 


(100) 


w—w,” 


pe 


o,2 sin?@ 


\¥ 


The electrostatic oscillations are characterized by the 
condition that ck > w. Thus the leading elements of the 
matrix are R,, and R,,, and to compute the corrections 
to the root thus obtained, it can be readily verified that 
it is sufficient to demand that the co-factor of R,, vanish. 
This yields the dispersion relation 


9 
o,* 
P 


2 y 214 | 
w Oy, aa 


o,2 ] 
3 242) sin?6| (101) 
t Doe” 2k 


The correction terms are small as long as m,”/m,,? and 
«,”/c#k? are both much less than unity. 

It may be remarked here that for transverse oscillations 
the correction terms to wm”, to lowest significant order, 
are of order w,/@, while for longitudinal oscillations 
the correction terms, to the lowest significant order, are 
of order (w,/@)?. 


3.15 Strong Magnetic Field 


We shall now consider the case of a strong magnetic 
field which implies that w,/m > 1. Retaining terms up 
to second order in m,/m, Eqs. (61) and (63) lead to (102), 
see bottom of page. In writing Eq. (102) we have used 


ow, 42 NMC ; m 
Zz wo? Be \1 4 M 


4a NM? <3 


Be gz’ (103) 
0 " 
and 
W,- Wy? 
pe pi (104) 
M, @ 
One of the roots of Eq. (102) is clearly given by 
c2 2 
ow (105) 


1 + (c2/.A2)’ 


this is the dispersion formula of Astrém for the extra- 
ordinary hydromagnetic wave. The ci acteristic fre- 
quencies of the other two modes are determined by the 
algebraic equation (assuming ¢/A > 1) 

co cz , 


ow! ——, w? (c?2k2 + e,,?) + 2k w,,? cos? 6 


A A 0. (106) 
The roots of this equation are given by 

oe —~ <2 k2+ w,,” 
-4 2 kw? cos* O]"'2. (107) 


In hydromagnetic waves w 
reduces to 

“4 

\Z 


this is the ordinary hydromagnetic wave of Alfvén. 


ck, and Eq. (107) 


pe 


AP; (108) 


a” —c? k? + ,,2 sin? 6 0 
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4. INSTABILITY IN A PLASMA CARRYING CURRENT 


We now consider the case of a uniform plasma com- 
posed of an equal number of electrons and protons in 
the absence of any external magnetic field. The ions will 
be assumed to be at rest, and the electrons, moving with 
a constant velocity #. We shall neglect the magnetic field 
which results from the electron current, a procedure 
which is plausible if the associated magnetic forces are 
much less than the electric forces considered. For the 
problem at hand, Eqs. (58) and (60) reduce to 


(2 — c? &) E +- c?k (k- E) 
; k- 
4rziwSeNi\v- u x) 0. (109) 
and 
oo. e E+! k XE 110 
’ mo k-u ow” )}- (110) 


We first consider electrostatic oscillations only and 
take k to be in the direction of the electron stream. 
Eqs. (109) and (110) then lead to the dispersion relation 


Ore o> 
(111) 


we (wo — ku)?’ 
where w,; and m,, denote the ion and electron plasma 
frequencies respectively. It is now convenient to intro- 
duce the dimensionless parameters 


o ku 7 112 
" ay ay Po ee ay 
The dispersion relation (111) then takes the form 
/ == I 113 
(x, y) (x y)2 T x2 . ( 5) 


A schematic plot of / (x, 7) as a function of x for 





) 1 is shown Fig. 7 for a hydrogen plasma. If the line 
fF 1 lies above the minimum of the curve, then it is 
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Fig.7 The plot of F (x,y) given by Eq. (113) as a function of 


for ) 1'for a plasma composed of electrons and protons. 
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clear that there exist four real roots of Eq. (113) and 
that all modes of oscillation are oscillatory. However, 
if the line / 1 lies below the minimum of the curve, 
then there are only two real roots and there must exist 
a pair of complex conjugate roots, one of which has 
a positive imaginary part and is, therefore, unstable. 
The minimum of the curve /’ (x, y) is determined by the 
condition that, for a given y, / (x, 7) be stationary, i.e., 
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Fig.8 The behavior of y, given by Eq. (116) as a function of e, 
the ratio of the negative to the positive ion mass. 


ra) 
5 xe P %) = 9. (114) 
This leads to 
€'ls ae . 
KT gpI en e<1. (115) 


The value of y for which this minimum is just tangent 
to the line / = 1 can be readily obtained by substituting 
this value of x into Eq. (113) and solving the resulting 
equation for y. Denoting these minimum values by x, 
and y,, we have 


Jo = [ells + (1 —e')9]" (116) 


The behaviour of y, as a function of € is illustrated in 
Fig. 8. If_y, >_y, Eq. (113) admits four real roots. Thus 
a cold streaming plasma is sfab/e to all perturbations 
having wavelengths smaller than the critical wavelength 
A, given by 


‘ ul ee a 
9) .s "/5 oy 
A. ane \! t5e (117) 
u es ie : 
~ anh (1—3e), fest. (118 


For values of y close to y,, by Taylor’s expansion, 
Eq. (113) can be written as (cf. Eq. (114)) 


(x aa x,)* Fos (x, >, Vo) 0 . (1 1 9) 
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This leads to 
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x= x, + ié'le Re ys (120) 
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it is clear that this expression is valid so long as 


: 
Je—J “ze 's. (121) 
Eq. (120) implies that 
dx - 
Im dy> to asy—y,. (122) 


For values of y close to the origin, and on the assumption 
that x/y <1, an approximation which is legitimate for 
the unstable roots (cf. Eq. (116)), we obtain 
x=tiehy. (123) 

Equations (122) and (123) give the behaviour of x near x, 
and the origin as a function of y. For the continuous 
range of values of y between 0 and y,, Eq. (113) has been 
solved numerically, and the results are illustrated in 
Fig. 9 where the real and imaginary parts of x are plotted 
as a function of ). 

The growth rate of maximum instability can be 
readily computed from Eq. (113) by observing that at 
the maximum 
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From Eq. (111) we obtain 
dx 1 
dy y i 125 
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Fig.9 ‘The behavior of the real and imaginary part of the fre- 
quency of oscillation (x = w/w,)as a function of the wave number 
of the disturbance (& = ,/u) for the unstable modes of oscillation. 
Since Im dx/dy vanishes at the point in question, it 
follows that the denominator occurring in Eq. (125) 
must be real; we take 

—Loaaetts!® geyony*, (126) 
where « is a real positive number of order €'/s, This 
choice of phase will be justified a posteriori. We use an 
asterisk to denote the value of any quantity correspond- 
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ing to the maximum growth rate of instability. Substi- 
tuting Eq. (126) into Eq. (113) and separating real and 
imaginary parts, we obtain 


F=1+ sp7 ) (1 + 20+ 3.0?) (127) 
and 
11+2a 


Soy ka (128) 


To the lowest significant order, Eq. (128) gives 
2 : 
Resi—) . (129) 


Substituting this value of x into Eq. (127), we obtain 
to the lowest significant order 


eo i (130) 


eo a ee 
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The maximum growth rate of instability, therefore cor- 
responds to (for a hydrogen plasma) 
w* 3'l2 [ €\ "ls ] 


Im )* Im =~ | ; 
' a oS 2 18 
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(132) 


Therefore the smallest e-folding time is 18 w,, 1. For 
a plasma with a particle density of 10" per cc, this time 
is of the order 10°-!° seconds. These instabilities grow 
so rapidly that it seems plausible to view them as provid- 
ing a mechanism for converting the organized energy, 
associated with the relative motion of ions and electrons, 
into the fluctuation energy of the oscillations. This 
provides an additional mechanism for the effective 
dissipation of the energy of organized plasma motion 
over and above collisions, and may consequently lead 
to enhanced resistivity of the plasma, compared to the 
corresponding value for a stable plasma. The behaviour 
of the real and imaginary parts of x* as a function of ¢ 
is illustrated in Fig. 10. 
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Fig. 10 The dependence of the maximum growth rate of insta- 
bility as a function of ¢ is shown by the solid curve. The correspond- 
ing real part of the frequency is shown by the broken curve. 
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It may be emphasized here that these results are valid 
only so long as the phase velocities of the waves are 
much less than ion and electron thermal speeds. When 
this problem is examined from the standpoint of kinetic 
theory in part II, it will be seen that the thermal correc- 
tions transform certain stable regions in the (#, w/&) 
plane into unstable regions, albeit with an exponentially 
small instability rate. Moreover, the region of the 
(4, /k) plane, corresponding roughly to the relative 
velocity u comparable to the ion thermal speed, is 
rendered stable. 

The case of two interpenetrating beams of electrons 
(or ions) can be deduced also from the analysis given 
here by setting ¢ 1 and by going over to a reference 
frame in which one of the beams is at rest. In this case 
one finds that the reciprocal of the smallest e-folding time 
is .866 m,, in contrast to (1/18) w,, as obtained earlier. 


5. PLASMA CARRYING CURRENT IN A MAGNETIC FIELD 


We now consider the case of a plasma, in which the 
electrons are moving with a constant velocity u relative 
to the ions, in the presence of an external magnetic field. 
We shall assume that the magnetic field which results 
from the steady state motion of the charges is negligible 
compared to the external constant magnetic field B,. 
This is plausible if the fractional change in the magnetic 
field induced by these external currents is small over 
the wavelength of the perturbation in question. The 
equations for the problem at hand are (cf. Eqs. (58) and 


(60)): 
(w? — c?&*) E + c?k (k- E) 


4ziw > eNv+, e(-v) 0. (133) 
p a \(o k- u) 
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(o k <u) (w,-k) —k(m—k- u)*| (u- E) 


k(u-E)}. (134) 

We shall first restrict ourselves to the case when 
B,, u, and k are along the z-axis. Then w, 
and Eq. (134) reduces to 


i (wo —k-u)w, X E—iw, 


k vanishes 


e7 1 i : 
moO? (w k-u) | [Per (wm ku)* c, 
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where n denotes a unit vector along the z-axis. From 
Eqs. (133) and (135) we obtain for longitudinal oscil 
lations 


a a? | (13 
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and for transverse oscillations 
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where 
E E,.+iE 


“2 ot uv 


(138) 


represents the electric vector for the right and left 
handed circularly polarized waves respectively. The 
upper sign in the denominator in Eq. (137) corresponds 
to £,, while the lower sign corresponds to £_. 

When £, does not vanish, we obtain for longitudinal 
oscillations the dispersion relation 


— ’ 
> oF (w — ku)? 1. 


It is to be observed from this expression that the lon- 
gitudinal oscillations are unaffected by the presence of an 
external magnetic field (cf. §4). Consequently, our 
earlier remarks are applicable in this case as well. For 
propagation along any arbitrary direction the present 
discussion of longitudinal oscillations is also valid if one 
replaces & by &), the component of k along the magnetic 
field. 

We shall now discuss the implications of the dispersion 
formula for the right handed circularly polarized wzve 
only; a similar discussion applies to the left hand 
circularly polarized waves also and will not be given 
here. We have for the dispersion relation (cf. Eq. (137)): 


~ . w—ku 0 
- 2 a,” , 
>? w—ku—o, 


Denoting the stream velocity for the ions and electrons 
by 4, and a, respectively, Eq. (140) takes the form 


(139) 
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-Wy2 {1+ }=0 141 
ve | ' @O— ku; —o, , (141) 


where, it may be recalled, wm, and w, denote the ion and 
eiectron cyclotron frequencies respectively. Introducing 
the non-dimensional parameters 


o ku mM 
x J »€ ’ 
— Doe M | 
a - ’ (142) 
, ¢ ck 
4 omy | 
Oye 4 Doe 
we can write Eq. (141) as 
F (x 1+ e—x® + ao? — & x 
(x) =14 | —— a7 
z 
= 0. (143) 
{+ ¥—Je 


We shall now discuss the implications of the relation (143). 
Let us first consider the case when there is no stream. 
Eq. (143) then reduces to 


1+ €—x? + 9? — e& —*_ — 


ex —x 


F (x) = 0. 


(144) 


The analysis of Eq. (144) is most readily effected by 
means of Fig. 11, where F (x) is plotted schematically 
as a function of x. Observe that for large x, F (x) >—x?, 
F (0) = «?, and there are simple poles at x — — x and 
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x = &€x. Thus the curve F (x) consists of three branches 
as indicated in Fig. 11 and there are clearly four real 
roots; consequently, all modes of oscillation are 
oscillatory and are stable. 
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Fig. 11 A schematic plot of F (x) given by Eq. (144) as a function 
of x. The singularities occur at the clectron and ion cyclotron 
frequencies. The curve F (x) intersects the x-axis at four points 
corresponding to four real roots. 


3.1 Stream Instability (Transverse Oscillations ) 


We shall now consider the transverse oscillations for 
the case where we have two beams of fast electrons 
traversing the plasma. Suppose the fraction 1 — f of the 
electrons is moving with velocity — u,, and the fraction p 
is moving with velocity u,. The net current due to the 
electron streams is 


J=eN| 


In order that the net current due to the electron streams 
vanish, we must demand that 


(1 — A) u, + Bua). (145) 


(1 — Pp) 4, = Pug. (146) 


The corresponding results for the problem of one beam 
can be obtained immediately from these results by 
setting #, and # equal to zero. For the beam of rela- 
tivistic electrons, the only modification in the dispersion 
relation is to replace the electron mass ™ by 


My 
[1 — Way 


where ™, is the rest mass of the electrons. The dispersion 
relation now becomes 


m YMo , 


(147) 
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—£ pe E : | Ve -0, (148) 
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where w,, now denotes the electron plasma frequency 
defined with respect to the electron rest mass, and, it will 
be recalled, that k and u are both along B,. Introducing 
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the dimensionless variable given in Eq. (142), we can 
write Eq. (148) in the form 


} ) 1 , 1 p > > x 
x - é x2 \2 —- <2 
\ \ P) Y1 V2 Ex —x 
l X/V1 xB | y-* 
(1 5) - = —7— 0. 
P) mx +t) = *+Je—(X/ V2 
(149) 


Firstly, we observe that if we have only one beam of 
electrons, Eq. (149) reduces to 


F (x) = 1+ e—x*+4+ a? — & =* = 
x 0. (150) 
yex—nt+kin) ° ™ 


It turns out that Eq. (150) has real roots, corresponding 
to oscillatory motions, for all values of the parameters 
of physical interest. 

We now consider the case of two beams of electrons. 
The first point which is to be observed is that the effect 
of adding another beam is to introduce a new pole as 
shown schematically in Fig. 12. As 4, increases from its 
zero value, the additional root moves to the left in Fig. 12. 
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Fig. 12 A schematic plot of F (x) given by Eq. (150) as a function 
of x. Note the additional branch of the curve introduced by the 
presence of interpenetrating beams of fast electrons. The curve F (>) 


intersects the x-axis five times corresponding to all real roots. 


When 4, is so large that the gap introduced in the curve 
of Fig. 12 straddles the x-axis, two real roots are lost 
and a pair of complex conjugate roots are acquired. 
One of these must have a positive imaginary part and 
corresponds to instability. This is illustrated in Fig. 13. 

When instability does occur it has a simple physical 
interpretation; namely, the condition that the gap 
straddles the x-axis is that 


x + Je + (2/Y2) = 0. (151) 


This can be interpreted as saying that, in its rest frame, 
the electron beam sees the circularly polarized wave as 
having a frequency equal to the electron cyclotron 
frequency, and the in-phase rotating component of the 
electric field is resonant with the particle motion. Since 
instability occurs only when the gap straddles the x-axis, 
it is clear that for this to occur the new pole must lie 
between the electron and ion cyclotron frequencies, 
i.e., between — z and &z. 
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Fig. 13 A schematic plot of F (x) given by Eq. (150) as a function 
of x for the case where the interpenetrating beams have velocitics 
corresponding to unstable modes of oscillation. Observe that J° (2) 
now intersects the x-axis only three times, corresponding to three 


real roots and two complex roots. 


}.11 Oscillations much below the Ion Cyclotron Frequenc) 


We now consider the consequences on the hydro- 
magnetic waves given by Eq. (78) of the introduction 
of beams of fast electrons such that the small gap straddles 
the real axis near the unperturbed root (denoted by 
subscript 1): 


1 "la 
m 59 
xy v(1 ea] : (152) 
Then for the case contemplated 
. x \-* 1 
F (x) = € —x*® + a? —e (1 — | —~(1— A) 
EX Y1 
1 1 x+y) ; 
— B—(1— p)—|1 —=-4) 
#2 1 4/71 
P : 0 (153) 
Yo" x Yo + (2/72) ei 
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x x 
F (2) - x2 + a? —| =" 
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- 0, (155) 

V2 + (% 2) | 

where in writing Eq. (155) we have used the fact that 
(1 — p) 1 = Pde (156) 


Making use of Eq. (152), we can write Eq. (155) as 
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For the case under consideration x/&z < 1 and, further, 
since we are interested in the case when x ~ x 


» we 
may set 
x42 — x2 ~ — 2x, (x — x). (158) 
Eq. (151) then reduces to (assuming x/zx < 1): 
l Vo 
F (x) -| 1 ai 2c, (xe — 2¢,) ++ BI = a 
A 72 a 
pz : 
0 159 
V2? [x —Jo+ (z/72)] = 
Letting 
p (1 <5} 22y 
fd Yo 
q P\, ai |, (160) 
72 v 
' p 
ee ee 
‘2 





Eg. (159) then leads to the following quadratic equation 
for x: 





' _ % q 
x2 + x (—x,-4 Yo— 
"Ve Pp | 
4 z r 
(+ 2 \(2—*)—— = 0. (161) 
' p \ ,) p 
The roots of this equation are given by 
7. _ q 
x r x, — Jet 
2 (sy, tet 5] 
] q 4r]'ls 
+ x +2—re+ =P | - (162) 
1 2 | 1 p 72 =| p 


Explicitly, we may write this as 


macanl 1 Lee M 4.12 [(1/ 42) — (Aaq/@?)| w,,, 

21 yn FT PO wif + (42 NM2/B,2)] 

1 {f mo , | 4 [(1/%2)—(a/@<)] o,, \2 

£5 {fo + St — kent BaF ae NAB 
2 Bw,,2m,/Y2" a P 


| [1 + (42 NMe?/By")j (163) 


It is clear from this expression that the maximum 
instability rate occurs when the term in the braces 
vanishes; this is indeed the approximate criterion given 
in Eq. (151). The spread in velocity of the electron beam 
about that given in Eq. (151) for which there is instability 
is given by (cf. Eq. (163)): 

2Pm,,2@,/)'2" 


el #2 (164) 


kA = Olli + 4aNMA/BQ)] 


5.12 Oscillations near the Ton Cyclotron I ‘requency 


Let us now consider the case when the physical 
parameters are such that in the absence of the beam the 
interesting root is given by Eq. (81) namely, near 
the ion cyclotron frequency. For the case under con- 
sideration, x < x, and we shall assume that J, is of 
order x. Denoting the root in question by the subscript 2, 
we have (cf. Eq. (81)) 


(165) 
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Equation (149) can then be written in the form 
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Eq. (167) leads to the quadratic 


i > 1 e \ ¢ 
4 x(—m— + £8) 4 (yp —2)(4 4 2) 
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Ay 
The roots of this equation are given by 
] 2 Gy 
os 7. + %3+Is—— +f 
ba 72 Pi 
4 og 2 4r.)2 a 
+5 |\%—rn+ — +4) -—+). (170) 
ad 72 Pl Pp 
The maximum growth rate of instability occurs when 
> ¢ " 
X_—Jet = ,f1 Uv; (171) 
72 Pi 


which is again similar to the approximate criterion 
stated in Eq. (151). The spread in velocity about that 
given by Eq. (151) for which instability does occur 
is given by 

: eo lee—oo 
4r, 4 Px (€x—X3) 


a V2 E 


kA ty (172) 

The instabilities due to the presence of streams of 
electrons (or ions) may be responsible for microwave 
noise in the laboratory plasmas and the radio emission 
from astronomical objects. 
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We shall present here a summary of the dispersion 
relations for the various modes of oscillation of a plasma. 
a) Static Plasma in the Absence of a Magnetic Field 

(i) Longitudinal electron oscillations (w? > &? O/m) 


6. SUMMARY OF SEcTIONs III2 


9 9 9 
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(173) 


(ii) Transverse electron oscillations (w? > £2 O/m) 


9 9 .2 k2 L? 0 ~~" i 1 74 
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(iii) lon oscillations (O; < 0,) 
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b) Static Plasma in a Magnetic Field, k || B, 
(iv) High frequency oscillations wm > o,. 
o,” < On. cz 2 
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(176) 
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(v) Oscillations near electron cyclotron frequency 
2 w.2 + c2f2 
(ao? > o, ck?) 


ne? lari 
o o,\1 o2 —(2# 4-o,,2) (177) 
e , 


(vi) Oscillations much below the ion cyclotron 


frequency 
2 L2 
we i= (178) 
€ 4 
1 + (c?/.A?) 
(vii) Oscillations near the ion cyclotron frequency 
(wo? > o,/ 7k?) 
9 
1 Opi (179) 
@ @ - F eure ‘ / 
: w? — (c?k2 +- w,,”) 


c) Static Plasma in a Magnetic Field, k | B 
(viii) Lower hybrid frequency (w < ck), 


OD; A. WD We \ Dre 


@ = (w, @,)'!: (180) 
(ix) Upper hybrid frequency, w > @,,o > ,, 
w= (w? + ,,?)'!2 (181) 


d) Stream Instabilities 
1. Longitudinal oscillations, in the presence of a 
stream of electrons with constant velocity u, are 
unstable to all perturbations (along u) having 
wavelengths greater than 


u 3 4 1 


L 29 
@,, \ 2M (182) 


Smallest e-folding time = 18 w,, 


This instability is unaffected by the presence of 


a magnetic field along u. 

Transverse oscillations, in the presence of two 
streams of relativistic electrons, are unstable near 
the Alfvén frequency and the ion cyclotron 
frequency, k and u being assumed to be along B. 
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7. HYDROMAGNETIC WAVES FROM CHEW, GOLDBERGER, 
AND Low THEORY 


We have shown earlier that, for a low density plasma 
in a strong magnetic field, the kinetic equations for the 
plasma reduce to equations for a single fluid if we 
neglect the heat flow along the lines of force and/or 
assume the low temperature approximation. We shall 
now discuss the propagation of hydromagnetic waves 
using this formalism. 

Consider, then, a spatially homogeneous plasma in 
static equilibrium in a uniform magnetic field, which we 
shall assume to be along the z-axis of a Cartesian system 
of co-ordinates, with the pressure tensor given by 


P=p,I1+(p —/,)nn, 


where n denotes a unit vector along the field lines. In the 
perturbed state, let the various quantities be denoted by 


v,0+00,p,+0),,p +0p,B+ 0B, (184) 


(183) 
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The equations governing the perturbed quantities are 
obtained by linearizing Eqs. (I1-135)—(II-140). These are 


ov 1 


0s, —-~V-dP4 aa < 0B) x B, (185) 
57? 0 oV-v, (186) 
") 
5B Vv x (v xB), (187) 
5 fa) B 00 
¢ hy lL p = RS 
p Pi Rte >: (188) 
and 
, , OB ; 00 . 
p <p Bt 3? °° (189) 


where 6B B- 6B/B. Introducing the Lagrangian 
displacement &,v 0€/07, Eqs. (186) and (187) can 
be readily integrated to give 


do =—oV-E (190) 
and 


6B=V x (E x B), (191) 


where the constant of integration, by definition, has 
been set equal to zero. It follows readily from Eq. (191) 
that 


dB V,-&,, (192) 


where we use a subscript | to denote the component of 
a quantity in the plane transverse to B. Thus 


V,=V 


The component of a quantity along B will be denoted 
by the subscript ||. We shall assume that all the perturbed 
quantities behave like exp / (— w¢ + k-r). Using the 
foregoing results, Eqs. (188) and (189) reduce to 


nn-V. (193) 


Op, pi (2V,-b,+ék &) (194) 
and 
OVp p (V,-€,+ 31k). (195) 
From the form of P, it follows that 
V-6P=V,6p,+7(pj)—p, FA On 
+- n(p -p,)V-dn ik Op. (196) 


The change dn in the unit vector along the lines of force 
is obtained readily from Eq. (191): 


on=—n-VE,+|(n-V)nn)-&. (197) 
Using the relation (197) into (196), we obtain 
V-OP = V, 6p, + (p,—p ) 8, 
+ik |(p, —p,)V,-b, +p jn. (198) 
From Eq. (191) it follows readily that 
(V x OB) x B= B(V,V,-€,— 47). (199) 


Using the foregoing results the equation of motion 
(185) splits up into the two equations 
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and 


(ow? — 3p k*)z k p 'k,-&). (201) 


On substituting for &, in accordance with Eq. (201) 
into (200) and forming the scalar product with k,, we 
obtain the dispersion relation 


. R2 r BB 
Ow k | p -p; ta) — 24:7 | A.+ =| 
L 2 
a (202) 
* 0" 3p k? 


If @ denotes the angle which the propagation vector 
makes with the z-axis, 


k =kcosO, k,=ksin0. (203) 


After some simplifications, Eq. (202) reduces to the 
following quadratic equation for q?: 

B 
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1 peal 1 + sin20 29) (204 
‘ae p cos is Pp, (1 + sin?) —p cos (204) 
kA 
re p,sin?0 cos?0 —0. 
The roots of this equation are given by 
. k2 ; B2 " _ 
w* 20 ca: +p, + 2p cos? + p, sin? | 
k2 {| B Sf af 2 _ 
 20\\|42 7 p, (1+ sin? 0) —4p cos?0 (205) 
+ 4p,2 sin®0 cos?) 
For propagation along the lines of force, 0 = 0 and 


Eq. (205) leads to the following two modes of oscillation 
(distinguished by the subscripts 1 and 2): 


ke) B 
mH (tn) ew 
and 
3 
wg =~ p (207) 


It is clear that the mode of oscillation corresponding to 
Eq. (206) is unstable if 

pi >aa th: (208) 
That is, if the distribution function for the ions and/or 
electrons is strongly peaked along the lines of force, 
the hydromagnetic waves propagating along the field 
lines are unstable. 


For propagation transverse to the field lines, # 2/2 
and we obtain 
> £2 / B2 
9 = { 
we? + p,). 209 
o | 82 Pi \ 


The corresponding phase velocity of the wave is given by 


@ 2P “le 
"> k | 0 ; (210) 
where P — p, + (B*/82) denotes the total transverse 
pressure. 
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For any arbitrary direction of propagation there exist 
unstable modes of oscillation for all @’s less than 6, as 
determined by the equation 


B - 9 p - 9 9 \ 
sat? (1 + sin? 0.) 3p sin? 0.+ p cos?@,. (211) 
As an example, if we take p, = B?/82 and p B?/22, 


we find from Eg. (211) that 6, — 26° 45’. 

It will be shown in part II of this survey that in the 
presence of a magnetic field, whenever the distribution 
function is anisotropic (i.e., P < P, or P, < P), we 
have instabilities, although their e-folding times tend to 
be rather small. 


8. HyDROMAGNETIC WAVES: GENERAL THEORY 


In our previous simplified discussion of hydro- 
magnetic waves we neglected the effects of the finite 
temperature, viscosity, electrical conductivity, thermal 
conductivity, and other transport co-efficients of the 
plasma. These shall now be taken into account in order 
to exhibit the coupling of the hydromagnetic and sound 
waves, and to discuss the damping of these waves due 
to dissipative effects. However, for the sake of sim- 
plicity, we shall assume that the ion gyration frequency 
is much smaller than the effective collision frequency 
so that all the transport co-efficients may be taken to be 
isotropic. Eqs. (II-94)—(II-97) then form the basis of 
our discussion. 

In the equilibrium state U—0O and B I 
where B, is constant and p and @ are uniform. For small 
departures from the equilibrium state let the various 
quantities be denoted by 


ViP+rA OTA, 


The equations governing the small departures are 
obtained by linearizing the equations of motion 
(11-94)—(II-97). These are: 


and B, ++ B,. (212) 


00, oV-v 0, (213) 
oy (214) 
I _— I 
-Vpy+ 4x (V x B,) x By + 4V?v 3 uVV-v, 
OB, ; 5 aa 
ry, V x (v x By) + in, ° B,, (215) 
and 
) 00 yA e 
“fs s2 sf +3 V-(4,V0,+»,D,), (216) 


where s = (yp/0)'': denotes the sound speed and D, is 
obtained from the linearized form of Eq. (II-89). Thus 
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0; 


Vv x B,—«V@,). (217) 
The temperature and the pressure are connected by the 
relation O = Mp/o; the linearized form of this gives 


M 
0, - (218) 
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Introducing &, the Lagrangian displacement of an 
element of the fluid, and assuming that all quantities 
behave like exp 7 (— wt k-r), after some straight- 
forward reductions Eqs. (213), (215), and (216) lead to 


0 oV -&, (219) 
B,=aV x & x B,). (220) 
and 
pi = 90, (221) 
where 
1 ik2¢2 
+ (222) 
x 470, 
and 
sie 
1+ 3/ A, —(¥', %,/0,) || A? M/(y wo) 
q A (223) 
1+ 3/ A, — (0, %,/0,) || M/(o w) 


It must be emphasized here that in writing the foregoing 
expressions we have assumed the various transport 
co-efficients to be isotropic. However, this may not be 
true, in general, for oscillating electromagnetic fields. 

Using the foregoing results in the equation of motion 
(214), we obtain 


” = 2 . a .. 
| w? + im = NE {sq 3/0 ° \kk -& 


+ aA? {kx [k x (§ x n)}} xX n, = (224) 


where n is a unit vector along B, and A is the usual 
\lfvén speed. Assuming that k lies in the (xz) plane 
making an angle 0 with the z-axis, Eq. (224) leads to the 
following equations for the different components of &: 








iu = 
| oo? x A? k2 cos? 9 + ; kw), 0, (225) 
) 
| u iu cad 
or — x Ark +i Ba— (2% op— f wk sin? 6|&, 
0 30 
Ou ee t 
\g 2 j—" lk sinfcosOé, —0, (226) 
30) “ 
and 
imu ‘ ; 
| ps? 3 - | & sin 9 cos 0 &, 
i wu iu 
jo + wk? —\q $ — 3,0 |AcostO £,=0. (227) 
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8.1 Absence of dissipation 


Let us first consider the limiting case of an ideal 
hydromagnetic fluid for which all the transport co-effi- 


cients vanish. In this limiting situation 1 = 1, ¢ 5, 
and « = 0. Eq. (225) then leads to 
@=+Ak cos 9; (228) 


these modes of oscillation remain unaffected by the finite 
plasma pressure. For Eqs. (226) and (227) to have a 
non-trivial solution, we must have (setting « 1, 
g lu 0): 


? —A?2 k — 52 k2 sin? 6 


—s2k2 sinf cos ow 


— 52 k2 sin cos 
— 52 k2 cos? 
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This leads to the following quadratic for m?: 


oo — a)? (A® + 5°) & 252k cos? 4—0. (230) 
The roots of this equation are given by 
i vy 
ms ( {2 1. 52) £2 
1 45? _A? cos? 0\' *| 
k2 (A? +- 5?) (1 23 
7 * ( 5) | (2 1442 jp (231) 
In the case of a strong magnetic field, s/A < 1, and 


Eq. (231) leads to the following roots (distinguished by 
the subscripts 1 and 2): 


@,? = k? (A? + 5? sin? 0) (232) 

and 
Me” 2 52 cos? f). (233) 
In the case of a weak magnetic field, A/s < 1, and 

Eq. (231) leads to 

@,? = k* (s? +A? sin? 6) (234) 

and 
oe k2A2 cos? A. (235) 


The principal result which stems from including the 
finite plasma pressure is that the hydromagnetic waves 
remain unaffected only if the motion of the material 
(i.e., &) is normal to the plane containing the magnetic 
field and the propagation vector; otherwise the disturb- 
ance breaks up into coupled hydromagnetic waves and 
sound waves. The chief application of this result is to 
waves originating in the surface layers of a sunspot 
where A may be comparable to s. These remarks may 
be of some significance also in laboratory plasma 
diagnostics. 

That the a? are real in the case of ideal hydromagnetics 
also follows from the self adjointness of the equations 
of motion which may, in general, be written in the form 


wob =F (E) - 


where F is a linear operator acting on &. The self adjoint- 
ness of F immediately leads to the conclusion that w” 
is real. Consequently, the criterion for stability is that 
there be no negative characteristic value @*. Also, it can 
be shown easily that a necessary and sufficient condition 
for this circumstance is that the potential energy 


(236) 


- \ dx E - F (&) (236a) 
always be positive. This latter condition is a statement 
of the energy principle governing the stability of ideal 
hydromagnetic system. This criterion, while derived on 
the assumption of collision-dominated plasma, can be 
demonstrated to give a lower bound on stability as 
compared with the conclusions of the theories based on 
the collisionless Boltzmann equation in the limit of a 
strong magnetic field. 


8.2 Effects of dissipation 


We now consider the effects of dissipation on the 
various modes of oscillation discussed in the preceding 
section. Eq. (225) gives 
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c2 ke 1 Zu 


Lg " k2 a, 237 
' 470, 0 ( ) 


wo — A? k2 cos? 041 
If c24?/4210,@ and &/0 A? are both much less than unity, 
to the lowest significant order Eq. (237) yields 

2 


é a7 
470, r 0 
We now consider Eqs. (226) and (227). In order that 
these admit a non-trivial solution, Eq. (239) see bottom 
of this page. 
On the assumption that ?k?/420,w, MA*/Om, and 
[A, — (%4%,/0,)] [42/9] are all much less than unity, 
we obtain in the strong field case (s/A < 1): 
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(238) 
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and 


w? = k? {stcos*d -ikRA E (1 + cos?) 
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Under the same approximations for the weak field 
case we obtain 
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D) 6 = 2 . @ 
w,? = £5* ! + —, sin? 0 


mi4un . 2H [, _%4\My ase 
Ss E 0 3 rate "1g, | 0 y | (242) 


and 
w,2 2 { A* cos? +E sk). (243) 


It is clear from Eqs. (240) and (241) that the damping 
of the hydromagnetic waves is 


ke 
expi— 5 


Since s/A < 1, we may write this approximately as 
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Baal be 7 
exp - ia Ale (245) 


This suggests that only the resistivity and the viscosity 
of the medium are significant for the dissipation of 
hydromagnetic waves. In order that these waves be 
realized in any physical situation, it is essential that the 
time required for their amplitudes to be reduced by 
a factor of e~! of their initial amplitudes be much larger 
than one period of oscillation of the wave. Thus we 
must have 

2|{¢ a : 

k \420, 7” 0 | kA’ (246) 


Since & 2a2/A, this relation implies that the wavelength 
P g 


. iu a 
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| 0 30 
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of the disturbance must be greater than the critical 
wavelength A, given by 
9 
i ce bh 
+ =}, 247 

P A \tze, 4 (247) 
IV. Oscillations of Bounded and 

Non-Uniform Plasmas 


1. OSCILLATIONS OF A CYLINDRICAL PLASMA 

We shall now investigate the oscillations of a cylin- 
drical plasma immersed in a uniform external magnetic 
field which will be assumed to be along the z-axis. 
In the equilibrium state the plasma is at rest and the 
electron charge density is neutralized by the ion charge 
density. We shall assume that the plasma is so dilute 
that direct interparticle collisions are negligible and that 
particles interact through long range electromagnetic 
interactions only. For the sake of simplicity we may, 
therefore, neglect the plasma pressure. For small depar- 
tures from the state of static equilibrium the equations 
of motion for the electrons and ions, assumed singly 
charged, are 


ov, 


Mm ¢ \E - 7 Ve > B, ) (1) 

Ov; 1 
M~ e(E4 ~¥; X By), (2) 

together with Maxwell equations 

ss 4% i 1dB 
VxB , J and VxE = ae? (3) 
where 

J=eN(v;—v,)- (4) 


We have neglected the displacement current as we are 
interested only in low frequency oscillations satisfying 
the condition w < ck. From the definitions of the mass 
density and the mean mass velocity, we obtain 


1 oj 1 1 aif 1 
Ne d+ | MT m |E V¢ | Mit v.| x By 
1 > * 
oe \E i By}, (5) 
and 
ov 1 
CB: ; J ya B,, (6) 


neglecting terms of order m/M. We now need a relation 
connecting v,, v, and J. From the definitions of J and v 
we readily obtain 

Ve~ Vv 


ne): (7) 


From Eqs. (5) and (7) we get 


9 dj | l l > 
42w,* ry E + ~ vx B, — = ved B,- (8) 
\¢ sf — a | A sin 0 cos 0 
. vt 0. (239) 
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\s we are interested in frequencies of oscillation much 
below the electron cyclotron frequency, we may drop 
the term on the left hand side of Eq. (8) and obtain for 
Ohm’s law: 
] ] 
E+-—vxB -Jx B=0. 9 
” cNe J . ( 


; 
From Egs. (3) and (9) we then obtain 


oB 


1 
op = VX \v x Bo—7,J X B, }. (10) 


Equations (6) and (10) form the basis of the charac- 
teristic value problem which must be solved subject to 
appropriate boundary conditions. We now assume that 
the ¢, 0, and z dependence of all quantities has the form 
exp 7 (— wt + mO + &z). Equations (6) and (10) then 


lead to 
“B= V x {UV x B) xn} xn +i2-(V xB) xa 
(11) 


where A denotes the Alfvén velocity and n is a unit 
vector along B,. On carrying out the normal mode 
analysis, we obtain the following equations for the 
various components of B: 


m_. ww w \dB, 
Bec 4 | ge L2 (2 z 
B, Aik| - 2, B, 4 (4 5 — 28) — |. 
(12) 
mk wo w dB 
Lo 2 L¢ z 
Bo A| — (A®&——a — D2) B,— kQ | (13) 
and 
1d 4dB,,({. ™) 
rdr’ dr * (7° ~ b,=0, (14) 
where 
1 w? \2 
| £2 _ #O2 5 
a4=(@—Za) —#2, (15) 
fa~Se— ee (16) 
A 2 — (0?/A) 2D 


and (2 = w/m,;. The solution of Eq. (14) which is 
regular on the axis of the cylinder is given by 


B, By] (zr), (17) 


where B, is an arbitrary multiplicative constant. Having 
determined B,, B,, and By are readily obtained from 
Eqs. (12) and (13). This determines the perturbation 
in the magnetic field inside the plasma. 


1.1 The Perturbation in the Vacuum Field 


The vacuum magnetic field satisfies the equations 


B=Vy and V2y= 0. (18) 


The solution of Eq. (18) which is bounded at infinity is 
given by 


Y = Ba Ky, (Ar), (19) 


where B, is another constant, K,,, (x) is the Bessel func- 
tion of purely imaginary argument (the notation is that 
of Watson), and we have suppressed the factor 
exp / (w@ + &z). 
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1.2 The Boundary Conditions 


The boundary conditions which must be satisfied at 
the perturbed boundary of the plasma are 


A(B-n)=0, (20) 
and 
B? 
A\=—) 0, (21) 
where A(X’) denotes the jump experienced by the 


quantity XX at the interface. The linearized forms of 
these equations are 


BL? = B. (r= R) (22) 
and 


Bry =Be (r=R), (23) 


where have now used the superscripts p and » to indicate 
the quantities referring to the plasma and vacuum 
respectively. Using these boundary conditions, and after 
some straightforward reductions, we obtain the dis- 
persion relation 


a : ow" TT ’, (% R) 
2 Q |. | &2 SPR?) te 
k A [m2 A2 \é A2 e A } (7) 
ARK’, (RR) 
; , (24) 
Re (AR) 


For a specified &, A, and, £2 Eqs. (15) and (16) can 
be solved to give m in terms of 7. When the @ thus 
obtained is substituted into Eq. (24), it admits an infinite 
number of discrete roots for y for prescribed values 
of k, A, and 22; denote these roots by y,. The disper 
sion relation then takes the form 


oo ow , we 1 ). 


wo? (1 ~ £2 _Az} \ A2 £2 4 ne (25) 


1.21 Low frequency oscillations 


Consider first the case when w/w; < 1. Then Eq. (25) 
yields the roots 
w? = k? A? (26) 
and 
w* = A? (k? + x,°), (27) 


which correspond to the ordinary and the extra-ordinary 
hydromagnetic waves. 


1. 22 Oscillations near the ion cyclotron frequency 


We now consider the case when w/AA < 1. Eq. (25) 
then yields 
wo? wo? 


k2 At Az (ke 73)" 


w = w?}1 (28) 
These oscillations are referred to as ion cyclotron oscilla- 
tions. 

In ion cyclotron oscillations the ion velocities are 
circularly polarized and the sense of rotation is the same 
as that for a free ion moving in a magnetic field. The 
electric field is elliptically polarized. One may decompose 
this electric field into two circularly polarized compo- 
nents, one of which will have the same sense of rotation 
as the ion motion in the wave. However, the ion velocity 
is 90° out of phase with that component of the electric 
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SPECTRA 


FOR ION CYCLOTRON WAVES 


growing spiral orbits and there is 
strong power absorption by the ions 
from the electric field. There are 
no waves at this very low density 
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condition all of these modes occur 
at magnetic fields considerably stron- 
ger than the cyclotron field B,. Once 
the ion cyclotron wave has been set 


STRONG POWER 
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1ONS FROM THE 
ELECTRIC FIELD 


HIGHER DENSITIES LOWER DENSITIES up, the ions do not absorb energy 

— | LONGER AXIAL WAVE-LENGTHS SHORTER AXIAL WAVE-LENGTHS |‘— from the electric fields. Therefore, 
; ‘ there is no damping of the wave. 

3 The second picture in Fig. 14 
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Fig. 14 


field which rotates in the same sense as the ion. Con- 
sequently, there is no absorption of energy by the ions 
from the waves. 

However, in the foregoing discussion we have not 
taken into account the thermal motion of the ions, 
and/or we have neglected also the collisions (short range) 
between the ions and electrons or the like particles. 
In either case damping of the ion cyclotron waves will 
result; in the former case we have Landau damping and 
in the latter case we have the usual collisional damping. 
When damping occurs there appears a component of the 
elliptically polarized electric field which is in phase with 
the ion motions. This leads to absorption of energy by 
the ions from the electric field of the wave. Conse- 
quently, this phenomenon is of fundamental interest 
in connection with the heating of the plasma. 

In the experimental situation one has the plasma in 
a magnetic field and one superposes on this a transverse 
electric field which oscillates at the ion cyclotron 
frequency. As the ion motion is in phase with the 
electric field, energy is fed into ion motion, and the 
ions will move in a spiral of increasing radius. 

At higher densities, say 10! ions/cm’, this single 
particle picture is no longer valid. Ion currents induce 
an electric field which has a decelerating effect on the 
ions, and heating a moderately dense plasma at the ion 
cyclotron frequency becomes rather inefficient. However, 
if we keep the frequency of the oscillating electric field 
fixed, but go to slightly stronger magnetic field, we find 
the existence of natural oscillations of the plasma. These 
natural oscillations, which we call “on cyclotron waves,” 
are, in fact, the short-wave, low-density limit for the 
extra-ordinary hydromagnetic waves described by 
Alfvén and Astrém. 

Fig. 14 shows the spectra for ion cyclotron waves. 
The resonances are indicated as lines and are plotted 
against magnetic field strength. We first look at the 
picture denoted by “no waves.” Here we have low 
densities and, as a consequence, single particle accelera- 
tions of the ions in the electric field. The ions move in 
32 


tron waves for somewhat lower den- 
sities. The various modes now occur at magnetic field 
values which are somewhat closer to B,. There is still no 
damping of the waves. In the third picture the density is 
quite low and the modes now occur at magnet ic fields very 
close to By. This is an intermediate case between the 
ion cyclotron wave resonance and the single particle 
resonance. For this condition we have waves, but the 
ions begin to move in growing spiral orbits and pick 
up energy from the electric field. In this case we have 
ion cyclotron waves which are slightly damped. 

The transition from the undamped waves to the slightly 
damped waves is particularly interesting because it 
provides a method for thermalizing ion cyclotron 
waves, that is, for transforming the wave energy quickly 
into heat. It turns out that it is equally valid to talk of 
wavelengths instead of densities, and it is the short 
wavelengths which are damped. To illustrate a plasma 
heating scheme, we make an analogy to ocean waves 
running up on a beach. A cross section of an ocean 
beach is shown in Fig. 15. The waves are moving in 
toward the shore and, because the water is getting more 
and more shallow, the wave length of the ocean waves 
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becomes shorter and shorter. Finally, in the shallow 
water the waves are unable to propagate and they 
“‘break;” the wave energy is transformed into heat. 
Similarly, one might use an induction coil to generate 
rather long wavelength ion cyclotron waves and allow 
these waves to propagate along the magnetic lines of 
force through a region where the magnetic field tapers 
off and becomes somewhat weaker. The wavelength 
becomes shorter and shorter and, finally, when the 
wavelength is sufficiently short the waves damp out and 
the wave energy goes into heating the ions. It woud be 
in this “magnetic beach” region that the real heating of 
the gas would take place. 


“< 


2. RADIATION BY PLASMA OSCILLATION 


We now wish to exhibit the possibility of radiation 
by plasma oscillations. Firstly, it may be remarked here 
that this problem is mathematically quite analogous 
to the scattering of electromagnetic radiation by a 
dielectric. Consider, then, a simple equilibrium config- 
uration in which the ions are held fixed with a uniform 
density and the electrons, whose thermal motions 
will be neglected, are distributed likewise inside a 
sphere of radius RK. From the linearized equation of 


motion 
02 
g eE 


x 
m 52 (29) 


and the linearized Maxwell equations, we readily obtain 


eVxVxE (cw? w,*) E, (30) 


where we have assumed the time dependence to be of 
the form exp (—/w/), and wm, denotes the electron 
plasma frequency. From Eq. (30) it follows that 


V -| (w W,*) E 0. (30a) 


Inside the sphere @, is constant; consequently, if 
w? +@,? Eq. (30a) leads to 


V-E=0 (31) 


for both inside and outside the plasma. 

Because of the symmetry of the problem, it is con- 
venient to seek a solution of Eq. (30) in terms of vector 
spherical harmonics: 


A=#,Y," (0,9,), (32) 
y 
rm (6 
B [A+ Ae (9, 9), (33) 
and 
] 
“m (0 ¢ 
c Le Dye? r Y," (0,¢)|, (34) 


where )"”, is the normalized spherical harmonic. A, B, 
and C constitute a complete set of basic orthogonal 
vectors on the surface of a sphere; namely, 


ffPm;9,9)-Q(m;6,9)d2=0, (35) 


where P and Q are any two distinct vectors A, B, and C. 
The normalization condition is 
(P43 9,¢~)-P (2, m3 9,g) d2=—0,, 6 (36) 
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The spherical harmonics satisfy the equation 


: Ba ; 
V2 Y," (9,¢) 2 Ki +- 1) Y," (0,9). (37) 
We now write 
E 2 A + - B-4 - a (38) 


where Ff, G, and H are functions of r alone. On substitut- 
ing Eq. (38) into Eq. (31) we obtain 


df 


y 


ar 


1 


Ki 1) | ‘G. (39) 





On substituting Eq. (38) into Eq. (30) and making use 
of the relation (35), we obtain the two independent 
equations 





12H I 1 

: jy2 ~» (w? — w,?) 5 l(/-A | H1i=—0 (40) 

dr* c2 73 

and 

dG ] [1 i - 1 i 

dr * (Ki+1)}'h [2 w,*) ya! (/ 1)| / 0 
(41) 


The problem, then, is to solve Eqs. (39)—(41) subject 
to appropriate boundary conditions. 


2.1 The boundary conditions 
From Eq. (30a) it follows upon integration across 
the interface that we must have 


(wo? — w,”) E,-n =o E,-n, (42) 


where the subscripts 1 and 2 refer to the plasma and 
vacuum sides respectively and n denotes the unit outward 
normal to the surface. Since only A has a radial com- 
ponent, it follows from Eq. (42) that 


(w? — w,?) F, 
It follows from ¢cV = E 0B/0¢ that the tangen- 


tial components of E are continuous at the interface. 
This requires that 


wow, atr=R. (43) 


G =G, r= A, (44) 


and 
M,=-f, atr= #. (45) 
Integrating Eq. (40) across the interface, we get 


dH, dH, 


, qr te R. (46) 

The conditions (45) and (46) may be combined to read 
d d a 

7, oA — In Hg. (47) 


2.2 Solution of the problem 


The solution of Eq. (40), which is regular at the center 
of the sphere, is given by 


H=H,r" J. (ns (48) 


where 


1 
kj : 3 (a w,”) (49) 
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and H/, is an arbitrary multiplicative constant. It can be 
seen readily from Eqs. (39) and (41) that F satisfies the 
same equation as 7; consequently, a similar solution 
holds for /. G is then determined by Eq. (39). 


2.3 The case of a perfectly conducting boundary 

An interesting situation which may approximate 
cavity resonators is the one where the plasma is enclosed 
by perfectly reflecting walls. The boundary conditions 
to be imposed in this case are that Eg and &, vanish 
at r = R. This requires that 


Tiss), (A R) = 0. 


Let the roots of Eq. (50) be denoted by /,,,. From Eqs. 
(49) and (50) we than obtain the dispersion relation 


(50) 


2 
(4 
2 2 2 5 
w w,* - - 51 
P F R2/ l,m. ( ) 
It may be remarked here that as R > 0c0, w > w,, the 


natural plasma frequency. 


2.4 The case of free boundary 


In the case of a free boundary, the inside solution 
must be matched to the one outside with due regard 
to the boundary conditions. The equations governing 
the vacuum field can be obtained readily from Eqs. (40) 
and (41) by setting w,? =0. Thus H satisfies the 
equation 


aH /7¢7+1) 7 

as + | 232 — ( 2 Ja - 0, (52) 
where &, = w/c. The solution of Eq. (52) is given by 
H : = r'ls[ Hy HS), (er) + Hoy H, 2), (ker)]; (53) 


where H,and H,® are the two kinds of Hankel func- 
tions whose asymptotic forms are given by (4,r > 1) 


em) ~ (= 


aA 


and 
2 \'"/2 ilz—(a/4)(2r 4 
H,2)(z) ~ (2) “e 55 
(9) ~ (25 (55) 
In Eq. (53) Hy and Hy, are constants of integration. 
We have assumed the time dependence of the wave to 
be of the form exp — ‘wt. Therefore, in order that at 
large distances from the sphere we have only an outgoing 
wave, we must set 7 0. Thus the permissible 


solution is 


01 — 
H = Hyr'hH,? (ker). (56) 
Applying the boundary condition (47), we then obtain 


d 


d | - 
ap Wn, y, (Ait) = Gln HS, (he). (57) 


We now restrict ourselves to the special case where 
k,R > 1 and &,R > 1. We may then use the asymptotic 
expansion 


J. (%) ~ ( 


\i 


yA 2 . I ~~ 
a cos | z 4 (er 1) 


. 


(58) 


and Eq. (54) for H,"". Eq. (57) then leads to 


a 
k, tan |k, R— = (/+4+- 1)] =shy. (59) 
Using the relation 
t l ‘<< 60 
arc tan 2 -In , 
~ 2i0 1—% a 
we obtain from Eq. (59) 
R a 
: (w? — w,?)" 3 (2n+/-+-2) 
1, w+ (w? —w,?)'): 
+ 5;In \ > =. (61) 
4i wW— (a —w,*)'s 


The notion of a radiative solution implies that, at a 
fixed point in space, the amplitude of the perturbation 
should die out asymptotically in time on account of a 
finite amount of energy available for radiation in any 
admissible initial perturbation of the plasma. This 
implies that m has a negative imaginary part. We con- 
sider oscillations close to the plasma frequency and set 

@ = W, +6- -iV, (62) 


~ 


where y is positive, and 6 and y are both assumed small. 
On substituting Eq. (62) into Eq. (61) and retaining 
terms to the lowest significant order, we obtain 
fy) ¢ . R2 1| 
= @ _ 
y Zw,R\ > A 


(63) 


and 


f R Cc / R2 2 
2y12—+ of —1| 
¢* 2mw2R\ ? & 


/\12 
: [n(x +1+5)| . (64) 
We now consider the physically interesting situation 
where the time taken by the wave to traverse the radius 
of the sphere is much larger than the period of a plasma 
oscillation. That is, 


R T, . 
= Wy 7 1, (65) 
where 
1 
T, -and T, (66) 
¢ @ 


Dp 


From Eqs. (63) and (64) we obtain to the lowest signifi- 
cant order 


y=, (2) E \» +1+ 2) (67) 
and 
o ; M, (=) a |» +1 +4 Alt (68) 





This establishes the possibility of radiation by plasma 
oscillations. However, it is not possible to determine 
the intensity of the radiation emitted within the frame- 
work of the linearized theory presented here and in the 
absence of any knowledge of the excitation mechanism 
for the oscillations. 

In the preceding calculations it was tacitly assumed 
that w? + w,”. There are, however, situations for which 
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w? = w,”. It will be recalled from §2 of Section III 
that within any region of a cold plasma of uniform 
density one can construct arbitrary wave packets of lon- 
gitudinal waves characterized by 


V x E=0 and V-E+0, (69) 


which are non-dispersive and which do not couple to 
the radiation field characterized by 


V-E=0. (70) 


The longitudinal oscillations, of course, can exist also 
in the interior of the sphere considered here, but clearly 
they do not radiate. 

Those radiating modes involving the functions F 
and G, discussed previously, are associated, however, 
with a fluctuating surface density and do radiate. 

The sharp transition between plasma and vacuum 
assumed here is merely a convenient mathematical 
device for representing a rapid, but continuous, tran- 
sition in density. Clearly the radiation associated with 
the fluctuating surface charge can be viewed as arising 
from a coupling between longitudinal and transverse 
fields induced by the density gradient. In fact, one can 
verify readily that for the case where the density is a 
function of position, the decomposition into solenoidal 
and irrotational waves is no longer possible. For the 
general case, however, the calculations become much 
more difficult and the radiation associated with charge 
separation has, as yet, been explored only incompletely. 


3. OSCILLATIONS IN NON-UNIFORM PLASMAS 


We now consider the longitudinal electron oscillations 
in non-uniform plasmas in which temperature and/or 
density gradients are maintained by some agency. The 
model we take for the steady state is one in which the 
electron velocity distribution is assumend to be iso- 
tropic against a fixed background of positive ions, 
resulting in an electric field which balances the pressure 
forces in the equilibrium state. The force balance con- 
dition is, then, 


Vp=—NeE, (71) 
where p (= N@) and E denote the pressure and the 
electric field in equilibrium respectively. It is also 


assumed that there is no magnetic field. If we assume 9 
to be constant, Eq. (71) can be integrated to give 


» » e@ 
N = Noexp = (72) 


where @ is the electrostatic potential given by E = — Vp 
and JN, is a constant. Eq. (72) is the well known density 
law for an isothermal atmosphere. On the other hand, 
if we assume JN to be constant, Eq. (71) gives the energy 
integral 


E= O — eg. (73) 
In general, it follows from Eq. (71) that 
p=f(N), (74) 


where / is an arbitrary function of its argument. 
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The equations governing the small departures from 
the state of equilibrium are (cf. Eqs. (II-47), (11-48), 
and (II-119)) 

r) 
t) 
ov 


mN ry; V -p, —#eE — NeE, (76) 


n > 
7+ V-(Nv) 0, (75) 


and 
r) —— ‘a 
“7 p,+V-(vp)I+pVv+p(Vv)?=0, (77) 


where # denotes the perturbation in the electron density 
and p,, the perturbation in the stress tensor. The per- 
turbation in the electric field E, for the longitudinal 
oscillations under consideration is given by the Poisson 
equation 

V -E, 47en. (78) 


Introducing &, the Lagrangian displacement of an 
element of the fluid, Eqs. (75), (77), and (78) can be 
integrated to give 


n -V -(NE€), (79) 
p, = —V-(€p) 1—p(VE)—p(VE)", (80) 

and 
E, = 42e NE. (81) 


Using the foregoing results in Eq. (76) we obtain 
1 1 
V-( | + —V > ( -V 
p ¥ (PR) +, Va(E- VP) 
4 , P , 
— nl ValEVN)+Y. pVE+p(VE)"| —wNo,7E. 


(82) 


2 
mN : . pV 


We now multiply Eq. (82) by —* and integrate over 
the volume of the plasma. Assuming that & vanishes at 
the boundary, Eq. (82) leads to 


| mv Edt \{ FE V -(p&) | F 1§-Vp | 
+(x) eva] 


2 


‘EVN | +2p(VE)sym: VE* + wNw,?\E thar. 
(83) 


Now all the terms except the second on the right hand 
side of Eq. (83) are positive definite. Thus the system 
will be unstable only if the second term is the dominant 
one on the right hand side of Eq. (83). 


3.1 Density gradients 


Let us first specialize to the case where there exist 
density gradients in a plasma at uniform temperature. 
Then f = NO, where @ is a constant and the second 
and the third terms on the right hand side of Eq. (83) 
cancel each other. Consequently, m? is always positive 
definite and the system is stab/e. 


3.2 Temperature gradients 


We now consider the case when 0 = O (x) and A 
is constant. Then the third term on the right hand side 
of Eq. (83) vanishes. For instability to occur the second 
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term must be capable of having a large negative value 
compared to the fifth term. In order of magnitude, if / 
denotes the characteristic scale of variation of 9, then 
for instability we must have 
: NO 
uN? < Rp (84) 
2 


Eq. (84) can be written in the form 


o] 


<—,=A,2. 
MO, 


hp (85) 

Thus, for instability to occur we must have significant 
temperature gradients over lengths smaller than the 
Debye length. However, at such small lengths the 
macroscopic theory is inapplicable. 

We may thus conclude that electron oscillations in 
non-uniform plasmas with arbitrary density or tem- 
perature gradients are macroscopically stable. 


V. Non-Linear Oscillations 


We now consider large amplitude electron oscillations 
in a uniform plasma and, for the sake of simplicity, 
restrict ourselves to one-dimensional oscillations only. 
Furthermore, we neglect the thermal motions of the 
electrons. Let x, denote the equilibrium position of an 
electron and (x) denote its displacement (not neces- 
sarily small) from the equilibrium position. The displaced 
position of the electron is then given by 


== ty + A lie): (1) 


In moving the distance X (x,) the electrons in the x, 
plane pass over an amount of positive charge eX per 
unit area, where NV denotes the equilibrium density of 
the electrons. If the ordering of the electrons in the 
x-direction is not changed, then all electrons which 
were originally on the positive side of a given electron 
remain on the positive side, while all those which were 
on the negative side remain on the negative side. We 
shall assume that the system is of infinite extent and that 
no charge enters or leaves the system at infinity. The 
equation of motion of an electron is given by 


aX 


m 
df? 


ek, (2) 


where, by Gauss’s theorem, the electric field is given by 
E=42NeX. (3) 


From Eqs. (2) and (3) we obtain for the equation of 
motion 


aX 


TP. w,?X. (4) 


where w,,” — 42.Ne?/m. The general solution of Eq. (4) 
is given by 


X = Xj(xp) sin@,t + X,(x9) cos @,¢, (5) 


where X, (x,) and X, (x») are some arbitrary functions 
of x». Thus each electron executes simple harmonic 
motion about its equilibrium position, independent of its 
amplitude and the motion of the neighboring electrons 
so long as the ordering of the electrons is maintained. 
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The condition that the ordering is maintained can be 
readily derived in the following manner. Let x, + 


denote the displaced position of an electron originally 
at x,, while x, Ix, + X + AX, denotes the 


displaced position of an electron originally at x, + Ax». 
In order that the ordering be maintained, it is clear that 
Ix, + AX > 0; we thus obtain 
ox , 
~~ (0) 
0Xo 


as the required condition. 
An interesting special solution to Eq. (4) arises when 
X, = Oand X, 


A sin kxo, (7) 


where A is a constant and & denotes the wave number 
of the disturbance in the x-direction. The solution of 
Eq. (4) is then given by 


xX =X + Asinkxp cos w,/ (8) 
and 


E=4aNeA sinkxgcos w,/. (9) 
At ¢ = 0 Eggs. (8) and (9) reduce to 


x =X + Asinkxo (10) 
and 
E=42NeA sinkxp. (11) 


We may solve Eq. (10) for x, as a function of x (this 
is most simply done numerically) and thus obtain / as 
a function of x. The ratio //E,,,, is plotted in Fig. 16 
as a function of x for various values of A. For small 
values of A the curves are essentially sine waves. For 
large values of A one gets a multiple valued curve which 
is, of course, meaningless. What has happened is that 
the ordering of the electrons has been destroyed. The 
crest of the wave falls into the trough and we have a 
wave breaking process somewhat similar to the breaking 
of water-waves, and we may expect the waves to be dis- 
sipated by this breaking mechanism converting the wave 
energy into the thermal energy of the electrons. 
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Fig.16 The behavior of E/E 
values of A at time ¢ = 0. 
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A convenient mechanical model which exhibits the 
breaking process is a series of elastic pendulums con- 
strained to move along a straight line. For small dis- 
placements the pendulums swing independently of each 


or 
b 
at 


| 


;” 


US 


other. However, if the amplitude of the displacement 
becomes large, the bobs of the pendulums strike one 
another, destroying their periodic motion. 

The breaking of large amplitude waves has been 
investigated by Dawson for the one dimensional case 
under various initial conditions. Firstly, consider the 
case when the initial displacement vanishes and the 
initial velocity is given by 


| V9 sin & Xo. (12) 


From Egs. (6), (8), and (12) it follows that the critical 
breaking amplitude is 


VY .=—w_/k. (13) 


In his numerical calculations Dawson investigated the 
case when the amplitude of the wave was taken to be 
7°, larger than the critical amplitude given by Eq. (13). 

Secondly, we consider the case when, in addition to 
the initial velocity given by Eq. (12), there is super- 
imposed on the particles a certain amount of random 
motion. The random motion was taken to be 


Vin = Vn exp (—V2/V,2) (14) 


whre I|”,, is a constant which is proportional to the 
thermal energy. The effect of thermal motions on the 
breaking of waves is exhibited in Fig. 17 for the case 
where the thermal energy is 4°, of the wave energy 
The lower curve in Fig. 17 shows the wave form for the 
wave at a time two plasma oscillations later. It is clear 
from this curve that the amplitude of the wave has 
reduced considerably and that the wave energy is being 
converted into the thermal energy of the electrons. 
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Fig.17 The behavior of F as a function of x when the amplitude 
of the wave is 7°,, larger than the critical amplitude at which the 
wave breaks, for different initial conditions. The crosses show 
the wave form for a cold wave just as it is starting to break. The 
squares show the wave form for a wave with initial random energy 
equal to 4°,, of the initial wave energy, also at the start of breaking. 
The circles show the wave form for a warm wave two plasma 
periods later. 


Fig. 18 shows the velocities of 270 pendulums after 
three oscillations for an initially cold wave. It is apparent 
that most of the particle velocities lie on a sine wave. 
However, there exist a number of particles with very 
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Fig.18 The velocities of the particles after three oscillations 


high velocities, both positive and negative. These are 
particles which managed to get in just the right phase 
with the wave so as to ride it, and they account for the 
major loss in the energy of the wave. 

In order to get a quantitative idea of the energy loss 
of the wave, we Fourier-analyze the displacement and the 
velocity in terms of the equilibrium position, i.e., 


> 

F 

/ 4 
’ 


NUX,; 
j 1 sin 15 
has 3 
n=l 
and 
\ 
-— NX; 
1 (x;) B, sin; (16) 


I 


where L denotes the linear dimension of the system. The 
Fourier co-efficients A,, and B, are given by 


’ 








\ 
? + ° 
z ' nax a 
A, ' X (x;) sin ’ 7 
N+1 4, 8 
n 1 
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2 " nx 
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Fig. 19 The ratio of the energy in the first harmonic to its initial 
value for a number of different cases. 
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For the Fourier analysis it is convenient to use the 
equilibrium position rather than the actual physical 
position since for such a choice there is no coupling 
between the different modes without breaking. The 
results of the Fourier analysis are shown in Figs. 19 
and 20. Fig. 19 shows the ratio of the energy in the first 
harmonic to its initial value (for a number of different 
cases) as a function of time. The cold wave loses very 
little energy during the first oscillation and a half, after 
which there is an almost continuous decrease in energy. 
As is clear from the diagram, the presence of a small 
amount of random motion causes the wave energy to 
drop to about 10% of its initial value over a period 
of one oscillation. 

The absolute values of the amplitudes of the various 
Fourier modes for the case of forty-five pendulums 
without any random motion is shown in Fig. 20 after 
2.25 oscillations. As is clear from the figure, there is no 
tendency to feed energy from the fundamental into any 
one of the higher modes. The curve is flat (with only 
small fluctuations superimposed on it) all the way out; 
only the fundamental has appreciably more energy than 
the rest of the modes. The energy of the wave seems 
to go into individual particles rather than into Fourier 
modes. 

One may expect the system to exhibit a Poincaré 
cycle and return to its initial state after a sufficiently 
long interval of time. Simple probabilistic arguments 
indicate that this interval of time will indeed be very 
long even for a small number of charge sheets. For 
example, such arguments indicate that a system of five 
sheets should spend approximately 1/400 of the time 
in a state for which the energy of the fundamental is 
more than 95%, of the total energy. A system of nine 
sheets would spend only one part in 1.6 « 10° of its 
time in such a state. 
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Fig. 20 The absolute values of the amplitudes for various Fourier 


modes for the case of forty-five particles per half wave length after 


2.25 oscillations. 


The recurrence time was investigated by Dawson 
following systems of five sheets and nine sheets for 
about sixty oscillations. The systems were started in the 
fundamental mode with amplitudes larger than the 
breaking amplitude. It was found that on the average 
there were five crossings per sheet for the five-sheet 
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case, and ten crossings per sheet for the nine-sheet case, 
in one-half an oscillation time. The state of the system 
was sampled for 120 values of the time. It was found that 
the five-sheet system returned to a state in which the 
fundamental had 95%, of the energy once during this 
time, while the nine-sheet system returned to a state in 
which the fundamental had half the energy once during 
this time. Fig. 21 shows the sum of 120 time-averaged 
velocity distribution functions for the nine-particle 
system. The smooth curve in Fig. 21 shows what would 
be expected if the motion was roughly ergodic. These 
results seem compatible with the probabilistic estimates. 

The foregoing discussion indicates that the breaking 
of waves can be an effective mechanism for the dissipation 
of wave energy into the thermal energy of electrons and 
that the presence of a small amount of random motion 
can greatly enhance this mechanism. 
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Fig. 21 The velocity spectrum of the particles. 
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RAYLEIGH-TAYLOR INSTABILITIES OF A MAGNETICALLY 
ACCELERATED PLASMA 


T.S. GREEN AND G. B. F. Nisietrr 


A. W. R. E., ALDERMASTON, BERKSHIRE, ENGLAND. 


High speed photography has been used to study the compression of a deuterium plasma by an axial magnetic 
field generated by a wide single-turn coil. In this experiment the inductance of the coil is small compared with 
that of the external condenser bank so that the plasma implodes at a constant acceleration of 5 X 10! cms/sec?. 
Flutes which develop on the outer surface of the plasma are interpreted as magnetohydrodynamic Rayleigh-Taylor 
instabilities produced by the inward acceleration of the interface separating magnetic field and plasma. The observed 
growth rates are in agreement with simple theory. It is shown that the ratio of the diffusivities of magnetic field 
and momentum determines whether the plasma resistivity or its viscosity is the dominant mechanism for damping 
the short wavelength instabilities. At the temperature and densities of this experiment the resistivity is primarily 
responsible for the damping and the wavelength of “maximum instability” calculated on this basis agrees with the 
observed value. 


Introduction constant acceleration and flutes which subsequently 
The rapid compression of relatively dense plasmas by develop on the ome surface of the plasma are — 
means of axial magnetic fields has recently received con- pasted as Rayleigh-Taylor instabilities produced by oe 
siderable attention [1]—[9] principally because this is a inward acceleration of the interface between the magnetic 
possible method of heating a plasma to thermonuclear field and po plasma. These instabilities -_ analysed ; 
temperatures. This paper describes a high speed photo- the basis * Taylors theory and the damping by resi- 
graphic study of such a compression process in which stivity and viscosity is discussed. It is shown that the 
the plasma is pinched to the centre of the discharge tube observed growth — and wavelength are in accord 
by a magnetic field generated by a wide single-turn coil. with simple theory. 
A stationary plasma confined by a uniform axial magnetic , 
field is expected to be neutrally stable against hydro- Experimental Arrangement 
magnetic disturbances [10] but in the absence of crossed A diagram of the experimental system is shown in 
magnetic fields there is no restoring force to prevent the fig. 1. A copper coil 21 cms wide surrounds a quartz 
development of interchange instabilities of the flute tube 8.0 cms internal diameter and is connected by a 
type [11]. In the present experiments the plasma col- single spark gap switch to a bank of capacitors sup- 
lapses towards the centre of the tube with approximately plying 3.5 kilojoules at 24kV. The total circuit in- 


QUARTZ DISCHARGE TUBE 


CONTAINING GAS AT LOW PRESSURE 
HIGH SPEED CAMERA 


























SIGHTED IN DIRECTION an “| MAGNETIC PROBE 
OF ARROW | LEADING TO AMPLIFIER 
AND OSCILLOSCOPES 
t m4 —-- ~ 
I oad 8cms ee 
SINGLE SPARK GAP SWITCH 
anne 
Fig. 1 
TOTAL CIRCUIT INDUCTANCE 0-14 yH CONDENSER BANK 
COIL INDUCTANCE ~ 0-O02pH (24kVv, 12mF, 35kJ) Schematic diagram of the experimental arrangement. 
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Fig. 2 

Streak photograph of a discharge 
in deuterium at an initial pressure 
of 80 microns. 
field is 13 kilogauss. 


ductance is 0.14 “H of which the coil contributes about 
0.02 “#H; thus about one-seventh of the voltage appears 
initially across the coil and since the inductance of the 
external circuit is so much greater than the coil induct- 
ance, variations of the latter during the discharge have 
little effect on the voltage and current waveform. A 
magnetic probe at the centre of the coil measures the 
field on the axis of the tube. 

Discharges in deuterium at initial pressures in the 
range 50—200 microns have been photographed axially 
as shown in fig. 1 using two high speed rotating mirror 
cameras. The streak camera records motion across a 
diameter and has a writing speed of up. to 20 mms per 
microsecond. The other camera is fitted with 30 separate 
lenses and gives a series of photographs of the whole 
tube cross-section at a framing rate of 8 X 10° pictures 
per second. The images are stabilised by using a Miller 
optical system and are 8 mms in diameter. 


Experimental Results 

Fig. 2 shows a streak photograph of a discharge in 
deuterium at an initial pressure of 80 microns. The first 
half-cycle in which the gas is ionised is not shown on 
the photograph; the gas current starts to grow at the 
beginning of the second half-cycle of the external cur- 
rent and this time is labelled ¢ = 0. The outer edge of 
the imploding luminous sheath is interpreted as the 


Peak magnetic f r 
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MICROSECONDS 


boundary between magnetic field and plasma and has a 
constant inward acceleration of 5 X 10™ cms/sec? lasting 
for about a microsecond. The inner boundary of the 
sheath is taken to be a shock wave which reaches a 
terminal velocity of 5.5  10® cms/sec. After the shock 
has reached the centre of the tube the plasma column is 
further compressed by the still rising magnetic field and 
then oscillates for about three periods before it expands 
to the wall when the external field drops. As the field 
changes sign and increases again a new implosion begins 
and the compression process is repeated. The magnetic 
probe on the tube axis becomes luminous when the 
plasma reaches the centre at about ¢ = 1.5 “sec. One of 
the flute instabilities can be seen as a luminous streak 
separating from the lower half of the plasma at about 
this time and reaches the wall at ¢ = 2.5 msec. 

Fig. 3 shows framing camera pictures of a discharge 
which is similar to that of fig. 2 except that the initial 
pressure is greater and there is no magnetic probe pre- 
sent. The upper row of photographs (up to 2.40 “ secs) 
show successive stages of one discharge whilst the lower 
row consists of photographs selected from two other 
camera runs. The discharges were repeatable in their 
characteristics so that fig. 3 shows faithfully the progress 
of a typical discharge during the first half-cycle and the 
beginning of the second. 

The cylindrical implosion shows up clearly and the 
plasma is seen to be concentrated in a cylindrical annulus 
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Framing camera photographs of a discharge in deuterium at an initial pressure of 150 microns taken using a high speed 


camera which viewed the discharge axially. The time in microseconds from the beginning of the first implosion is shown under 
each photograph. The exposure time of each frame is 0.12 microseconds. 
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with an inner core of trapped axial magnetic field. The 
framing camera has a depth of field of about 10 cm and 
was focussed on the mid plane of the coil. For the first 
two microseconds the plasma is in sharp focus but be- 
comes increasingly blurred towards the end of the half- 
cycle; this is because the plasma is being ejected from 
the ends of the coil. During the new implosion at the 
beginning of the next half-cycle the plasma is again in 
focus. 

Irregularities in the plasma cylinder develop during 
the implosion phase and these subsequently become 
well-defined on the outer surface of the plasma and are 
believed to be instabilities of Rayleigh-Taylor type. The 
irregular inner surface of the quartz discharge tube prov- 
ides nonuniform starting conditions and these perturb- 
ations grow in response to the inertial forces produced by 
the inward acceleration of the interface between plasmaand 
magnetic field. The growth time of the instability is a 
fraction of a microsecond and the typical wavelength, 
or wavelength of “maximum instability,” is about 
1.5cms. The shape of the instability—sharp flutes of 
dense fluid separated by round-ended troughs of the 
lighter fluid—is characteristic of the later stages of 
Rayleigh-Taylor instabilities. 

The temperature of the plasma in fig.3 has been 
estimated in a number of different ways. The shock 
velocity gives a temperature of 3 X 104 °K but this is a 
lower limit since at these temperatures joule heating is 
comparable with shock heating. Moreover estimates of 
the temperature from shock velocity are imprecise in 
this range where the temperature is a sensitive function 
of the shock velocity. The temperature has also been 
estimated from the measured damping of the radial 
oscillations [12] assuming that this damping is purely 
resistive; this method gives an estimated temperature 
of 6 X 104°K. Pressure-balance calculations yield a tem- 
perature of 7 X 104°K and a particle density of 1017 cm’, 
Thus none of these estimates is precise but they indicate 
that the temperature is of order6 104° K and the particle 
density of order 1027 cm”. 


Growth Rate 


The instability of an accelerated fluid interface was 
first studied by Rayleigh [13] and the theory subsequently 
developed by Taylor [14]. Experimentally, observations 
on liquid-liquid interfaces have been made by Lewis[15] 
and Allred and Blount [16] among others. If @, and @, 
are the densities of the lighter and heavier fluids respec- 
tively and a is an acceleration applied normally to the 
interface, then an initial disturbance of amplitude A, 
and wave number & (A = 22/A where / is the wave- 
length) is amplified in the form 


A= Ayer (1) 
where 
C2 — O1 ,|2 
n= + 0.) 2 | (2) 


The interface is unstable only when the acceleration is 
directed from the lighter to the heavier fluid and equa- 
tion 2 shows that the growth rate increases indefinitely 
with increasing wave number (decreasing wavelength). 
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Kruskal and Schwarzschild [17] were the first to dis- 
cuss Rayleigh-Taylor instabilities of an accelerated inter- 
face between an electrically conducting fluid and a 
magnetic field and they described in detail the way in 
which the inertial forces produce charge separation and 
particle drift which drive the instability. They showed 
that though the magnetic field produces the plasma 
motion the gross nature of the instability and its growth 
rate are the same as for ordinary fluids. 

We now apply Taylor’s theory to the observed in- 
stabilities though this comparison can only be ap- 
proximate since Taylor’s theory is a linear theory for an 
incompressible fluid of infinite extent and none of these 
idealizations apply in our case. Since we can regard the 
magnetic field as a fluid of zero density equation 2 be- 
comes 
22a \: 


a 


In our experiment is 1.5cm and ais5 < 10!" cms/sec?. 
Thus the characteristic growth time is 


n= (ak)” (3) 


1 A \3 3 
i=——} = 2 X 10-* secs (4) 

n 22a 

which is of the order of the observed value. 

Conditions in this experiment are particularly favour- 
able for the growth of Rayleigh-Taylor instabilities be- 
cause of the constant acceleration for a period of one 
microsecond, and since the instability grows like e”' and 
n is proportional to a’:, the shorter the time taken to 
reach a final velocity i.e. the greater the impulsive nature 
of the acceleration, the less the growth of the disturbance. 
In the limit of a pure step-function in velocity there is 
no growth at all. 


Damping Mechanism 


Taylor’s theory ignores effects such as surface tension 
and viscosity which will eliminate or damp out the 
shorter wavelengths. However, subsequent investi- 
gators [18]—[22] have considered these effects and shown 
that surface tension and viscosity act as damping proc- 
esses which produce a wave number &,, of “maximum 
instability.” The effect of viscosity, for example, in- 
creases with decreasing wavelength because the shearing 
forces increase, and the wavelength of maximum in- 
stability is essentially that wavelength at which the rate 
of doing work by inertial forces is equal to the rate of 
energy dissipation by viscous forces. 

In a plasma the instabilities will be damped by the resis- 
tive diffusion of magnetic field into the plasma as well 
as by viscosity, and the inertial, viscous and resistive 
contributions to the growth rate # can be written: 

Inertial term: (a k)’!s 

Viscous term: vk? 

Resistive term: k/4 206 
where y is the kinematic viscosity and a is the electrical 
conductivity of the plasma. The inertial term represents 
the rate at which the instabilities grow in response to 
inertial forces whilst the other terms are the reciprocals 
of the characteristic times for momentum and magnetic 
field to diffuse across a wavelength and are obtained by 
observing the similarity of the diffusion equations for 
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momentum and magnetic field, a similarity pointed out 
by Batchelor [23]. The relative importance of viscosity 
and resistivity as damping processes is given by the 
dimensionless parameter P where 


P= 4x00 = K_/R, (5) 
P is the ratio of the diffusivities of momentum and 
magnetic field and is also the ratio of hydromagnetic to 
hydrodynamic Reynolds number. 

The conductivity and viscosity of a fully ionised plasma 
are given by Spitzer [24] and Marshall [25] and we may 
write for deuterium, 


1-16-2°/: (AT,) ( 
oO " er (6 
sm, cz ln A 
Le 5(k T;) : ~ 
) i 1 4 (/) 

nm, 87')nm;leAlnA 


where 7, and 7; are the electron and ion temperatures, 
m, and m, the electron and ion masses, # is the ion density, 
e is the electronic charge, ¢ the velocity of light, & is 
Boltzmann’s constant, /« the dynamic viscosity and In A 
is the logarithm of a collision integral tabulated by 
Spitzer which is approximately seven for our conditions. 
These equations show that the electrical conductivity 
depends on the electron temperature and the viscosity 
on the ion temperature though at our temperatures and 
densities relaxation times are so short (much less than a 
microsecond) that 7, and 7, are effectively the same. 
Moreover collision mean free paths are short compared 
with the Larmor radii so that it is permissible to use a 
simple scalar conductivity and viscosity. 
Inserting numerical values for o and » gives 


3.2 x 10-5. 
P=420Y [4 (8) 
n 
showing that P is very sensitive to temperature. The 
damping effects of viscosity and resistivity will be ap- 
proximately equal when P 1 and for deuterium at a 
particle density of 107 cm’? this corresponds to a tem- 
perature of 2 x 10°°K. At temperatures above 2 x 10°°K 
the instabilities are damped primarily by viscosity and 
at lower temperatures by resistivity. Thus at our esti- 
mated temperature of 6 107 °K it is resistivity that 
determines the wavelength of greatest instability. 
Equating the inertial and resistive terms in the growth 
rate gives 
1 


A (2/2)* a *o * (9) 


m 


showing that /,,, is insensitive to a. Since a is proportional 
to 7,"2, 4, is inversely proportional to the electron 
temperature. Putting /,, = 1.5cms, a=5 
and InA = 7 gives 7, 5.8 x 104 
with the estimated value. 

\t higher temperatures where viscosity provides the 
damping we should expect 


10!? cms/sec” 
K, in agreement 


A,, =22a sys (10) 


so that the wavelength is again insensitive to the ac- 
celeration but since y» varies as 7;*/:, A,, should vary as 


/, sand should rise rapidly with increasing temperature. 
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Conclusions 

Analysis of high speed photographs of the instabilities 
which develop during the compression of a plasma by a 
rapidly rising uniform axial magnetic field suggests that 
these are produced by inertial forces resulting from the 
inward acceleration of the boundary between magnetic 
field and plasma. The estimated temperature and density 
of the deuterium plasma is6 & 104°K and 10!7 cm™3 and 
under these conditions the instabilities are damped by 
electrical resistivity. The agreement of the observed 
growth rate and the observed wavelength with simple 
theory tends to support this conclusion. In more power- 
ful discharges the instabilities should be damped by 
viscosity and this increases rapidly with temperature. 
The growth of the instabilities would also be limited by 
a more impulsive acceleration. 
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UN CRITERE NECESSAIRE DE STABILITE HYDROMAGNETIQUE 
POUR UN PLASMA EN SYMETRIE DE REVOLUTION 
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FRANCE 
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Pour ce travail, on établit une condition nécessaire de stabilité pour un plasma en symétrie de révolution. Ell 


correspond a des déplacements localisés au voisinage des surfaces a pression constante sur lesquelles les lignes d 


5 ae 


force magnétique sont fermées. Ce critére se réduit a celui de Suydam en symétrie cylindrique. 


Introduction 
En géomeétrie cylindrique B. R. Suydam [1] a donné un 
critére de stabilité qui est en fait un critére vis-a-vis des 
déplacements limités au voisinage d’une surface 0 
tante, surface qui dans ce cas est un cylindre. 
En effet nous exprimons le déplacement & sous la 
forme 


cons- 


E cit), Cott. $7), et 


On peut dans le principe variationnel donné par Bern- 
stein et autres [2] intégrer sur # et 


x et minimiser algé- 
briquement sur &, et ¢, et obtenir 


3 (rfé’ +2) 
SW onde S 1 (f2__ p) = 
Ou ben \ 2rdr l ye 1. k22 ( f )é 
| / f= KB,+ (mr) B, 
ou AB, (m|r) Bo 
|) b= (2/r) J, Be 


k et m étant convenablement choisis, f(r) s’annulera au 
voisinage d’un cylindre quelconque dans le plasma; soit 
r = ace cylindre. Supposons un déplacement localisé a 
son voisinage. Nous pouvons alors remplacer r, g(r) et 
h(r) par a, g(a), h(a) et f(r) par (r —a) f 


en posant 7 r a 


(a) et écrire 


i 


eae 
dW = 2a \ |‘ | A 


m +- k2 a 


g(a) &) ] 


a) & P - 
i) (a) ¢* 


, 
at 





perturbation étant limitée entre les 2 cylindres de 
rayon 7 Hs? 4) OUD est tres petit. 

0) W peut encore s’écrire en posant / 
grant par partie le terme en & &’: 


x et en inté- 


J 
? 73 f’ 


, 2a f2(a)y - 
Ol = \c x2 £2 4. 482) dx 
m= +- k a 
1 
ivec 
] , 2p’ 
1 > (2? —azf h (m* +- k2 a*)) 


aB? (wip) 
u Bs (rB,) 


j 


Pour obtenir la derniére forme de 
de la relation d’équilibre du plasma: 


p’ + B,B’, + BoB’s 


A ona tenu compte 


- (Be /r) 0 


Or on peut démontrer (voir annexe) que 
I 1 
1 &°) dx (A + 1/,) \ &2 dx 
a). 
1 1 
et qu’il existe des déplacements qui permettent d’appro- 
cher d’aussi prés que l’on veut l’égalité. 
\insi pour avoir OW > 0, il faut 


c’est l’expression de Suydam. Elle fournit ainsi un critére 
nécessaire et suffisant pour les déplacements localisés 
dans le cas de la géomeétrie cylindrique. 

En géomeétrie de révclution la svmétrie du probléme 
étant moins forte, une seule intégration dans |’expres- 
sion de dW peut étre faite facilement. Pour continuer 
nous ferons alors ’hypothése des déplacements localisés 
aux voisinages de surface a égale pression de type parti- 
culier correspondant au choix de Suydam en symétrie 
cylindrique. 

Nous obtiendrons finalement un critére qui généralise 
le critére de Suydam, évidemment nécessaire car cor- 
respondant a des déplacements possibles, mais probable- 
ment pas suffisant en général pour la stabilité vis-a-vis de 
tous les déplacements suffisamment localisés autour 
d’une surface 4 pression constante. 


Les équations fondamentales 
L’équilibre d’un plasma a l’approximation hydro- 
magnétique est déterminé par les équations suivantes: 

divB=0 
B x rotB grad p 
En symétrie de révolution, il est commode d’écrire: 


. a: <;) Ce 
B 1 grad y (1) 
r r 


qui est la solution générale dans ce cas de divB = 0 


r distance d’un point a l’axe de révolution 
z coordonnée le long de lV’axe de révolution 
ey vecteur unitaire perpendiculaire au plan méri- 


dien passant par ce point 
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p pression du plasma (définie a 4 7 prés) 
Tet y = sont des fonctions quelconques. 

La deuxiéme équation d’équilibre fournit les renseigne- 
ments suivants: 7 et p ne sont fonctions que de y, elle- 
méme définie par |’équation suivante: 


dp 1dT? 
, 2 ? 
“7 , dp 2dy (¢) 
avec 
2 2 
H 0 10, 9 


or? rdr | dx? 

Dans ce travail, nous étudions la stabilité des équilibres 
régis par ces équations. 

Le point de départ sera le principe d’énergie développé 
par I. B. Bernstein et autres [2] que nous prendrons sous 
la forme: 


. 1(° : 2 
OW, =; \ dr{LO+ (m-&) (j x n)P 
+ yp (div &)? — 2(n-&)?(j x n)-(B-V)n} (3) 
ou O = rot (§ < B) 


n = est le vecteur unité normal a la surface p constant 
E est un déplacement arbitraire, satisfaisant 4 la con- 
dition n- & = 0 4 la limite du plasma. 
La stabilité pour les déplacements qui ne déforment 
pas la surface du plasma sera obtenue si dW, reste tou- 
jours positif. 


Le critére de stabilité 


Nous emploierons le systeme de coordonnées ortho- 
gonal suivant, introduit dans [2], adapté a léquilibre 
dont on recherche la stabilité: Les surfaces de coordon- 
nées seront: 

y = constantes,  — constantes, 7 = constantes 

Les surfaces 7 = constantes sont choisies normales 
aux surfaces y et @ de telle sorte que (y, 0, 7) forme un 
systéme orthogonal direct. 

L’élément de volume s’écrit: 


dt = Jdpyd0dy 
Dans ces coordonnées, l’opérateur gradient s’écrit: 
n) 1 0 1 my) 
. N af. 
V Bes, rs murY) | JB “x dy 
B: longueur de la projection de B sur la plan méridien. 
De l’équation d’équilibre et de l’équation (1) on tire: 


J = rotB=jes+ BT’e, 


avec 
r Oo 
; “ 2 
J Joy l™ (4) 
et 
a fe i 
p= r r2 (5) 


p’ et 7” signifie a! et “7 
dy dy 


Un déplacement & est déterminé par ses composantes 


(Fy So ¢,) 
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En posant 


rBi,=— X,&/r=Y,&,/B =Zet Y—(T/r*)Z=—U 
et Popérateur 
iste i a i 
L’équation (3) s’écrit: 
OW, ; | dy dOdy | = plDXP 
aa, J 


rioUu . 2d /JX\}? aU 
r | [B L / 
Sex . oy 7) id E op ' rR, 
ou 
06 


yp rs) E i)? mee 
oa ‘. A. 3D 2 
ae ie (JX) +J ' Dz| 2/KX (6) 


avec 


K IT’ d(logr) / a(log B/) 


r2 oy r oy 
IT’ dflog(BJ/r)] ,a(log BJ) 
r2 oy oy (7) 


Remarquons que K < 0 est une condition suffisante 
de stabilité. 


Cette condition, dans le cas de symétrie cylindrique 


(B= Bir), J/=1/B, p= y(r)) vécrit TT’ <0 


condition bien connue dans ce cas. 


Posons: 
X = XX," cos mO + X4”™ sin m6 
m 
U = (U,"/m) sin mO 
m 
Z=XZ," cos mO + Z," sin mb 


m 


(U,™/m) cos m@ (8) 


On obtient pour chaque # (mw +0) 
- x ( aA,. jf, 
(OW x) im 7 \4vax i| ce. + ms x,) 


Xe 
+(° * ot x,) | 


dy 1) | 72 Re 
1100, . 0 J X,)* 
ae E 0x oy | 
7i1dU, , dd /]%X,)? 
Lae FE yi | Oy P 





aX, fs ~]? 


+ BU, ay + Rw | 


+ JB [Us + 8 + Fax] 
+ : so (JX,)+J/U,+ i + mi 2a] 
2) (XP + X¥){ = (log r) oe £ ‘= (log B/ )\t 


L’indice m a été supprimé sur X,X, U,U, 7,2, pour 
simplifier. 


Posons Kg = Ay" +i XS 
U T p m 4. il ” 
, 4 - 2,” 4+ tZ” 





de) Se 
tw we 


your 
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L’indice m de X UZ sera supprimé chaque fois qu’il 
n’y aura pas de risque de confusion. On peut écrire 


. x j 1 aA . st 
(OW) m > \ d a B dy ims X 
ri0U 8d ([JX\* oie ” | i . 
"Jom ox vn dy r )| +JR\U+ op | r BB 
yp| 2 JT Ze 


. vo . <Siae 
(J/X)+ Jl +5, ims 


+ J dy | 
) TT’ jd ny 
2 or - a , 
2I\au (log r) r2 r oy (le ¢ BJ) | P (9) 


Dans la suite nous utiliserons le résultat suivant relatif 
aux solutions périodiques d’une équation différentielle 
du type suivant: 


2... . 

+ia(zy)Z7+ b(y)=0 
07x ? 
On suppose que a(z) et b (7) sont des fonctions de 
période 7p. Si 

uv ho 
$a(y) dy \a(x) dx +2Ka(K entier) 

zx 


on a une solution périodique, de période 7,, et une seule 


1 Xr ho 4 
Z() x) exp/ x (10 
(7%) 1 expi §a(x) dx \ d(x) exp: | a(u) dudx (10) 
x % 
En effet, on peut vérifier directement que c'est ane 
solution. F:lle est périodique. En effet: 


» 
2) x 
ho a 


1 


i 
29+ ie §a(x)dx ) &(=expi | a(x) duds 
J . 


2 
Xo Lhe 


et faisons le changement de variable 


xXx=A +t 
uv Xot+v 
Z(x + Ze) 
1 x 1 xo : 
1 exp i a(x) dx \ b(t) exp? \ a(v)dv dt = Z(y) 
x x 


Enfin il n’y a en pas d’autre car on obtient la solution 
générale en ajoutant 4 cette fonction périodique parti- 
culiére une solution générale de l’équation homogéne 
(0x7/d7) +7a(z) x= 0, Cest-a-dire 


Z, : 
z Cant i \ 0 (x) dx\ 
; \ ; J 
vA) 
solutions qui ne sont jamais périodiques. 
Revenons 4 (6 Wr) ,, et minimisons d’abord sur Z. 
En posant 


f . sp JX) + JU+ 3 im! z) 
’équation d’Euler s’écrit: 
| Ps # 
dy 7m 72 Sj 
x 


Soit f= f(y) exp 





im \ (/T/r?) az| 


0 


f doit étre périodique en y (période y,). Ceci ne peut 
se produire que si 


§ (JT r)\dy=2kn 
§ indique l’intégrale prise sur une période. 


Sauf pour un ensemble de points particuliers (de mesure 
nulle) cette relation n’est pas vérifiée. Done / (y) = 0. 
La minimisation impose donc que f = 0. 

Encore faut-il trouver une solution Z convenable c’est- 
a-dire périodique en y déterminée par / Q. Or f/ 0) 
est de la forme étudiée ci-dessus 


oZ . JT, ee : 
¥* im 2Zt+/l bay) 0 
avec 
a(x) m fT |r 
. 0 
b(z)=Jl tay! 


fonction périodique de période 7, et 


$a (x) dx m § (/T/r?) + 242 en général 
La solution cherchée qui donne le minimum par rap- 


port a la variable 7 est donc 
z 1 

° 1 exp im § | JT/r*\dy 
i+ Xo x 
VPs em ee 
\/e + oy) *) exp im\ (J1 r2\du 


x 


a 


dx (11) 


e 








® 


INTERPRETATION DE § (JT /r*) dy. 

La ligne magnétique tracée sur la surface y = cons- 
tante et qui passe en un point M de cette surface a 
pour équation: 

‘i rdO fr 0 
Jdy B, Tr! 
Soit: 
do (JT/r?) dz 
A § (JT |r?) dz (12) 
Ainsi la quantité § (JT |r?) dy nest autre chose que 
langle @ pris entre deux point M et M’, M’ étant le 
point ot la ligne magnétique recoupe pour la premiére 
fois le paralléle passant par M. Cet angle est indépendant 
du point M. 

Les cas particuliers ou 9, = 2 4a signifient que la 
ligne magnétique se referme. En général 0, +242 et 
la ligne magnétique ne se referme pas. 

Posons alors 


ox r 
1} dy - 
0X iajt x (13) 
ox r2 
ay OO a i — 
= maxX oy \;3 x r x 
Entre les » ¢ & il existe la relation: 
07) de ae 
= < ) 
dy dy iméE—=( (14) 
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Moyennant cette relation, on peut tirer univoquement 
NX et la partir 4 e § périodiques en 7 en utilisant une 
fois encore le théoréme (10) appliqué a l’équation: 

oX imjJT ., 
r Xx é 
07x i 
On peut alors écrire: 


in a (° ( jel* ‘sae Pe 
OW pp, > \4vdz ape + let 


2 J 4 >) 7 \9 an ~ 
+B | + JX] 2)\X 2K) (15) 
avec: 


1 exp im $ (/T/r*) dz 
Z 


ho 





\ €exp im \ (JT /r?)\ du | dx 


vA vA 


Remarquons alors que si on pose: 


Z 
9 (7, a) \ (J T/r?\du 
on a: 
(4+ xo. 4) = 9(z,a) + % 
donc: 
0(4,a) = p(z,4) I z/Z0 
avec: 
PLL > for4) = P\Z>4) 


ce qui permet d’écrire: 


exp [¢”p(7,4@) imQe7 Zol 


- 
is exp(— /mp) 


\ €(x) ¢ imp(r,a) e—? MO," tz dx (16) 


Comme € (x) ¢ ’ est une fonction périodique en x 
(période 7,) on peut développer cette expression comme 
sult: 

Ls Bug, CAE Ke (17) 
k 


€(x) e imp ra 


et en reportant, on trouve: 


- % Akm 9 
i a : . pimp ya) pinky) yx g 
X te» ixk m0.) é 0 (18) 


k 


Remarquons que le changement du point de référence a 
se traduirait par un simple changement des a Pour 
la suite, nous écrirons p (7) . 

En analogie avec l’hypothése de Suydam [1], nous 
étudierons les déplacements (4, ) tels que 227k — m0, 
est voisin de 0. ; 


im* 


Soit une surface y = constante telle que 4, /2M/ ne 
soit pas rationnel. On peut toujours choisir les entiers 
& et m de telle sorte que 217k — mf, soit aussi voisin de 0 
que l’on veut. 

On peut supposer que ces déplacements sont parmi 
les plus dangereux. Ils coincident d’ailleurs avec les 
déplacements de Suydam si l’on se place dans le cas 
particulier de la symétrie cylindrique. 
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Aussi allons-nous étudier les déplacements définis par 


é Ape CLM? RAK Ke 
(19 
g vA) 
\ h ‘A imy Zinky zx, 
“> Pkm 7 (2 mQ,) ‘ x’ € z 
On constate que: 
4 ‘ho é “ 
ts 22k — mb), V 
avec: 
22k — m0, 
y 
Zo 
Posons alors: 
X oe t ivX ivoe™ 
Donc: 
0 = Oly) 
, oe A - - 
q mp (7) 2a k : C (yp) mi 4,4) vy+C(p 
(0 
C () étant une constante arbitraire. 
Introduisons la fonction / telle que: 1 AX 
01) ee imJT 
dy v7 AA + VDILV. a 4 } 
10h im]l 
a + ivh + J h| 
07 an 
OF .[. 0” . 0 og] 
X ji jy y 
oy l oy 0 OY} 
D’aprés la relation (14): 
_— O¢ 07) .[. ov J og oh 
ims X |/ L 7p y 
oy dor oy 0 oy ox 
im][l 
ivh J | 
JT’ y y ]T’ 10% vo’ ivdg 
: r2 * ; re mow Ym 0 m oy 
i oh vy. Pe 
- h 
MOY W r= 
J , , ] 
1 > Be X | Be) 
D’out: 
af [fT riyl’ 1dr. ivdg 
SW lwdyo* spe . t 
Ow Pm 2 \ PALE Ur 32 jir mow m OW 
7dh oY : yo’)? [ 2 
prs OP poly 4 A] _27x1 
MOY mM r2 pn | r Be j 
(20) 


Les calculs qui ont permis de passer de l’expression (9) 
de OW’,,, 4 Pexpression (20) sont valables en dehors des 
surfaces y = const. sur lesquelles les lignes de force 
sont fermées. Pour un ~ donné ces surfaces correspon- 
dent a des valeurs discrétes de y. L’équivalence entre (9) 
et (20) est donc établie si l'on exclut de l’intégration 
sur y un petit intervalle autour de chaque valeur singu- 
liére de y. Mais comme les quantités 4 intégrer dans (9) 
et (20) sont continues en y, les expressions (9) et (20) 
sont aussi égales sans restriction sur l’intégration en y. 
L’elimination du terme provenant de y p (div &)? neces- 
site une étude plus détaillée an voisinage des surfaces 
singuliéres de y. 


ow 











UN CRITERE DE STABILITE HYDROMAGNETIQUI 


Nous allons chercher un critére local, c’est-a-dire que 
nous nous limiterons aux déplacements localisés au 
voisinage des surfaces y 
de 0. 


Wo et tels que vy soit voisin 


ov 
\y vo) (5) 
7 


(y — Yo) ¥’o (21) 
0 

Ce choix de déplacement est analogue au choix fait 
par Suydam. En effet, le critére de Suydam est en fait 
un critére local au voisinage des surfaces r —r, ot v 
cylindrique = & + mB,/rB, (voir introduction). 

Ce y cylindrique correspond exactement au y défini 
ci-dessus en symeétrie de révolution. Les cvlindres r = r, 
définis par Suydam se généralisent selon les surfaces 
y y, 0U les lignes de champ magnétique sont fermées. 

Ainsi dans 6’; nous remplacerons y par (y Wo) Vo 
et les quantités telles que r, B etc. ... relatives 4 la solution 
d’équilibre seront prises constantes et égales a la valeur 
qu’elles prennent pour y 
en général de 7. 


yo: Elles dépendront encore 


= 2 - /1 Vo. Vo’ 0° 
™ ' 2 \ / ’ . Ze \/ | rot my! \Y Yo) WO 
Zdh JT 2 [ j 12 
P d + 32 Ty 4 ? 4 7 
MOY id JB? YI Ba 2/ Ky (22) 








Nous limiterons encore les déplacements en supposant 
la fonction / réelle. Dans ce cas le seul terme contenant 


ae , 1 /ds\2 
04/07 dans lintégrant est de la forme —, | ) et il est 
’ . me 07 


alors évident que l’expression sera minimum pour 
Wl 0; OW... ne contiendra que la fonction / ce qui 
permet de minimiser par rapport 4 elle sans difficulté. 


On trouve: 


” ae B2o? [r2 (JT Vo 
OW, > \“ pdy B72] J | 2 = 
e 2 
vga? 2T/JT . 
\4 Yo) mn oO cs r B on 
Vo’ 0’ ]j?T? cil 
\P — Yo) wy 6 }4 r BA | aja +25) 


En intégrant par partie le terme 9 0 par rapport ay 
on peut écrire OW’;.. sous la forme OW, C*R 


, 


R \dyl(y Yo)? 0 24 10} (24) 


avec 


Ld 


_ alter BB. 
C* 2\ lp ] T2 9% 





bo 
C\| nq 
~ 


Or, on a démontré (voir annexe) que pour o (1) 
o (— 1) = 0et oo’ continues 
I 
— “i l o7 
\ (20? 1o*)dy>|A t q) \ eax 


et qu’il existait des déplacements qui rendent la valeur 
du premiére membre aussi voisine que l’on veut de alle 
du second. 








On peut appliquer ce résultat al’expression OW’... et on 
obtient le critére suivant: 
»? 9 , jl ] Yo . 
® | B2 2} | 7 ~i-.-4 
O ) r Bil ) I TR a 
> | eee) fa 
En tenant compte de la formule (4), QO... s’écrit encore: 
> »” 
r2 Be si . 
/*B2 8d 7 v, \2 


0 | 27K 


4 dy| 7? 
<< $47! 7| po a 2 dy\ J 2 


Dans cette expression 


Vo l 0 E 7 = 
W dy oy r2 
en posant il / : 
r2 
: Pay J 7 § ud VA 
u“ > i“ 
0 y j d Z 


et en tenant compte des équations d’équilibre du 


systéme, on peut réécrire 0... sous une autre forme: 


> 0/ ° Jdy 
4 »* Ay ? 4/2 
v / f vy 4 / | Bb |1 (1 2/r? B*)| 
1 . Jdy 


4 ru Ys at 


(/2/r2 B?)| 
— Jdy 
T2p'u {) a .1 spe 26 

red fF BPufl (72/r? B2)| 

Dans ces expressions, il reste encore une fonction 
arbitraire, c’est le paramétrage des courbes y = cons 
tante, c’est-a-dire que l’on peut choisir arbitrairement 
pour 7 une fonction quelconque de 7». 

Z i \ Zo) 

Yo Ctant tel que: 


Jdyz lod xo aid 


Il nous reste 4 minimiser sur cette fonction / soumis 
a la restriction: 


A $dyz § (Fo) 4 Zo 28 
Dans l’expression (26) seuls les deux derniers termes 


varient dans la transformation 7 f (7). La minimi 


sation est algébrique et la valeur minimum pour ces deux 
termes s’écrit: 

1 
49 [1 + (72/r? B2)] (J ir?) dy 


[ofan ae r§ i 


dps P 2 B2 
anna s Jdyz 
pT? D 52 Bey (72/r2 B2)] 
L’expression finale minimum de Q ... s’écrit: 
1 S272. e Jd7/? 
C ; ; da; a) 
Y 4 . 1 72 Taz ay9 re A ie l . r2 2 
0’ f OJ a. an ( Ls a 
j 7 a y az j 7 32 & 


ds! B, s designant l’abscisse curviligne 


Notons que / dz 


des courbes 4 pression constante. 
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Le critére obtenu s’exprimera donc par: O ., > 0. 
Si nous nous placons en symétrie cylindrique on a / 
1/B,, T =rBy, dx devient dz, 0/op devient d/dr 
et #¢ = /T/r? = B,|(rB,). Toutes les fonctions ne dépen- 
dent que de r. On obtient dans ce cas le critére de Suydam 
en utilisant ’équation d’équilibre 
2 Dn'y2 
Lt <p a 
ob > 0 30 
4 7 rB2 rs ( ) 
ip du 


a 
‘et mw’ signifiant alors —~ et 
rs § dr dr 


Norte 
Ce travail terminé, nous avons recu un article de 
B. B. Kadomtsev [3] dans lequel l’auteur généralise 
aussi le critére de Suydam pour une géométrie torique 
queléonque mais en imposant 4 l’avance une forme 
particuliére de déplacement localisé. Le résultat obtenu 
revient dans notre méthode 4 supposer / — constante. 
Dans ce cas nous retrouvons en géométrie de révolution 
son résultat. L’expression (29) est différente et plus 
sévére car nous avons supposé seulement que / est 
réelle et minimisé dans ce cas. 
Le critére de B.B.Kadomtsev s’écrit avec nos 
notations: 


1 1 

4 *Jdyz $(/T2/r2) dz 

f r2 (1+ § JB dy 

_ *Jdye z , 
6 er 2 PJ dy 
Ope y= . § / B dy 
,po] ,. (9 J dz)? 

+S sy4%—P §/B8dy~ 

Annexe 


Il s’agit de démontrer que la quantité 
+] n 
; fx 72 dx 
5 (9) 1 
jp2dx 

est toujours supérieure 4 1/4 pour toutes fonctions » 
s’annulant 4 + 1 continues ainsi que leurs dérivées, la 
valeur 1/4 pouvant étre approchée d’aussi prés que I’on 
veut par une fonction de cette classe (C). 

Nous indiquerons les grandes lignes d’une démonstra- 
tion basée sur le développement de _y en polynomes de 
Legendre, orthogonaux sur (— 1, + 1). Prenons pour 
une combinaison linéaire des N premiers polynomes 


Phe) 


N 
F p3 ay P,,(x) (A 1) 
1 
Nv 
y(—1) = > (—1)" 4, = 0 
1 


D’ou 
N N 
>4n 0 S i. 0 
n i,m 8. n 0, 2,4 


En utilisant les propriétés intégrales des polynomes 
de Legendre on peut écrire: 





N (n m fe ee 
5 : Pe = (# 3) a,2 +2 2 a ayn (n+ 1) 
2 2 a,2/(2 + 1) 
1. Ry 
2 


N 
{Ry n’ayant que des # pairs 
R, = R,,, + R,, sires a, 
” | y; n’ayant que des » impairs. 


La suite traite R,,,, le calcul pour R,, est identique. 


En tenant compte de > a, = 0 on peut poser: 
ig = Ne 


ag = Xy— Xe 


ra 


ay Xn +2 Xn 
ay Xy 


Ry écrit alors: 
\ 


N—2 
. 2n+1 :. , 
Ry Xn | Xn oa a Xn+2) TT on Xy* 
. = n 
aVEC A, = 3et 
(2 — 3)? 
XK, = 2n—1 52 
en : = Xn 2 
ce qui fournit: 
2n+1., 1 
Xn ? “5 1 
= ? | 1 
2(1+ $+... x3) 


X, est donc toujours positif ce qui prouve que Ry lest 
toujours aussi. 

Ainsi: S> } pour tout yde la forme (A1). S’il existait une 
fonction de la classe (C) telle que S (9,) < } on pourrait 
trouver un de la forme (A 1) tel que [5 () — S (,)] soit 
arbitrairement petit, donc tel que S (y) < } ce qui est 
impossible. 

Il en résulte que S (7) > 4} pour tous les y de la 
classe (C). 

Pour montrer que l’on peut approcher d’aussi prés 
que l’on veut la borne 4, choisissons pour les x,, la 
forme suivante: 


1 5) 
Xn ( ] )n/2 n? TP Ln 
avec P > 0 et 


t 2 2n+1 
= - - + avec /, > 0. 
n 


a 
Alors 
N--2 N—2 
1 1 s t,2 
2 a D2 n 
Ry > Kn by | n? (n + AF] <8 »2 mes 1 


carn, <2n 


D’autre part: 


? 
2 e < 
D,, ay 2n 1 


D 


D 


O 
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N\ 2 
x2 Ni 2 42 Ainsi pour P > 0, & (NV) tend vers zéro lorsque N > 4 
~~ * £2), * 
n 2 R, ‘ 1 ) p2 l 
1 1 ‘n2 1 1 \]° D,~ ” v—2 
is L — nena 1 
n' Nt t, \(a+ 2)! N! 1 K(N){S t,2 neh +is 
V—2 n=2 
25°87 11 1 \? me 
Dy 3 - | ny car 2n+-1< 3n Donec Ry/Dy pour N assez grand sera inférieur a 
. Nh l . , ’ , - 
a disons 24 P? et l’on peut donc déterminer des fonctions ) 
-— — telles que Ry/Dy soit aussi voisin de 0 que lon veut 
? "42 > = 72 en choisissant P > 0 assez petit. 
) to. n = n 
Dy> 3 > nv NP np>l 
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RAPID COMPRESSION HEATING OF A PLASMA IN THE LINEAR Z-PINCH 


H. A.B. Boptn, A. A. NEWTON AND N. J, PEACOCK 


A.W. R.E., ALDERMASTON, BERKSHIRE, ENGLAND 


The electrical characteristics of discharges with high and low initial rates of current rise are briefly described and 


the electrical measurements are used to determine the work done on the plasma by rapid compression. The particle 


and the magnetic energy in the discharge are discussed for both high and low initial rates of current rise (d//d? (0)) and 
for discharges with and without an axial magnetic field (/7_). It is shown that the gas is heated much more efficiently 


both with and without H/, at high d//d¢ (0). For d//dt (0) ~ 1.5 + 10! amps/sec (H, 


0) 15°, of the initial bank 


energy goes into heating the gas and 10°, into the magnetic field at the first maximum contraction; for d//d¢ ~ 10"! 


amps/sec (H, 


0) 11°,, goes into the field and only 6.5 


into the gas at this time. The energy input to the gas is 


found to become limited well before peak current, and the initial implosion is not followed by much of the 


subsequent isentropic compression heating which might be expected, so that the temperature never exceeds about 
300 eV; it follows that most of the observed neutron yield (~ 106) is not thermonuclear. Data obtained from the 


electrical waveforms and from high-speed photographs are presented and analysed, and it is shown that large currents 


build up outside the central plasma column after the first maximum contraction. Further evidence suggests that 
these currents flow in wall material boiled off by radiation from the central column. These wall effects account for the 


absence of further isentropic compression heating and also contaminate the central plasma column; they are such 


as to prevent any significant gain in heating by increasing the power input to the discharge. It is concluded that 


similar effects probably account for the failure to achieve a thermonuclear temperature in many high power pinch 


experiments. 


1. Introduction 


The rapid heating of a plasma is desirable for a number 
of reasons. For example, if a cold gas is heated up it 
passes through a stage when it is only partially ionized 
during which time it loses energy very rapidly by transi- 
tions involving bound states*. It is desirable to minimise 
the duration of this phase by making the power input large 
compared with the radiated power. Also, the con- 
tainment problem is simpler at high temperatures. 
Heating by rapid magnetic compression (see, for 
example, Andrianov et al. [1], Bickerton [2], Fitch [3]) 
is one of the most promising methods of obtaining a hot 
plasma. The simplest geometry for this approach is 
the linear z-pinch. 

In earlier experiments on fast linear z-pinches carried 
out at A.W.R.E. the neutron emission was studied at 
high intial rates of current rise (d//d¢(O)); this work 
showed that even from high-power discharges the 
neutron emission was not thermonuclear (Bodin et 
al. [4]). There followed a more detailed investigation 
of such discharges, with both high and low values of 
dl/dt(O), (hereafter described as “fast” and “slow’’), 
both with and without an axial magnetic field (/7,). 
The main qualitative features of fast and slow discharges 
and a few of the measurements on rapid compression 
heating, were reported by Bodin et al. [5]. 

In this paper, which is divided into two main parts, 
the experiments on rapid compression heating are 
described in detail. The first part discusses the energy 
put into the gas and the magnetic field up to the time 
of the first maximum contraction. It is then shown that 
in high power discharges the energy input to the gas by 


* 


This has also been pointed out recently by J. L. Tuck in Phys. 
Rev. Letters, 3, 313 (1959). 
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compression never significantly exceeds its value at that 
time. The second part of the paper presents and analyses 
results which show how the heating is limited by wall 
effects. The consequences of these effects are discussed. 


2. Apparatus and method of measurement 

The experiments were carried out using the Maggi 
condenser bank, described by Fitch and McCormick [6]. 
The capacity used was 100 “wf and most of the measure- 
ments were made at 20 kV. The total inductance, L,, 
external to the gas discharge was 5 = 10°* henries for 
fast discharges and 200 = 10°* henries for slow dis- 
charges, giving d//d¢(0) in the gas of 1.5 = 10! and 
1011 amps/sec respectively. A clear quartz discharge 
tube, 34 cms long, 15 cms bore and 8 mm wall thickness 
was used, with plane copper electrodes. Weak preioniza- 
tion was always used in the experiments described. 
Streak photographs were taken and the voltage (1°), 
current (/) and d//dt were recorded. 


3. The electrical characteristics of discharges with 
high and low initial rates of current rise 

Examples of the voltage and current waveforms are 
shown in figure 1 for fast (A) and slow (B) discharges. 
A (above) shows a typical current oscillogram for a 
discharge in which the bank inductance is comparable 
to the initial tube inductance, but much less than the 
final tube inductance. In this case the current rises very 
rapidly to a high value (360 kA in 0.6 sec), and then falls 
to a lower value (208 kA) at the first maximum contrac- 
tion (first pinch) of the plasma. This form of the variation 
of current with time is desirable because (see Fitch [3]) 
at the first maximum contraction more energy can go 
into the gas than into the magnetic field associated with 
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Fig. 1 
A. dl/dt (0) ~ 1.5. 


10" amp/sec, I’g = 20 kV, py = 320 Dy 


the discharge. (See section 4.) The voltage is not 
sensitive to changes in the tube inductance when the 
bank inductance is small, and the oscillogram in this 
case (A) dces not show any of the rapid fluctuations at 
the time of the pinch which are characteristic of a 
high-inductance bank (see B). The high frequency 
oscillations seen on the voltage oscillogram, and to a 
lesser extent on the current, are due to cable reflections 
in the bank (Fitch and McCormick [6]). The current 
in the slow discharge (B) is predominantly determined 
by the large series inductance and shows only slight 
evidence for pinching. In this case the current at the 
first maximum contraction is 190 kA. 

The application of 7, results in several apparently 
stable oscillations of the plasma, as found by Bezba- 
tchenko et al. [7] The oscillations were seen on streak 
photographs and also on the current waveform for a low 
bank inductance (Fig. 2). The peak current is larger 
for H, 0, contrary to predictions based on simple 


inductance considerations. (See also section 6.) 








usec 
Fig. 2 Tracings of typical current waveforms with and without 
H,. py = 350 w Dg, Vy = 20 kV, dl/dt (0) ~ 1.5 - 101? amps/sec. 
A. H, =0 B. H, 1480 oe 


4. The particle and magnetic energy in discharges 
with high and low initial rates of current rise 


THE METHOD OF CALCULATING THE DISCHARGE ENERGY 
The total amount of electrical energy put into the 
discharge tube up to time, /, is 


> 


\V lds (1) 


Tracings of typical current and voltage 


RAPID COMPRESSION HEATING OF PLASMA 


272kA 


CURRENT 





20kV 
VOLTAGE 





A. 4 4 4 


msec 


waveforms at high and low d//dt (0), H 0 


B. dl/dt (0) ~ 10" amps/sec, I’, 350 4 Dy 


where | is the voltage across the tube and / the current. 
Using equation (1) it can be shown* that at time / 


(\Vldt=3LI+3\ i (dL/ds)dt (2) 


where L tube inductance. 

By definition | LJ? is the energy in the azimuthal 
magnetic field, and therefore the last term must represent 
the total mechanical work done on the system as a result 
of the motion (due to the magnetic forces) of current 
carrying conducting material, i.e. as a result of changes of 
inductance. The total mechanical work 
discharge is 


done on the 


V =} \ [2(dL/dt)\dt 3 

Equation (3) is independent of the geometry since Hf 
includes the work done due to the motion of all the 
currents in the tube. Where all the currents flow in 
a central plasma column all of the mechanical work is 
done on this column. 


For discharges with no axial magnetic field W goes 
into particle energy W’, and 
Y= (4 
For discharges with 17, 
Y.=WVW—wv’ 5) 


U 


where W’ is the energy 
trapped field. 

The total magnetic energy, W’,, 
for unstabilized discharges, 


required to compress the 


, in the discharge is, 


W ILE (6 


H HO 2 
The total magnetic energy in the discharge for dis 
charges with H, is 


Vv, —W W (7) 


H 


~ 


HO 


where W,. is the energy in the trapped axial field. 


42 


Resistive heating, which can be shown to be small in times ot 
interest, is neglected. 


uw 
uw 
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THE PARTICLE AND MAGNETIC ENERGY IN THE DISCHARGE 
UP TO THE TIME OF THE FIRST MAXIMUM CONTRACTION 
Discharges without axial magnetic field 

Figure 3 shows the particle and magnetic energies 
in the discharge as functions of radius up to the first 
maximum contraction for the fast and the slow discharges. 
In the fast discharges it is found that whereas WW, rises 
steadily to a maximum at the first maximum contraction, 
W,, levels off so that W/W, = 1.5 at this time, i.e. more 
energy is put into the gas than into the magnetic field. At 
the first maximum contraction about 15°, of the total 
energy initially in the bank (20 kJ) has gone into the gas, 
and about 10°, into the magnetic field. 

In the slow discharge W,, is always greater than W, 
and both increase steadily as the compression proceeds. 
At the first maximum contraction 6.5°, of the energy 
in the bank (20 kJ) has gone into the gas and 12°, into 
the magnetic field. In the slow discharge most of the 
work is done at small radii, in contrast to the fast 
discharge in which a substantial amount of energy is 
put into the gas during the early part of the compres- 
sion, when the current (Fig. 1) is high. 
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Fig.3 Total particle energy (W,) and magnetic energy (W’,) in 
the discharge as functions of radius at high and low d//dt (0) 
without axial magnetic field. 


There is a gain by a factor 7 in mean power input and 
by a factor 2.3 in W, for the fast discharge. An increase 
in W, by a factor of only 2.3 for a factor 15 increase in 
the inital d//d¢ is not unexpected when d//dt (0) is changed 
by varying the bank inductance, using a relatively low 
voltage. It arises because the very high initial d//d¢ in the 
low inductance case is only sustained for a small fraction 
of the compression time near the start while the tube 
inductance remains comparable to the bank inductance. 

At the first maximum contraction the plasma radius 
in the fast discharge is larger by a factor 2.1 than that 
for the slow discharge; this is a consequence of the fact 
that more energy is put into the gas in the former case 
during the earlier part of the compression, and the 
plasma pressure limits the contraction at a larger radius. 
Another example of this is seen inthe streak photographs in 
fig. 9, for H, = 0, which show that the plasma radius at the 
first maximum contraction increases with the voltage. 
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Fig.4 Total particle energy (W,) and magnetic energy (W’,,) in 
the discharge as functions of radius at low and high d//d¢ (0) with 
axial magnetic field, H,. 


Discharge with axial magnetic field 

Moderate axial magnetic fields were applied to reduce 
the fastest growing modes and improved stability was 
found (Bodin et all. [5]) for fields ~ 1000 oe. Fig. 4 
shows the particle and magnetic energy (including that 
in the trapped field) in the discharge as functions of 
radius up to the first maximum coniraction for the fast 
and the slow discharge with H, = 1480 oe. It is found 
that with an axial magnetic field W, is always less than 
W,, in both cases, although only by 13°, for the fast 
discharge. Comparing figures 3 and 4 it is seen that the 
application of H/, does not have a marked effect on W, 
and most of the magnetic energy is still in the Hy field. 


lon temperatures 
The values of W, and ion temperature, 7, at the first 
maximum contraction are shown as functions of initial 
pressure for the fast discharge in fig. 5. For H, = 0 the 
concept of temperature is fully valid in most cases and 
is always approximately valid, as ion thermalization will 
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Fig.5 Total particle energy (W,) in the discharge and ion tem- 
perature (7) as functions of initial pressure at high d//d¢ (0) with 
and without axial magnetic field, H,. 
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occur in ~ 0.2 “sec in the longest case given (7°, ~ 300eV 
at intitial pressure ~ 254). For H, = 1480 oe the 
trapped field will delay the thermalization unless, for 
example, there is a shock wave, and values of 7, are not 
given for this case. 7; is calculated from W, assuming 
no losses and that all the energy is in the ions, i.e. 


where V/ — number of deuterons in the tube. 

It is seen that the highest temperature which could 
have been achieved at the first maximum contraction is 
300 eV at 25 mw, the lowest pressure at which a discharge 


can be initiated with a bank voltage, I,, of 20 kV. 


THE WORK DONE ON THE GAS UP TO PEAK 
CURRENT-MAXIMUM TEMPERATURES ACHIEVED 


Figs. 6 and 7 show W,, and W,,, as functions of time 
for H,=0O and 1480 oe respectively*. In the fast 
discharge W,, rises to a maximum at the time of the 
first maximum contraction, and then falls as the plasma 
expands against the field; as the plasma contracts for a 
second time W, begins to rise again but it is found that, 
for H, = 0, it never exceeds the value reached at the 
first maximum contration. For 7, = 1480 oe (see fig. 7) 
it only exceeds the early value by 20°,, during the later 
stages. Thus, for H, 0 there is no further compres- 
sion heating after the first maximum contraction, while 
for Hl, = 1480 oe there is only a little subsequent 
heating. The observation that the plasma temperature 
never significantly exceeds that achieved at the first 
maximum contraction is general in the present experi- 
ments and is discussed in the next section. 
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Fig.6 Total particle energy (W’,) and total magnetic energy (W’,,) 
in the discharge as functions of time at low and high d//dt (0), 


H, 0. 


* W, after the first maximum contraction includes a small 
contribution due to currents which are found (see section 5) to 
grow outside the central plasma column after this time. At earlier 
imes all W, goes into the central plasma column. 


RAPID COMPRESSION HEATING OF PLASMA 


The thermonuclear yield expected from the calculated 
ion temperature in these experiments (see preceding 
section) is only 10®—104, which confirms the conclusions 
reached in the study of the neutron emission (Bodin et 
al.[4]) that most of the observed yield (~ 10°) was not 
thermonuclear. 
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Fig.7 Total particle energy (W) and total magnetic energy (W,,,) 
in the discharge as functions of time at low and high d//d? (0) 
with axial magnetic field. 


5. The heating limitation in high-power 7-pinches 


The processes limiting the heating were revealed 
when data from the streak photographs and the elec- 
trical waveforms were compared. 


THE METHOD OF ANALYSIS OF THE STREAK 
PHOTOGRAPHS AND THE ELECTRICAL MEASUREMENTS 
The inductance L (/) of the discharge can be obtained 

from | and J by the equation 


| (d/dt) (LT) (8) 


which can be written 


L(t) 0 (9) 


Resistive terms have again been neglected because, 
except in the early sheath formation stages and in the 
later stages when impurities from the walls enter the 
discharge, the current and voltage waveforms are con- 
sistent with a purely inductive impedance of the discharge. 
Furthermore, if resistance is included the main conclusion 
(that /, can substantially exceed 7, — see below) would 
be strengthened. 

The plasma radius r(/) can be found from the induct- 
ance using the relation 


L(t) = 2/log [R/r (7)] x 10°-* henries (10) 


where R radius of outer coaxial current return and 


/ = tube length, on the assumption that the current 
flows ina thin sheath. (The inductance for a thick sheath or 
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for a uniform current is larger.) The radius of the current 
sheath found in this way from the electrical measurements 
is denoted by 7,. The radius of the outer boundary of the 
central luminous column (r,) is determined directly 
from the streak photographs. Values of the current 
corresponding to 7, and r, are, respectively, /, (measured 
directly by the monitor) and /, (determined from the 
streak measurements, using equation (10) to find a value 
of L.(¢) from r, and inserting this value and the observed 
voltage into equation (9)). 7, only equals /, if all the 
current Jows on the outside of the luminous central 
column; in this case r, and 7, are also equal. 


Results obtained from analysis of the streak. and 
electrical measurements 


Fig. 8 (A, above) shows /,, /, and also the current 
given by theory as functions of time for three different 
discharges; fig. 8 (B, below) shows the corresponding 
values of r,, 7, and theory. The initial gas pressure is 
arranged so that in all three cases the mass of gas, M, 
per cm length of the tube is approximately the same, the 
voltage is 20 kV in each case so the pinch-time scaling 
factor (1“*/M/)j{ is the same for the three discharges to 
a good approximation. The theoretical curves were’ 
computed for the present experiments by Taylor and 
Nisbet [8] from the Rosenbluth [9] snowplough equa- 
tions, taking into account the effect of condenser-bank 
inductance and /7.. 

It is found that /,, /, and theory are in fairly good 
agreement in the early part of the discharge in all cases 


WITROGEN DISCHARGE 4, =0 












DEUTERIUM DISCHARGE H,~ 0 





(see fig. 8, A). The biggest discrepancy during this time, 
for deuterium with H, 0, can be explained by 
uncertainties in determining r, during the very early 
stages when the whole curve is very sensitive 
values of 7, just after the plasma leaves the walls. 

The agreement between /, and /, in the early 


to the 


stages 
is sustained for more than two successive contractions 


in the case of deuterium with //, 1480 oe. For 
deuterium with 7 0, 7, begins to exceed /, soon 
after the first maximum contraction; /, falls as the 


central column contracts for a second time (see observed 
radius time curve) but /, continues to rise. For nitrogen, 
/, becomes much greater than /, before the first maximum 
contraction and the first dip in the observed current 
occurs well before the first maximum contraction of the 
central column. When the discrepancy appears the tube 
is presenting a lower inductance to the bank than it 
would if all the current flowed in the central region. 
This lower inductance represents a larger “‘effective” 
radius for the current channel and the radius/time curves 
bear this out. For nitrogen r, is about 5 cm while 7, is 
less than 0.5 cm at the first maximum contraction, and 
most of the current must flow in the outer regions of the 
tube. The observed peak current for a nitrogen discharge 
is 800 kA, while for deuterium it is 600 kA, although the 
radius of the luminous central column is smaller for 
nitrogen. 

The evidence from the current measurements could be 
consistent either with all the current flowing at radii in 
the neighbourhood of r,, or with a current /, flowing on 
the outside of the central column together with other 


DEUTERIUM DISCHARGE H, - 1480 o¢ 
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Fig. 8 A (above) Comparison between the current/time variation measured directly (/,), and calculated from the streak measurements (/ 


and according to theory for three different discharges. 


B (below) Comparison between the radius/time variation determined from the streak photographs (r,) 


, calculated from the observed 


current (r,), and according to theory for the same three discharges. 
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RAPID COMPRESSION HEATING OF PLASMA 


H,= 1480 Oe 
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Fig. 9 Radial streak photographs at different values of bank voltage with and without /7, (p ~ 


2R) 15.4cm. r, is measured from 


currents at radii >r,. However, the agreement between /, 
and theory and r, and theory during the whole of the im- 
plosion phase in all three cases shows that the dynamical 
behaviour of the discharge is consistent with a current / 
flowing in the central column, even though, for the 
nitrogen discharge, the total current becomes much 
greater than /, before the first maximum contraction. 
\fter the first maximum contraction of the radius time 
behaviour of the central column is still largely consistent 
with a current /, flowing in it, and the most likely ex- 
planation of the discrepancy is the growth of additional 
current outside the luminous central column. 

The observation that /, can exceed /, is a general one, 
of which the cases shown are typical. The current dis- 
crepancy is observed to depend on the energy input to 
the gas, the compression ratio and the presence of high-7 
impurities in the gas. The growth of current outside the 
central region has 


et al. [1]. 


also been observed by Andrianov 


Qualitative observations from the streak photographs 

\t about the time current begins to grow outside the 
central volume an increase in the background luminosity 
is observed on streak photographs, consistent with a 
marked increase in luminosity near the discharge tube 
walls. Similar effects have been observed by Burganov 
et al. [10] and Komelkov [11]. The onset of this light 








4.- SECS 


300 


350 ww). Diameter of tubc 


} 


the outer edge of the luminous sheath. 


varies with ratio and the 
presence of impurities in the gas in the same way as the 


current discrepancy. 


energy input, compression 
In nitrogen discharges faint back- 
ground light appears during the initial implosion and 
becomes very intense at the first maximum contraction. 
In deuterium discharges at pressures > 100 yw this light 
appears at or soon after the first maximum Contraction 
and is very much brighter than the central colum, which 
becomes completely indistinguishable. Figure 9 shows 
streak photographs in deuterium at different bank volt 
ages (I) both with and without an axial magnetic field. 
The middle two pictures correspond to the curves for 
deuterium in fig. 8 and it is seen that with 7 0) the 
onset of background light and the current discrepancy 
both begin at 1.5 usec; with 7 1480 oe the light 
comes later, and /, and /, are in agreement for much 
longer. 

At initial <= 100 @ the central 
column is invisible, but background light still appears 


deuterium pressure 
very brightly at or soon after the first maximum con- 
traction (determined electrically). Even in discharges 
with 99.99°, pure deuterium the same general phenomena 
are observed on the streak photographs. The addition 
of 21/,°,, of nitrogen to a deuterium discharge radically 
alters its character, and it behaves more like a pure nitro 
gen discharge than a deuterium discharge in the respects 
described above. 
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Interpretation of the results obtained from the streak 
photographs and electrical measurements 


The observations decribed above can largely be ex- 
plained by the following qualitative theory. At some 
time during the discharge the surface of the tube wall is 
vapourised by radiation from the central plasma, and 
current (in the same direction as the central current) 
builds up in the gas produced in this way, which becomes 
highly luminous due to the presence of oxygen and 
silicon. This model is consistent with the observation 
that wall effects are more pronounced and occur sooner 
if the energy input is raised, or if the content of impurities 
is increased or if the compression ratio becomes larger. 

A discharge in wall material, once initiated, will grow 
rapidly and the radiation need only provide a trigger. 
As the temperature of the wall surface depends on the 
radiation flux, the total intensity of radiation required 
can be quite small. It is easy to account for large radiation 
fluxes in heavy gases, but in deuterium the bremsstrah- 
lung can be shown to be negligible and the radiation 
could either come from impurities drawn in from the 
walls in the initial stages or from partially ionized 
deuterium, or possibly from some kind of turbulence. 
It has been shown (see, for example, McWhirter [12] and 
Kogan [13]) that the radiation from a deuterium plasma 
can be greatly increased by the addition of very small 
amounts of impurity. The wavelength of the radiation 
must be short, as wall effects have been observed when no 
visible light is detected from the central region; further- 
more, quartz is transparent to wavelengths > 2000 A. 

Wall effects can also be initiated by plasma which 
escapes from the central region and hits the walls. In 
some slower unstabilized discharges in light gases there 
is evidence for plasma jetting out and causing a highly 
localised increase in luminosity at the walls (see, for 
example, fig. 9 for 1’, = 15 kV, H, = 0). For d//d 
1.5 x 10! amps per sec the onset of wall light is usually 
symmetrical (see fig.9 for ¥, = 20 and 25 kV and 
H, = 0 and 1480 oe), and occurs in times rather short 
to be explained by plasma hitting the walls; in these 
conditions the radiation mechanism is much more likely. 

Further evidence for a discharge in wall material is 
the appearance of a second, highly luminous, sheath 
which contracts slowly in some conditions at low power; 
this sheath begins to form when the background of light 
from the vicinity of the walls appears. In other conditions 
the luminous region appearing near the walls is seen to 
move in slowly without forming a distinct sheath. The 
electron temperature in the discharge in vapourised wall 
material will be low and resistive effects important, lead- 
ing to a thick current sheath. 

These wall effects can explain many anomalies; for 
example, the appearance of the first “dip” in the current 
waveform well before the first maximum contraction in 
nitrogen discharges; the observation that the peak cur- 
rent for discharges in nitrogen is 30°, greater than for 
those in deuterium; the decrease in peak current for 
deuterium discharges when H, is applied. 


The influence oj wall effects on the compression heating 


The production of conducting material at the walls 
and the subsequent growth of currents in this material 
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screens off the inductive electric field from the central 
column and the rate of current rise in this column falls. 
The current in the central plasma is therefore limited and 
there is little further compression heating. At later times 
the central plasma becomes heavily contaminated, as 
observed by Andrianov et al. [1], Komelkov [11] and 
others. 

If attempts are made to increase the energy input to 
the gas by, for example, increasing the voltage, then 
wall effects appear earlier, and the heating becomes 
limited sooner, with negligible gain. Fig.9 shows 
evidence for this from streak photographs at I”, = 15, 
20 and 25 kV; it is seen that the background light comes 
sooner at higher voltages. Increasing the axial magnetic 
field or reducing the voltage both reduce the energy 
density in the plasma at the time of the first maximum 
contraction and are observed to delay the onset of wall 
effects. In fact it is found that over a wide range of con- 
ditions there is a limit to maximum energy input to the 
gas. There is evidence for this in fig. 6 and 7 which show 
that the maximum value of WW, is almost the same with 
or without H7,. With H, = 0 more energy is put into 
the gas at the first maximum contraction, and there is 
no further heating. With H/,, the energy input to the 
gas is less at the first maximum contraction, wall effects 
are delayed, and there is some further heating —suf- 
ficient to increase W, to about the value achieved with 
H, = 0. In general the results suggest that, in the 
present experiments, as soon as the energy density in 
the central plasma exceeds a certain value, wall effects 
appear and prevent further energy input to the central 
column. In some experiments using very high voltages 
(120 kV) Andrianov et al. [14] found that large cur- 
rents flowed at the walls right from the start of the dis- 
charge. 

6. Conclusion 

In low-power z-pinches insufficient energy is put into 
the discharge to achieve a thermonuclear temperature, 
even if the initial implosion is followed by further com- 
pression heating; in this case such compression takes 
so long that the heating is limited by instabilities, im- 
purities from the walls and other cooling processes. In 
very fast pinches wall effects limit the heating at a much 
earlier stage and so in no case is a thermonuclear plasma 
produced. However, the ratio of particle to magnetic 
energy in the high-power case, before the heating be- 
comes limited, is very favourable for further compres- 
sion if wall effects can be eliminated. It seems probable 
that the type of effects observed in the present experiment 
are responsible for the failure to achieve a thermonuclear 
temperature in many of the high-power, fast-pinch ex- 
periments so far carried out. 
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NUCLEAR FUSION 1, 62-63 (1960) 


ENERGY DISTRIBUTION OF NEUTRONS PRODUCED 
BY A THERMONUCLEAR REACTION* 


W.R. Faust ANp E. G. Harris *?* 


U.S. NAvAL RESEARCH LABORATORY, WASHINGTON D. C., U.S.A 


The energy distribution of neutrons in a thermonuclear reaction is calculated. It is assumed that the reacting 
mixture of gases is in thermodynamic equilibrium and that the differential cross section for neutron production is 


isotropic in the center of mass system. 


The energy distribution of neutrons produced by a 
thermonuclear reaction has recently become of interest 
because of the possibility of producing such reactions [1] 
in the laboratory. Some calculations [2] have been made 
previously but were completely different in approach 
from the methods presented herein. 

In order to calculate the neutron energy distribution, 
it is assumed that the reacting mixture of gases is in 
thermodynamic equilibrium at a temperature 7’ and that 
the equilibrium is not disturbed appreciably by the 
nuclear reactions which are taking place. It is also as- 
sumed that the differential cross section for neutron pro- 
duction is isotropic in the center of mass coordinate 
system of two interacting particles. 

It follows [3] from the assumptions made above that 
the number of pairs of ions with velocities in the range 
de,, de, is 
VV (My My)’ !s 


Nye (22kT)8 *P | me 


aap mc2)| dey de, (1) 

Boltzmann’s constant 

7 = temperature 

C,, C, represents the velocities of ions of mass 
m, and m, respectively 

> Me numbers of ions per unit volume of 
mass #,, Ms respectively 


where k 


{lif m, + m, 
Y12 17 jf 
\2 If M7, Ms 

The differential rate for neutron emission into the 
solid angle d{!2 = 22 sin@ d@ is obtained from (1) by 
multiplying by ga(g) where g c, —c¢,| and a(g) is 
the appropriate reaction cross section, namely 
Vy Vq (My My)"! note 
Yio (22RKT)PS 


2X 2 al 4 ] dQ ) 
-- PET V1 + Moly I ite ee (2) 


The energy distribution of neutrons observed in a co- 
ordinate system at rest relative to the gas may be obtained 
from (2) by transformation of coordinates and subse- 
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quent integration over variables extraneous to the prob- 
lem. The transformation involves a change of variables 
from c,, Cc, and cosf to a new set g,G and £ the 
neutron energy. It is convenient to perform the trans- 
formation in two separate steps; the first [4] 


. My 
c¢=G = 
My 
mW 
, 1 , 
c= G (3) 
. My 
cos 4 cos f’ 


where #%/ = ™, + mz 


The second transformation is 


gs=—s 
G’ G 
2 
cos 4’ m E (sy —" (4) 
2 G uv 


In the last equation of (4) & is the neutron kinetic 
energy measured in the laboratory coordinate system 
while # is the neutron velocity in the center of mass co- 
ordinate system and ™ is the neutron mass. The Jacobian 
transforming the original set of differentials to the new 


set dg, dG, df: is 


d(C, , Cy, cos 9) 
0(g,G,£) 


d(c,,c,cos9) d(g’, G’, cos 4’) 
0(g’,G’,cos#’) d(g,G,F) 


(5) 


Chapman and Cowling [5] have evaluated the first 
determinant on the right hand side of (5) and show it to 
bequal to unity, while the second determinant is 


d (g’, G’,cos 0’) 1 


d (g,G,F) (6) 


muG 
where G G|. 

Limits of domaine for g, G and F are obtained from 
the original limits of the old variables subject to the 
restrictions implied by the transformations, thus it is 
found that 
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NEUTRON ENERGIES IN THERMONUCLEAR REACTIONS 


Since the limits of G are determined by F and 4a, the 
G integration must be performed first. By utilizing the 
Jacobians found above, the new differential distribution 
can be written, 


Vy Vg (Wy Mg)! , 
sya (2UkT)o° 4) 


f 1 > , WyMe .\\ dt: 
exp } rk (m0 G2 5; *)( 4G 48 Fa (8) 


oO 
~~ = 


Integration of (8) over all directions of G and g yields 


Vy Yq (Wy Mg) '2 f 1 - Mm, My ,\) 
= 7\q £5(Z) exp | tte, G2 4 = o2 
Viet (AT)? * /*P\— 2a7T \” my © Ih 
G2dG 2 dedE 
(9) 
muG 


Next integration with respect to G and g in that order 
vields the desired energy distribution of neutrons, namely 
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For convenience the integration variables will be 
changed to £, (m, M/2 my) g , and with the substitu- 
tion 
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F:,)| where x \ 1 |, equation (10) 
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\n excellent approximation to the actual form of the 
neutron energy distribution function (11) may be ob- 
tained by noting that E ~ x @ while E, = AT is very 
much less than /. Under these condition equations (11) 
can be approximated as 


VV; 2 md F: f =] |?) 

> O-trFE 
U(AT)* 49) 2mm, mx O~*P | mi T|' E Ia(Q+E£,) | 
\ dl ,o (£:,) exp ( E. kT) (12) 


It follows from equation (12) that the width of the 
neutron energy distribution at e~! of its maximum value 
is (to terms of first order) 
wre 
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Ol: QO (13) 
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ABSTRACTS IN ENGLISH 


Plasma oscillations, (I), I. B. BERNSTEIN (Project Matterhorn, Princeton University, N.J.,U.S.A.) AND S.K. TREHAN 
(New Delhi, India) Nuclear Fusion 1, 3 — 41 (1960) 


This article is the first part of a review of the theory of wave phenomena in plasmas. The basic kinetic theory is developed, taking 
into account coulomb collisions. From this, magnetohydrodynamic equations are derived. It is demonstrated that it is legitimate to 
deal with a closed subset of these equations either in the limit of collision-dominated phenomena, or in the limit where the effective 
phase velocities of the phenomena of interest are much greater than thermal speeds. In case of collision-dominated plasmas the 
theory of transport coefficients is discussed. 

These equations are then applied to an extensive treatment of small-amplitude wave phenomena in plasmas. A discussion of the 
dissipative effects on hydromagnetic waves is given. Hydromagnetic waves are also considered from the Chew, Goldberger and 
Low theory. Longitudinal and transverse oscillations in current-carrying plasmas are also discussed. 

Oscillations of a cylindrical plasma are considered and the phenomenon of ion cyclotron resonance is discussed. The possibility of 
radiation by plasma oscillations by a uniform sphere is exhibited. Some general results on the stability of longitudinal electron 
oscillations in non-uniform plasmas are given. A brief treatment of large amplitude electron oscillations is given and the breaking 
of these oscillations as a dissipative mechanism for the organized plasma motion is discussed. 


Part II of this paper, to appear in a future issue of this journal, will be devoted to the discussion of plasma oscillations directly 
from the kinetic theory. 


Rayleigh-Taylor instabilities of magnetically accelerated plasma, T. S. GREEN AND G. B. F. Nisierr ( Asomic 
Weapons Research Establishment, Aldermaston, Berks., England) Nuclear Fusion 1, 42 — 46 (1960) 


High speed photography has been used to study the compression of a deuterium plasma by an axial magnetic field generated by a 
wide single-turn coil. In this experiment the inductance of the coil is small compared with that of the external condenser bank 
so that the plasma implodes at a constant acceleration of 5 - 10! oms/sec?. Flutes which develop on the outer surface of the plasma 
are interpreted as magnetohydrodynamic Rayleigh-Taylor instabilities produced by the inward accelaration of the interface sepa- 
rating magnetic field and plasma. The observed growth rates are in agreement with simple theory. It is shown that the ratio of 
the diffusivities of magnetic field and momentum determines whether the plasma resistivity or its viscosity is the dominant mecha- 
nism for damping the short wavelength instabilities. At the temperatures and densities of this experiment the resistivity is primarily 
responsible for the damping and the wavelength of “maximum instability” calculated on this basis agrees with the observed value. 


A necessary condition for hydromagnetic stability of plasma with axial symmetry, C. MERCIER (Group de 
recherches de 1’ Association Eurotom-C.E..A., Fontenay-aux-Roses, Seine, France) Nuclear Fusion 1,47 — 53 (1960) 


In this work a necessary condition is derived for the stability of a plasma with axial symmetry. This condition corresponds to 
displacements which are localized in the neighborhood of those surfaces of constant pressure on which the lines of force are closed. 
This criterion reduces to Suydam’s criterion in the case of cylindrical symmetry. 


Rapid compression heating of a plasma in the linear, z-pinch H. A. B. Bopin, A. A. NEwTON AND N. J. PEACOCK 
(Atomic Weapons Research Establishment, Aldermaston, Berks., England) Nuclear Fusion 1, 54 — 61 (1960) 


The electrical characteristics of discharges with high and low initial rates of current rise are briefly described and the electrica 

measurements are used to determine the work done on the plasma by rapid compression. The particle and the magnetic energy in 

the discharge are discussed for both high and low initial rates of current rise (d//d¢(0)) and for discharges with and without an 
axial magnetic field (H,). It is shown that the gas is heated much more efficiently both with and without H, at high d//d7(0). 
For dI/dt(0) ~1.5 - 101" amps/sec (H, = 0) 15%, of the initial condenser bank energy goes into heating the gas and 10% into the 
magnetic field at the first maximum contraction; for d//dt ~ 10" amps/sec (H, = 0) 11°, goes into field and only 6.5% into the 
gas at this time. The energy input to the gas is found to become limited well before peak current, and the initial implosion is not 
followed by much of the subsequent isentropic compression heating which might be expected, so that the temperature never 
exceeds about 300 eV; it follows that most of the observed neutron yield (~10*) is not thermonuclear. Data obtained from the 
electrical waveforms and from high-speed photographs are presented and analysed, and it is shown that large currents build up 
outside the central plasma column after the first maximum contraction. Further evidence suggests that these currents flow in wall 
material boiled off by radiation from the central column. These wall effects account for the absence of further isentropic compression 
heating and also contaminate the central plasma column; they are such as to prevent any significant gain in heating by increasing 
the power input to the discharge. It is concluded that similar effects probably account for the failure to achieve a thermonuclear 
temperature in many high power pinch experiments. 


Energy distribution of neutrons produced by a thermonuclear reaction, W.R. Faust anp E. G. Harris 
(U.S. Naval Research Laboratory, Washington, D.C., U.S.A.) Nuclear Fusion 1,62 — 63 (1960) 


The energy distribution of neutrons in a thermonuclear reaction is calculated. It is assumed that the reacting mixture of gases is in 
thermodynamic equilibrium and that the differential cross-section for neutron production is isotropic in the center of mass system. 
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RESUMES EN FRANCAIS 
Oscillations du plasma (1), |. B. BERNsTEIN Er S. K. TREHAN (Project Matterhorn, Princeton University, Princeton, 
N.]., Etats-Unis d’ Amerique) Fusion Nucléaire 1, 3 41 (1960) 

Cet article est la premiére partie d’une étude d’ensemble portant sur la théorie des phénoménes ondulatoires dans les plasma. Les 
auteurs y développent la théorie cinétique fondamentale, en tenant compte des chocs coulombiens ect parviennent ainsi a établir 
des équations magnétohydrodynamiques. Ils démontrent que l’emploi d’un sous-ensemble complet de ces équations constitue une 
approximation légitime pour l’étude de phénoménes dans lesquels soit les chocs prédominent, soit les vitesses effectives de 
phases sont trés supérieures aux vitesses d’agitation thermique. Ils cxaminent la théorie des coefficients de transport dans le cas 
de plasmas ot les chocs prédominent. 
Ces équations sont ensuite appliquées a une étude approfondie des phénoménes ondulatoires de faible amplitude dans les plasmas. 
Les effects de dissipation sur les ondes hydromagnétiques font objet d’un examen critique. Les ondes hydromagnétiques sont 
également étudiées du point de vue de la théorie de Chew, Goldberger et Low. On examine aussi les oscillations longitudinales 
et transversales dans les plasma traversés par un courant. 
Les auteurs étudient les oscillations d’un plasma cylindrique ainsi que le phénoméne de résonance cyclotronique des ions. Ils 
montrent que les oscillations d’une sphére uniforme du plasma peuvent produire un rayonnement. Ils donnent quelques résultats 
généraux concernant la stabilité des oscillations longitudinales des electrons dans des plasmas non uniformes. 
L’article renferme une bréve étude des oscillations de grande amplitude des electrons et de l’interruption de ces oscillations en 
tant que mécanisme de dissipation pour le mouvement du plasma organisé. 
La deuxiéme partie de l’article, qui paraitra ultericurement dans un numéro de cette revue, sera consacrée a |’étude des oscillations 
du plasma, directement fondée sur la théorie cinétique. 


Instabilités de Rayleigh-Taylor dans un plasma soumis 4 une acceleration magnétique, T.S. GREEN ET 
G. B. E. Nisterr (Atomic Weapons Research Establishment, Aldermaston, Berkshire, Royaume-U ni) 

Fusion Nucléaire 1, 42 — 46 (1960) 

Les auteurs ont cu recours 4 la photographie ultra-rapide pour étudier la compression d’un plasma de deutérium par un champ 
magnétique axial créé par un enroulement large a spire unique. Dans cette expérience, inductance de la bobine est faible, comparée a 
celle des condensateurs extéricurs, de telle sorte que le plasma «implose» a une accélération constante de 5 10!" cm/sec”. De l’avis 
des auteurs, les stries qui apparaissent a la surface extérieure du plasma correspondent aux instabilités magnétohydrodynamiques de 
Rayleigh-Taylor, produites par l’accélération centripete de linterface qui sépare le champ magnétique ct le plasma. La vitesse de 
croissance observée concorde avec la théorie. Les auteurs démontrent que, le rapport entre la diffusivité du champ magnétique et celle 

du moment determine si c’est la résistivité du plasma ou sa viscosité qui contribucnt surtout 4 amortir les instabilités de courte 
longeur d’onde. Dans les conditions de température ct de densité de l’expérience, l’'amortissement est dd surtout a la résistivité 

ct la longeur d’onde «d’instabilité maximum», calculée a partir de cette hypothése concorde avec la valeur effectivement observée. 


Un critére nécessaire de stabilité hydromagnetique pour un plasma en symétrie de révolution, CLAup 
MERCIER (Groupe de recherches de P association F:uratom-C. f:,.A., Fontenay-aux-Roses, Seine, France) 
Fusion Nucléaire 1,47 — 53 (1960) 
Pour ce travail, on établit une condition nécessaire de stabilité pour un plasma en symétrie de révolution. Elle correspond a des 
déplacements localisés au voisinage des surfaces a pression constante sur lesquelles les lignes de force magnétiques sont fermécs. 
Ce critére se réduit 4 celui de Suydam en symétric cylindrique. 


Chauffage du plasma par compression rapide dans la striction linéaire avec champ magnétique axial, 
H. A. B. Bopin, A. A. NEwron Et N. J. Peacock (Atomic Weapons Research Establishment, Aldermaston, Berk- 
shire, Royaume-Uni ) Fusion Nucléaire 1, 54 — 61 (1960) 

Les auteurs décrivent bri¢vement les caractéristiques électriques des décharges avec accroissement initial du courant rapide ou lent, 
les mesures électriques suivant a determiner l’effet de la compression rapide sur le plasma. Ils étudient l’énergie corpusculaire et 
lénergic magnétique inhérentes a la décharge pour des accroissements de courant initiaux (d//d/(0)) rapides et lents, et pour des 
décharges avec ct sans champ magnétique axial (#7,). Ils montrent que le chauffage du gaz est beaucoup plus efficace lorsque la 


valeur de d//d#(0) est élevée, qu’il y ait ou non un champ magnétique axial. Pour d//d#(0) ~ 1,5 & 10" A/s (HH, = 0), 15°,, de 
énergie initiale fournie par les condensateurs servent a4 chauffer le gaz et 10°(, sont communiqu¢s au champ magnétique au moment 
de la premiére contraction maximum; pour d//¢ ~ 101! A/s (#7, = 0),11°,, sont absorbés par le champ magnétique et 6,5°,, seule- 


ment par le gaz, au moment considéré. On constate que la quantité d’énergie fournie au gaz se trouve limitée bien avant que le 
courant n’ait attcint sa valeur maximum, ct que |’échauffement par compression isentropique qui suit l’implosion initiale est beaucoup 
plus faible qu’on pourrait le prévoir, de sorte que la température n’est jamais supéricure 4 environ 300 cV; il s’ensuit que la plus 
grande partie du nombre observé de neutrons émis (~ 10*) n’est pas d’origine thermonucléaire. Les auteurs exposent et analysent 
les données obtenues par I’étude des formes d’ondes électriques et des images photographiques prises 4 grande vitesse; ils montrent 
que des courants importants se forment a l’extéricur de la colonne centrale de plasma aprés la premiére contraction maximum. D’autres 
données semblent indiquer que ces courants passent dans la mati¢re constitutive des parois chassée sous l’action du rayonnement 
provenant de la colonne centrale. Ces effets de paroi empéchent le chauffage par compression isentropique de se poursuivre et 
contaiminent également la colonne centrale de plasma; leur nature est telle qu’ils rendent impossible tout accroissement significatif 
de la température par une augmentation de la quantité d’énergie communiquée a la décharge. 
Les auteurs concluent que c’est probablement a cause d’effets analogues que l’on n’a pas pu obtenir de températures thermonucleé- 
aires, au cours de nombreuses expéricnces de striction du plasma mettant en jeu de grandes énergics. 
Distribution énergétique des neutrons produits au cours d’une reaction thermonucléaire, W. R. FAusr £1 
E. C. Harris (U.S. Naval Research Laboratory, Washington, D.C., Etats-Unis d’ Amérique) 
Fusion Nucléaire 1, 62 — 63 (1960) 


Les auteurs de la monographie font le calcul de la distribution énergétique des neutrons au cours d’une reaction thermonucléaire. 
Ils partent de ’hypothése que le mélange gazeux réagissant se trouve en état d’équilibre thermodynamique et que la section efficacc 
différentielle relative 4 la production de neutrons est isotrope dans le syst¢me du centre de gravité. 
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AHHOTAIIMUM TIO-PYCCKU 


Ocunanaumu nua3zmpr (I), UW. B. Bepnurtreitu u C. K. Tpexan (MarrepxopucKkui npoexr, IIpuxncTrouckui yHu- 
sepcureT, IIpuncron, How Jacepcu, CLIA) fAlnepHpmit cunte3 1, 3-41 (1960) 


JlaHHaaA cTaTbA ABNAeTCA NepBoit uacThbio OG630pa TeOpum BONHOBBIX ABIeHMit B MnasMax. Pa3spa6oTaHa oneann anna 
KMHeTM4¥eCKaH TeOpMA C YYETOM KYJIOHOBCKMX CTONKHOBeHMM. U3 9TOM TEOPMM BbLIBOAATCA MarHuTOrMapOAMHaMn socked 
ypaBHeHua. TloKa3aHo, 4TO MOXKHO OOpalaTbcA C 3AMKHYTOM NOACKCTeMOM 9TMX ypPaBHeHMM MIM B npezenax —e 
JOMMHMpYIOUMX CTONKHOBeHMM MIM B Cny4Yanx, KOrZa scpceKTMBHbIe Ga3z0Bpie CKOPOCTM MHTepecyIoulMx aBneH? 
3HAaYMTeNbHO 6ONbUIe TemnoBbIX cKOpocTeii. IIpH pacCMOTpeHMM JOMMHMpyIOUIMX MNa3MeHHbIX CTOJKHOBeEHH O6CYK- 
waeTcA TeEOPMA TpaHCNOPTHBIX KO9ICCMUMeHTOB. 

3aTeM 9TM ypa8SHeHMA MpMMeHAIOTCA K MHTCHCMBHOMY McCIeqOBaHMIO BONHOBbIX ABNeEHM C MaJIOM aMMIMTyAON » 
ma3Max. IIpMBogqutTcaA paccMoTpeHMe Bompoca O BIMAHMM pacceAHHA Ha YTMAPOMArHMTHbIe BOJIHbI. FuAPOMATHNTHEIC 
BONHbI pacCCMaTPMBaIOTCA Tax7Ke C TOUKM 3peHMA TeopuM Ubi, TonpaGepra u Jloy. OScyxAalwTcCA Take MPOAOJNIbHbIE 
M monepeyHbie OCLUMINAWUMM B TOKOHeECyUIMxX Mma3Max. ee 

PaccMaTpMBaloTCA OCLIMANAUMM UMNMHApMYeCcKOM Mna3sMbI M OGcyxRAaeTCA ABIEHMe MOHHO-LMKNOTPOHHOTO PesoHaHca 
TloKa3aHa BO3MOXHOCTh NONYYeHMA M3TyYeHMA B Pe3YNbTAaTe OCIIMAAAUMM MNa3MbI B OAHOpOAHOK cchepe. IIpuBoAATCA 
HeKOTOpble OGUIMe Pe3yYNbTAaTbI CTAGMABHOCTM MPOMONbHbIX YICKTPOHHbIX OCLIMIIAUMM B HCOAHOPOAHbIX M1a3Max 

KpatTKo o6cyxaeTCA TPaKTOBKa 9JICEKTPOHHbIX OCLUMIIAUMI BhICOKOM AMMIMTYAbI a TakxKe MepepbiB ITMX OCIMI- 
NAUMM KaK Me€XaHi3M pacceAHMA ANA OpraHM30BaHHOTO LWBMKeHMA [a3Mbl. 

YactTb II 9Tout cTaTbM, HAaMe4YeHHaA K ONyONMKOBaHMIO B NocnesyouleM HOMepe AaHHOTO RypHata, GyAeT NocBaAUIeHAa 
paccMOTpeHMlO OCLIMINAUMM Mna3MbI HeMocpesxCTBeEHHO C TOYKM 3peHMA KMHETMYeCKOM TeOpMit, 


Heycroiunsoctu Peitnu-Teiinopa B nmua3Me, ycKopaemoi B MarHuTHOM nose, T. C. [pun u Ix. B. ©. HuOnerr 
UccaedosaTervcKkuti weXTP NO aTOMHOMY Opys9cur;, OnrdepmacTon, Bepxwup, Ax22UAh) 
fAlyzepHprii cuntes 1, 42-46 (1960) 


TIpu M3y4eHMM CMATHA neutTepuesBon Mmia3sMbiI C MNOMOLUIbI aKCHasbHOro MarHMTHOrO MOA, o6pa30BaHHOrO UWIMpOoKON 
KaTYUIKOM C OHMM BMTKOM, GbINa MCMNONb3OBaHa BbICOKOCKOpOCTHaA chbotorpaqdua B faHHOM 9KCMepMMeHTe MHAyK- 
TMBHOCTB KaTYUIKM HacTONbKO He3HAadMTeNbHa MO CPaBHeEHHIO C MHAYKTMBHOCTbIO BHeLWIHeH 6GaTapeu KOHAeCHCATOPOB, 
4TO Tma3Ma Te4YeT C MOCTOHHHbIM yCKOpeHMeM B 5X 10** cm/cexK*. 

3aBuxpeHna, o6pa3yembie Ha BHeLIHeEM TIOBCPXHOCTM MN1a3Mbi, PaCCMAaTPMBAaIOTCAH Kak MarHuTorMNAposAMNHamMurtecKan 
HeycToMumBocTs Pelinu-Telinopa, BbI3BaHHaA BHYTPCHHMM YCKOPeHMeM TIOBeEPXHOCTH, passe Aoulen MarHMUTHOe ro 71e 
Mu na3My. Ha6ntwnaemoe BO3spacTaHMNe CKOpocTei COOTBeETCTByeT MpocTon TeopHn. TloKa3bIBaeTCA, 4YTO OTHOWICHKEe anucb- 
cdy3nn MarHMTHOTO MONA M MMITy1bCa ONpesenneT, ABNHETCA NM CONPOTMBNeHKMe MWasMbI MIM ee BASKOCTb JOMHMHHK- 
PYIOULMM MeCX@aHM3MOM JNA TauwlieHNMH KOPOTKOBOJIHOBOM He yCTOHMMBOCTH IIpu remnepatypax MUM MAIOTHOCTAX 9ToroO 
SKCNePMMeHTa COMPOTMBNeCHNe ABNAeTCA OCHOBHOM TIPMYMHON rauwieHMA U ANIMHAaA BOJIHbI MaKCMMaJIbHOHU HeycTon- 
4YMBOCTM», MOACUMTAHHAA Ha B9TOM OCHOBe, COOTBETCTByeT HaGmonaeMOl BeNN4MHe. 


HeoOxoxMMbBIit KpMTepHii THAPOMATHUTHOM cTAGMALHOCTH 1A 11a3MbI C OCceBO CHMMeTpHei, K 02 Mepcbe (Mccre- 
OosaTenvcKaa zepynna Accouuauuu Esparom — KowmuccapuatT no arTromuot sxnepeuu — MonrTene-o-Po3, 
Ppanuuag). AyepHpm cunte3 1, 47-53 (1960) 


B 9Tol pa6ote ycTaHaBNMBaeTCA HEOGXOAMMOe YCNOBMe cTaGMNbHOCTM ANA MNasMbI c OCceBOH cuMMeTpuel. OHO COOT- 
BeTCTBYeT MePeCMeCLIICHMAM, NOKAaJIM3OBaHHbIM TO COCeACTBY C NOBeEPXHOCTAMM TIOCTOAHHOTO aBNeCHMA, Ha KOTOPbIX 3a- 
MbIKaIOTCA CHJIOBbI€e MAarHMTHbie JIMHMM. DTOT KPMTePpMU CBOAMTCA K KPMTePHMio Cainema ANA cny4aa WUMNMHApNYecKOoN 
CMUMMeTPpHK. 


BaicTpoe cmaTHe HarpeBaeMOi NNa3MbI B WMHeEMHOM z mnMHYe, X. A. B. Bonun, A. A. Hproron u H. &K. Iimucox 
(AccaraedosaTervckut weHTp NO aTOMHOMY OpyacuI0, Onrdepmacroun, Bepxwup, Ax2zaun) 


AizepHpm cunte3 1, 54-61 (1960) 


B foKnafe AaeTcH KpaTKOe ONMCaHMe 9NeCKTPM4eCKMX XapaKTePMCTMK Pa3PANOB C BbHICOKO/ M HM3KOi MepBOHa4aNb- 
HOM MHTCHCMBHOCTbIO POcTa TOKa M 9EKTPMYeECKNX M3MeCPeHHii, MCNONb3yYeMbIX ANA ONpexeneHuA BLINONHeEHHOM pa- 
GoTbI MO GbICTPOMy C*aTKIO Mia3MbI. PaccMatpuMBaetca Bonmpoc 06 9HeEprMm YaCTMIL M MarHMTHOTO Nona B pa3paze 
aA Cny4aes C BbICOKOM M HM3KOM NepBOHaYaNbHOM MHTEHCMBHOCTSIO pocta TOKa ( d//d(0)) mu ANA pa3spanoB mpm Hann- 
4UM M B OTCYTCTBMM aKCManbHOTO MarHMTHOrO Nona (// TloKka3biBpaetca, 4TO ropa3fz0 Gonee sdcdbekTMBHOe HarpeBa- 
HMe ra3a MpOMCXOAMT Kak NPM HaNM4nu, Tak M B OTCyTCTBHe (H,) mpm GonbUoit BenuunHe (d//dt (O 1a dijdt (0) ~ 1,5x 10" 
ammep/cek, (H,= 0)15% NepBoHavanbHol sHeprMu KOHZeHCaTOpOB MyeT Ha HarpeBaHMe raza mu 10% Ha MarHuTHoe mone 
mpi Me€pBOM MaKCMMaJIbHOM COKpaljeHMM; fA d//dt~ 10" amnep/cex (H,=- 0.) 3a ToT xe mepuogx BpemeHu 11% uaer ua 
none M 1MuIb 65% Ha ra3. Tlepefmaua sHeprmu ra3a OKa3.1BaeTCA OrpaHM4eHHO! 3amqonro no nOCTMaKeHMA MaKCMMAaJIb- 
HOTO 3Ha4e€HMA TOKa, M NepBOHAaYabHOe peskoe vwARaTMe HE COMpOBORAaNoch NOCNeAYIOUIMM HarpeBoM M309HTponnH- 
“ecKOrO CaTMA, KOTOPOTO MOXKHO Gbi0 GbI ORKMNAaTb, TAK YTO TeMMepaTypa He MpeBbana NopagzKa 300 9B. Orciona 
cnefnyeT, “TO GONLUIMHCTBO Ha6mIoMaeMOro BbIXOZa HeMtTpoHOB (~ 10°) He TepmMOoxZepHOrO mpoucxoxmpeHMA. TIpuBo- 
AATCA UM aHaNM3MpyOoTCA WaHHble, MOAYYeHHbIe MO POPMaM 9neEKTPM4eCKMX BONH UM MO BbICOKOCKOPOCTHBIM choTorpaduamM;: 
MOKa3bIBaeTCA, YTO Mocre MepBOrO MAKCMMAaNbHOTO COKPalleHMA BOSHMKAaIOT CMIbHBIC TOKM BHe WeHTpanbHoro uwIHypa 
nla3Mbl. JatbHeliume aHHbie NMOATBepxKAaloT, ATO ITM TOKM MpoTeKaloT B MAaTepMane, BLITAaPeHHOM M3 CTeHOK M3qy- 
4eHNeM WeHTpasIbHOTO WIHypa Mia3MbIl. STO BAMAHMe CTEHOK OGbACHAET OTCYTCTBMe AanbHeliuero M309HTPONM4eCKOrTO 
HarpeBa CKaTHeM M NPMBOAMT Take K 3arpA3HEHMIO LEHTPaNbHOTO Mma3MeHHOTO LWIHypa; OHO TaKOBO, 4TO MmpexoTBpa- 
ulaeT m1060€ 3aMeTHOe yBenNM4eHMe HarpeBa MpM pocTe Normoulaemoi sHeprum B pa3pane. B 3aKnioY¥eHMe roBOpMTcA, 
“TO NOAOGHbIe SGqeKTbI, BEPOATHO, ABNAIOTCA NPHYMHOM Heytau AOCTM*ReHMA TeMNepaTypLI TepMoAReEpHBIXx peakuni 
BO MHOrTMX 9KCNe€PMMeHTaX MO NMHYY BbICOKOM MOUIHOCTM. 


Pacnpeyenenue 9Heprun HeiiTpOHOB, O}paz0BaHHLIX B pe3yILTaTe TepMoayepHOK peakuuu, Y. P. Daycr nu E. 
ix. Xappuc (Mopexaan uccredosaTeravcKaa aabopatopusx CIA, Bawuneton, Oxpye Koaymbua, CUIA) 
SinzepHpm cunte3 1, 62-63 (1960) 


IIOACUNTbIBaeTCA pacmpexeneHue 9HEprMu HelTpOHOB mpi TeEPMOAREPHON peakunu, Tpexznonaraetca, 4To pearupyw- 
ulaH CMe€Cb Tra30B HaXOAUTCA B TEPMOAMHAMMYeCKOM PaBHOBeCHMH M 4TO AMcbchepeHUMaNbHOe MonNepeyHoe ceyeHne 1A 
reHepauMit HeMTPOHOB ABIACTCA M30TPOMHbIM B CMCTeMe LIGHTpa Macc. 
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RESUMENES EN ESPANOL 


Oscilaciones del plasma (I), I. B. BERNsTEIN y S. K. TREHAN (Proyecto Matterhorn, Universidad de Princeton, 
Princeton, New Jersey, los Estados Unidos de Norte America) Fusién Nuclear 1,3 — 41 (1960) 


Este articulo constituye la primera parte de un estudio sistematico de la teoria de los fendmenos de tipo ondulatorio en los plasmas. 
Los autores desarrollan la teoria cinética en que se basan, teniendo en cuenta las colisiones couldmbicas, y derivan de ella ecuaciones 
magnetohidrodinamicas. Demuestran que cabe emplear un sistema cerrado de ecuaciones, tanto hasta el limite en que predominan 
las colisiones como hasta el limite en que las velocidades de fase efectivas de los fenémenos de interés son muy superiores a las 
velocidades térmicas. Comentan la teoria de los coeficientes de transporte para el caso de plasmasen los que predominan las colisiones. 
A continuacién aplican estas ecuaciones a un tratamiento matematico muy extenso de los fenémenos ondulatorios de pequefia 
amplitud que tienen lugar en un plasma. Analizan el efecto de disipacién en las ondas hidromagnéticas, y estudian éstas también 
desde el punto de vista de la teoria de Chew, Goldberger y Low. Discuten las oscilaciones transversales y longitudinales en plasmas 
portadores de corriente. 

Seguidamente estudian las oscilaciones de un plasma cilindrico y discuten el fendmeno de resonancia ciclotronica de los iones. 
Muestran la posibilidad de que una esfera uniforme del plasma emita radiaciones por oscilacién. Dan algunos resultados de caracter 
general sobre la estabilidad de las oscilaciones longitudinales de electrones en plasma no uniformes. 

También comentan un tratamiento abreviado de las oscilaciones de gran amplitud de los electrones y la interrupcion de estas oscila- 
ciones como mecanismo de disipacién del movimiento de un plasma organizado. 

La segunda parte del trabajo apareceré en otro numero de esta revista y tratara del estudio de las oscilaciones del plasma aplicando 
directamente la teoria cinética. 


Inestabilidades de Rayleigh-Taylor en un plasma acelerado magnéticamente, T. S. GREEN y G. B. F. NIsLetT 
(Atomic Weapons Research Establishment. Aldermaston, Berkshire, Reino Unido) Fusion Nuclear 1,42 —46(1960) 


Los autores han estudiado la compresi6n de un plasma de deuterio por un campo magnético axial producido por un solenoide 
de una sola espira de diametro grande. En este experimento la inductancia del solenoide es pequefia comparada con la de bateria 
de condensadores exteriores, de manera que la implosién del plasma se produce con una aceleracién constante de 5  10!? cm/s?. 
En la superficie exterior del plasma se producen irregularidades que se consideran como inestabilidades magnetohidrodinamicas 
de Rayleigh-Taylor, producidas por la aceleracién centripeta de la superficie de separacién entre el campo magnético y el plasma. 
Se ha observado que la velocidad del proceso responde a suposiciones tedricas sencillas. En el articulo se demuestra que la relacién 
entre la difucibilidad del campo magnético y la cantidad de movimiento determina si la resistividad del plasma o su viscosidad 
es el mecanismo dominante para amortiguar las inestabilidades de corta longitud de onda. En las condiciones de temperatura y 
densidad en que se efectué el experimento, la resistividad constituye el principal factor de amortiguamiento y la longitud de onda 
de la “‘inestabilidad maxima” calculada sobre esta base coincide con el valor observado. 


Un criterio necesario de estabilidad hidromagnética para un plasma de simetria axial, CLAUDE MERCIER 
(Groupe de Recherches de Association Euraton — C. E. A. , Fontenay-aux-Roses, Seine, Francia) 
Fusién Nuclear 1, 47 — 53 (1960) 


En la memoria, el autor establece una condicién necesaria fara la estabilidad de un plasma de simetria axial. Esta condicién corres- 
ponde a desplazamientos que quedan localizados en las proximidades de aquellas superficies a presi6n constante en las que las 
lineas de fuerza son cerradas. En el caso de una simetria cilindrica, este criterio queda reducido al de Suydam. 


Calentamiento de un plasma por compresién rapida en la constriccién lineal con campo magnético 
axial, H. A. B. Bonin, A. A. Newron yN. J. Peacock (Atomic Weapons Research Establishment, Aldermaston, 
Berkshire, Reino Unido) Fusién Nuclear 1, 54 — 61 (1960) 

Los autores describen sucintamente las caracteristicas eléctricas de descargas efectuadas con indices altos y bajos de crecimiento 
inicial de la corriente, y examinan las mediciones que se utilizan para determinar los efectos que tiene sobre el plasma una com- 
presion rapida. Se estudia la energia de las particulas y la energia magnética resultantes de la descarga con respecto a indices altos 
y bajos de crecimiento inicial de la corriente (d//dt(0)) y con respecto a descargas con campo magnético axial (H,) y sin él. Para 
valores altos de di /dt(0) cl calentamiento del gas es mucho mayor, tanto con un H, como sin él. Para un d/jd¢(0) de aproximadamente 
1,5 X 10!" amp/s (campo H, = 0), el 15 por ciento de la energia inicial de Ja bateria de condensadores pasa a calentar el gas y un 
10 por ciento al campo magnético al producirse la primera contraccién maxima; para un d//dt ~ 10" amp/s (H, = 0), el 11 por 
ciento pasa al campo magnético y sdlo un 6,5 por ciento al gas. La energia absorbida por el gas queda limitada bastante antes de que 
la corriente alcance su maximo, y la implosi6n inicial no va seguida de un calentamiento por compresi6n isentrépica de la magnitud 
que podia esperarse, de manera que la temperatura no rebasa nunca los 300 eV aproximadamente; de ello se deduce que la mayor 
parte de la cantidad de neutrones emitidos (~ 10®) no es de origen termonuclear. Los autores exponen y analizan datos obtenidos 
de las formas de onda y de fotografias ultrar4pidas, indicando cémo después de la primera contraccién maxima se forman grandes 
corrientes fuera de la columna plasmatica central. Otras pruebas sugieren que esas corrientes pasan a material de las paredes des- 
prendido por la accién de las radiaciones de la columna central. A este efecto de las paredes se debe el que no se produzca un 
ulterior calentamiento por compresién isentrépica, asi como la contaminacién de la columna plasmatica central; su naturaleza es 
tal que impide todo aumento importante del calentamiento cuando se aplica mayor energia a la descarga. Los autores llegan a la 
conclusién de que es a efectos andlogos a los que probablemente se debe la imposibilidad de alcanzar una temperatura propia 
de las reacciones de fusién en muchos experimentos realizados sobre la constriccién aplicando elevada energia. 


Distribucién energética de los neutrones producidos en une reaccién termonuclear, W. R. Faust y 
E. G. Harris (U.S. Naval Research Laboratory, Washington, D.C., Los Estados Unidos de Norte America) 

Fusion Nuclear 1, 62 — 63 (1960) 

Los autores calculan la distribucién energética de los neutrones producidos en una reaccién termonuclear. Al efectuar el calculo, 


se supone que la mezcla reactiva de gases se encuentra en equilibrio termodinamico y que la secci6n eficaz diferencial correspondiente 
a la produccién de neutrones es isotrépica en el sistema del centro de gravedad. 











